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ABSTRACT 


In this thesis, the problem of designing the layout of integrated circuits is examined. The 
layout of an integrated circuit specifies the position on the chip of functional components and wires 
interconnecting the components. We use a general model under which components are represented by 
rectangles, and wires are represented by lines. This model can be applied to circuit components 
defined at any level of complexity, from a transistor to a programmable logic array (P1.A). We focus 
on the standard decomposition of the layout prablem inte a placement problem and a routing 
problem. 


We examine problems encountered in layout design from the point of view of complexity 
theory. The general layout problem under our model is shown to be NP-complete. In addition, two 
problems encountered in a restricted version of the routing problem -- channel routing -- are shown to 
be NP-complete. ‘The analysis of heuristic algorithms for NP-complete problems is discussed, and the 
analysis of onc common algorithm is presented. 


The major result presented in this dissertation is a polynomial time algorithm for a restricted 
case of the routing problem. Given one rectangular component with terminals on its boundary, and 
pairs of terminals to be connected, the algorithm will find a two-layer channel routing which 
minimizes the area of a rectangle circumscribing the component and the wire paths. Each terminal can 
appear in only one pair of terminals to be connected, and the rectangle used to determine the area 
must have its boundaries parallel to those of the component. If any of the conditions of the problem 
are removed, the algorithm is no longer guaranteed to find the optimal solution. 


Thesis Supervisor: Ronald L. Rivest 
Title: Associate Professor of Computer Science and Engineering 


Key words: VLSI (very large scale integrated) circuit layout, component placement (rectangles), 
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Chapter 1: Introduction 


The research reported in this thesis is an investigation of algorithms for the layout of integrated 

Seana Integrated circuits are formed on silicon chips by créating layers of different substances (e.g. 
metal, polysilicon) in geometric patterns on’ the chip through a variety of fabrication techniques. 
Electronic components are formed by the interaction of regions in the different layers. ‘Wires connecting 
components are simply regions betwcen two components ona layer. Laying out a circuit consists of 
determining the patterns for each layer on the chip to create the desired components and 
interconnections. For example, in nMOS/FET technology;.a designer creates a transistor by drawing a 
sain for the polysilicon layer and a region for the diffusion layer which‘ cross when the designs for the 
two layers are superimposed [Mec80]. In general, designing a layout requires knowledge of the 
interactions between layers for the technology being used and limitations of the fabrication process being 
used. The goals of the designer are to put as much circuitry as possibic in as small an area as possible, and 
to have it work correctly and as fast as possible. A good example is the layout of a microprocessor, where 
the amount of information which can be processed, the number of functions which can be performed, 
iid the spehekne prcceeiagars kuputea 2 = 
The layout problem as described above contains a huge number of variables and leaves much 

room for cleverness by the person designing the layout. It does not lend itself well to an algorithmic 
approach where a well-defined model and set of operations are employed. However, standard layouts can 
be used for each component needed in a circuit. These components need not be simple electronic 
_ components such as transistors or resistors, but may be logic gates or even higher Icvel circuit subsystems. 
Given these components and the interconnections necessary to realize the desired circuit function, the 
layout problem consists of allocating a proper-sized region of the chip for each component and 
determining the pattern of wires forming the interconncctions on each layer. This is the layout problem 


as we will mean it. We have lost the flexibility of tailoring the layout of each component to the particular 
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application, but have greatly simplified the problem. , 
In general, we will wish to place a set of components on a:pessibly multi-layered planar surface. 
The components will coatain: (ermminala, which ase points to-which wifes can connect.. The desired 
interconnections will be specified. by giving disjoint sets of.terminals, Each set of terminals, called a net, 
should be interconnected. The interconnections will te. sade by wires which definc paths betwoen 
terminals in the layers. ‘There will be constraints which: the layout must satisfy, such ax a minimum 
separation between uncaanected. wines. - These constraints.eng called design rules. | 
The major. motivation for developing algorishms se anive the jayout problem is the:complexity 
of the integrated circuits being denigeed. A chip may-new-contain tens of thousands of transistors. Hand 
layout of these integrated circuite, even. with the aid of eonsputerized. graphics, is very costly, time 
consuming, and error prope, Standardized componenis, auch as logic gate odlls, are already employed in — 
the industry to simplify circuit. design.and layout, and layout by computer‘has tea iciptennensed {Fe76] 
| {Per77]. -The tradeoffs are similar % those in. computer.:programming. | Programming .in assembly 
language gives a programmer the freedom to devise clever waysto-iake:a program run faster or require 
dogs storage. However, when meny: large, complicated progsams nced to. be: written, the savirigs. iw time 
and the, conceptual simplification gained by using a high level language and compiler are worth the loss of 
flexibility. Just.as one-might write an. often used subroutine in: amsentbly language, one might do the 
layout for. a circuit. subsystem to be used in many circuits by hand. However, ‘for. most projects, the 
savings in.time and. the added faith in the correctness of the design due to the simplified structure make 
the component approach preferable. 
__ Although working automated layout systems do exist, she problem of designing algorithms to do 
integrated circuit layout is not solved. Many of the existing algorithms plact further restrictions on the 
problem. ‘Typically, they require that all-components be the same size in one dimension. ‘Such an 


algorithm places the components in rows, forming an array. The wiring.(calted. roxting) is done-ia the 
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spaces between rows. Also, although several algorithms have been: designed and implemented, little 
analysis of the algorithms has been. performed. . Empirical evidence “is often cited in the literature. 
However, papers describing algorithms which may not find optimal solutions rarely present mathematical 
analyses of the quality of the solutions found. For example, one does not find statements of the form 
“this algorithm always finds a solution within 50% of the optimal.” ‘In this dissertation, we discuss layout 
algorithms from a complexity theory point of view. We focus on the’ performance of algorithms, both in 
terms of the quality of the solations they preduce and the running times they require. We work with 
‘subproblems of the fayout problem since they are-easier-to approach and are-commonly encountered. An 
example of such a subprobiem is the routing problem when. components have been placed in rows. 
The thesis is organized as follows. ‘In Chapter 2, we-review the techniques used by cxisting 
automatic layout systems. In Chapter 3, the model we wilt use for layout problems is presented. We:also 
describe a second model ~- the graph thepretic model — which has proven: very useful in characterizing 
the area required by various intercompection patterns. In Chapter4,:we discuss the complexity of a 
mumber of versions and subproblems of the layout. problem. A review. of previously: known 
NP-completeness results is given. We prove the. NP-completeness of'a rectangle placement problem and 
two problems encountered in channel routing: We analyze 2 previously ‘known heuristic algorithm for 
one of the channel routing problems. In Chapter 5, a new algorithm is presented for a special case of the 
routing problem: In this problem, terneinals lie ot the boundary of one rectangutar component. Pairs of 
terminals must be interconnected. The algorithm finds a minimum .area routing for a-channel fouting 
modcl and has running time Oct’), where t is the number of terminals: By developing this algorithm, we 
have shown that there are non-trivial routing probloms which: are ‘not NP-hard. ‘Most routing problems 
are cither known to be NP-iard or are:so closely ‘rélated to NP-haid problems as 10 Be considered 
intractable without an actuat proof of NP-hardness. In Chapter6, we discuss properties of the algorithm. 


In-Chapter 7, we summarize and present open problems-and directions for fiture sancti: A review of 
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basic definitions and notation used throughout this thesis is presented in the appendix. 
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Chapter 2: A Review of Layout Attomation Techniques 
2.1 Classical Approach to Layout 


The problem of placing components on a surface and making the required interconnections in 
one or more layers is not new. Research on the layout problem was initially done for printed circuit 
boards in the 1960's. ‘The layout problem for printed circuit boards is closely related to that for integrated 
circuits -- components are placed on a board and interconnections are made by conducting strips (wires) 
printed on two or more layers of the board. Wires on different layers are insulated from cach other, but 
wires on the same layer must not cross unless a connection is intended. Depending on the manufacturing 
technique, a conducting path may be able to change layers only at fixed positions on the board (called 
fixed vias) or anywhere on the board (called floating vias). Some researchers have used models for the 
layout problem which they intend to apply to both printed circuit boards and chips [Han76]. However, as 
we shall discuss below, the objectives and assumptions for printed circuit boards and integrated circuits 
are different. 

Traditionally, layout of circuits has been divided into two phases -- the placement phase and the 
routing phase. Separate algorithms have been designed for each. In the placement phase, components 
are assigned positions; in the routing phase, the paths which the wires will use are determined. For 
printed circuit boards, the goal is usually to minimize the total length of wire used. Generally, the board. 
is divided by a rectilinear grid and components can be placed only at certain locations on the grid. 
Terminals where wires must attach are at fixed positions on cach component; these positions match 
locations on the grid. Automatic layout systems for integrated circuits have borrowed the algorithms and 
models from printed circuit board research and expanded on them. Most systems use standard cells of 
uniform size which are arranged in rows and columns, leaving a grid of horizontal and vertical streets in 


which connections can be made. Examples of standard cell systems can be found in [Fe76], [Per77]. 


. li = 

For placement algorithms, components are gencrally modeled as points and the board or chip as 
an array of locations on which the components must-be-placed. For-printed circuits, components are 
usually standard size packages 20 that size is not a facter,, Far igtegrated circuits, a size parameter may be 
associated with each point and a capacity. with cach lacation. For example, in [Sc76] cach location 
represents a row of standard. cells and has a capacity. which-is the length of the row. In the majority of 
placement problem formulations, the objective is to find a placement which most facilitates the routing to 
follow.. A functian of the placement is aoasa as an objective function.to be optimized. The function is 
supposed to be an indication of the difficulty of the routing given, the placement. Most often the function 
is an estimate of the total wire length used for routiag. . Various estimates of the. wire length needed to 
route one net (set of terminals to be interconnected) are. used, @g. the half-perimeter of the snalien 
rectangle enclosing all terminals of the net; the: length of the:shortest.spanaing wee of the net... The 
estimate must be eaay to compute. 

There are two types of placement algorithms -- constructive initial placement and iterative 
improvement. Most automated layout syatems. use. both, although either can be used alone: iterative 
improvement can be used with a random initial. placement. Constructive initial placement algorithms 
place components one by one based on, their connectivity. to.companents previously. placed. In many 
systems, a designer may choose and place the first components. When: the first component is chosen by 
the algorithm, a special component such as a bonding pad may-be chosen, or an arbitrary component may | 
be chosen. Given an initial placement, iterative improvement algorithms move compencnis to improve: 
the objective function. Thesc algorithms are local optimization algorithms, One. common technique is to - 
exchange pairs or rearrange larger sets of camponents and test whether the valuc of the objective function 
has been improved. In another method, called Force-Directcd Relaxation, a component is moved to a 
point where the "forces" due to its connections to other components are balanced. A description of 


various placement algorithms of both types can be found in Chapter 5 of [Des72]. An expcrimental. 
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comparison of several of these algorithms is presented in fHati76} 

The routing portion of the layout problem has received much attention in the past ten or fifteen 
years. Many algorithms have been developed to find “near optimal” stlutioes to the: routing protilem 
when sets of points which must be interconnected are given as input. Most of the algorithms attempt to 
minimize total wire length. These algorithms almost exclusively use the’ rectilinear (also known as 
Manhattan) measure of distance. 

Routing can be approached a number of different ways. Nets may be connected using Steiner 
trees, i.e allowing wires to branch at points between terminals; ov sy Be feaniichead to spanning trees, 
where no branching wires are allowed. (Sec Figure 2.3) Covance pith ‘may be allowed to change 
layers at arbitrary points, at fixed points, or not at afl. The tree of contections for cach net may be 
determined before the actual paths are found -- called wire fist determination -- or the trees may be 
determined as paths are laid out. When each path must lie totally on one tayer, the laycr assignment, i.e. 
determining which paths will He on which layer, is offen done firsi. Ali paths on a layer must be routed so 
that they do not cross. However, many of the test recent routers for integrated circuits and printed 
circuits use the restriction that horizontal wire’segiticnts ‘and vertical wite segincnts are’ bri'sepiurate ryt 
(Horizontal and vertical are the perpendicular directions of the rectitinear metric.) Paths are composed of 
horizontal and vertical segments. Fach change of direction is Scares of layer. Typicatly, two layers are 
used -- one for each direction. In the projection of the layers on one plane, two paths may intersect within 
perpendicular segments without being electrically connected. “When paths can change laycrs only at fixed 
points, called vias, the vias may be assigned to perticutar conricctions before routing. For printed circuits, 


the number and location of vias is often restricted due to the fabrication technology. For integrated 


1. Under the rectilinear metric, points (x,,Y,) and (x,,¥,) in a cartesian coordinate system are distance 
bk -x51+ly, Yl apart. _~™ 


-3- 
circuits, vias are contact cuts [Me80}. They are usually allowed anywhere, although it is desirable to 
minimize the number used since they require extra area. 

Most algorithms to da the actual routing of wire paths fall into one of the following four 
categories: maze routers, line routers, cell routers, and che conan The first algorithms used on 
printed circuit boards were mace routers. Maze-routers find one path at at time. They dre based on Lee’s 
algorithm [L.ce61] for finding te shortest path between ma nodes in a graph. In fart, a path can be 
found between methbers of two sets of nodes rather than between two specific nodes. The graph used for 
routing is a grid graph containing. fgbidden regions. The slgoithm is a breadth first. scsi of the usable 
grid points. (See Figure 2.2.) Multiple layers may be modeled by using one two-dimensional grid for 
each layer. Each grid point (for unrestricted vias) or each of a special set of grid points (for fixed vias) on. 
one planar grid is adjacent to corresponding grid points on other planar grids. 7 

Cellular routers. [Hit69] algo divide. up the rousing. plane into: a dental nemagemdnt of cells. 
Again, each path is found, one at a time, by a breadth first search 0 over al cell Bach ceil represents a 
region of the routing surface, including all layerk: This represesaton Sr obtined by projecting all layers 
on one plane, then partitioning the plane into pieces of uniform size and shape, These pieces are the cells, 
Each cell is age enough to fit mere than one wire width on: each layer. The rovting alghrithm must 
define the entrance and exit points ofeach routing path inthe bet of cells traversed by the path From a 
given point on a.cell boundary, the rentable pointe Side the cd's boundaries are ealculated taking into 
account wires which have, been previously routed and have ‘cut off parts of the cell from other parts 
is alice The determination at the actual paths used — oll is left for-a ‘second algorithm. 
Cellular reniterside-not roquirc: ‘the sera arumnber of grid pein wih must be rpraenie when using a 
maze router. : 
- Line routers [Hi69] do got partition the plane. These routers work directly with horizontal and 


vertical line segments. They build up a path betwcen twd. points out of these segments. A-line router 
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Figure 2.1: Steiner tree councctions versus spanning tree connections. — 
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Figure 2.2: Lee's Algorithm. 
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finds paths only in one plane. Paths area determined one at a time. Given two target points to be 
connected, line segments are extended from the targets until they hit obstacles. Points along these 


segments from which a nage be extended _ wil cui an anna Beviously hit are found. In 


Bu 'e 


this way, “escape lines” are encrad unil a path been the two target pois can be constructed from 
segments of these escape lines. (See Figure 2.4) of e 3 4 


The channel routing technique las?) wag developed for routings with horizontal ‘id vertical 


wires on separate layers. The area. state! for 08 Ns divided into —— and vertical streets. 


Each horizontal or vertical street ¢ can. contain . umbe ae horizontal vertical path segments, 
respectively. Paths are first found through the recta without regard confi within each street. After 
ail connections have heey pobally. routed. sheoigh she. sirecis, the. regret within each, street are 


arranged so that no two overap. ‘The number of pare! lanes or chimel in each street is minimized, 


i ‘be weed from each iptminal out into the street. 
ents te 16 allow a path, Sefenent, to ‘change channels. 

(See Figure 2.5.) The global outing problem is actualy path finding, problem on a graph where more 
that one path may use the same la. The local routing problem (or chianel assignment problem) 
resembles a packing problem rather than a path-finding problem. These problems will be described in 
more detail in Chapter 4. eee a 
The algorithms described above are those most ofen used by layout automation systems. They 

do not encompass all algorithms. ‘A survey of routing algrithms canbe fend in Chapier 6 of [Des72] 
and in [Hi74). A study of router performance is cece i err. ‘None of the ‘algorithms are 
guaranteed to find an optimal eon pica they are alowed to ie. up and reroute -- exhaustively 


searching all posible ‘routings. Furthermore, ‘if the amount of area . which c can be used is bounded, the 


1. The use of the term channel has become confused in the literature. It is used by some authors to 
denote a strect, by others to denote lane within a street. We follow the terminology of [Has71]. 
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Figure 2.3: Cellular routing. 
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Figure 2.4: Line routing. 
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algorithms are not guaranteed to route all connections, even though such a routing exists. In practice, 
failed connections are completed manually. Many systems use a combination of algorithms, beginning 
with faster algorithms to route mogt connections and = slower but more successful algorithms to route 
the final connections (e.g. see {Po79). Some gor ‘ee aanele saceraible to the problem of failing to 
complete all connections than others. When cach point-to-point connection is assigned to a single layer 
and then routed, the paths found first may encircle a terminal ‘not ye connected ad anything. It is 
impossible to find a path te this terminal without changing Sayers oF or spina wp ciety routed path. 
(See Figure 2,6.) For routers of this: ype, the-order. tf-which the paths are ansigned f for a b probe can 
make a great difference in whether the algorithm i | aces all the connections. Given this, 
it is ieermins that Abel [Ab72] caubctadon from his entpiricat-stisdy that, overall, the peformanc of such 
routers is fiat significantly affected by-various ordering criteria: (The router used is a Maze router with an 
added heuristic so that paths do not run next to a row of terminals at minimum spacing from the 
terminals. The paths being avoided would lead to a large number of blockéd terminals.) The channel 
routing technique is less susceptible to the problem of failure to route-all: eannections ‘siaee the exact - 
jonition of each segment is determined ater all paths are globally assigned. In fact, if area is not limited, 
the streets can contain an arbitrarily large number of unos and ve routing can always be obtained. 

Several characteristics of the models used for past research on the layout problem are very 
restrictive, especially when considering the layout:of very large scale integrated (VLSI) circuits, where up | 
toa million transistors are packed on one chip in a very. diac. arrangement. The use of cells which have 
one or both dimensions fixed limits the type of scaupontet which ean be used. One may imagine having 
in the same circuit a very large component which is an array of registers and a number of small 
components realizing a spccial function. Once the components being used are of varying size, it is a waste 
of space to place them in an array scanned by “streets” fina the length andl width of the chip. The 


size of each component and the way the components fit together on the chip are important. 
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When an array of locations is used for cell placement, it is teasonable to use total wire length as 
a measure of the worth of a particular layout. Different arrangements of components within the array 
only affect the area of the layout insofar as they affect the interconnecyon paths taten. The total wire | 
length is an easily computed approximate mesinire ofthe pmount of routing sofae which has been used 
at any point while routing is being done. This, in turn, é an acne mans of the congestion of 
wires in routing areas. Congestion may affect t the rota of interconnections still to be made. When 
the area used for routing between components is fred, which is usually true for printed: circuit boards 
too much congestion may result in falture to route some of the interconnections. When the area used for 
routing can be expanded by moving the components (While keeping their relative positions fixed), as is 
more likely for integrated circuit design, too #tuch Conjptstion may result in a larger overall layout size. : 
When components are not of relatively eroiricg! ns mer placement has great effect on total 
layout size. Wire length is no longer a good prnaeeeae to layout size. This is illustrated in 
Figure 2.7. The placement of components and routing of wires interact in a much more complex manner 
to determine the total size of a layout. The most recently developed layout systems no longer treat 
component size as a parameter which is of secondary importance. Components are modeled as rectangles 
-- giving them shape as well as size -- rather than as points with a size parameter attached. When 
_ rectangles must be placed on a plane, the way they fit together influences the area used by the placement 
and the shape of the spaces left for routing. Preas and Gwyn [Pr78] have retained a constructive initial 
placement based on connectivity but place rectangles in a plane rather than points on a grid. Their 
iterative improvement phase tries to minimize the area of the circuit by selecting as the candidate for 
placement modification a component one of whose dimensions contributes to the widest part of the 
layout in that dimension. This component may be rotated, reflected, shifted, or exchanged with another 
component. A modification is accepted if the area is reduced. In [La79}, the initial placement is 


produced by dividing a square of area the same as the total area of the components into rectangles of the 


Figure 2.7: Tradeoff of wire length versus area. 
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same area as the individual components. This gives an approximate placement which ts modified to fit 
the actual components ipescremests are ee to the placcment by using rotations and shifts. Routing 
is vabes into account in the improvement phase. ‘The width of each street is estimated based on global 
routing and included in the area calculation Brinkmann (876) lo uses the technique of. dividing a large 
rectangle into smaller ones to. find an  spprimate placement. 7 

Rciiting programs for the most recent system aloo ry © ve veal area rather than wire length. 
Chansel routers are used in pr and par Ca router are eaay to use when arca is the parameter 
to be ¢ optimized since local routing anomie street, with. Svea can be allowed to shrink or expand a8 
neeee to complete the routing In [L079 component are stowsd bo eid positions above the routing 
are independent of one anther 0 tat aca nok wasted unnecessarily by aligning the edges of the 
compcaens. | ian 

Wire ‘brett remains, an maperiaet Varsiucioa for a yout because it diet altos the quality 
and speed of signal in me eee Ina simation such a tat shown in a Elec 2 where both size and 
total wire length cannot be minimizad at the same time, a tradeoff must be made. In this conten, it may 
i aaa ta dedi eagle mot anportant requiem ofthe eahing circuit sie At other times, only 
a maximum | wire length mien be  impored « on a certain Ineconnecton in the circuit. To make oe 
more compl, we — Ss a a) in which the Les cpeneeay was that two interconnections have 
approximately the. same length say when two ‘ontpem of one ‘component are the inputs to another 
component In shor the dered measur for the worth of» layout canbe made very complicated if 
enough factors are considered. Paya quatieies may Seri. 6 on a other layout properties such as the 
density of wires in an area. In [Agu77] and [Ru77}, a he igereg ss | deserived in which total power is 


minimized and timing constraints are observed. In INo7, panei la dela must be minimized 


can be designated “critical” and treated bea 


2.2 Alternate Approaches to Layout Automation . 


The approaches described in Section 2.1 separate the placement and routing phases of layout 
The input is a set of components -- either specified as points or rectangles - and a set of interconnections 
to be made either among the components as pciats o or ‘among terminal point on the commacinente ‘There 
’ are alternate approaches in the literature where the ponological spect of ‘apes are ended using uti 
theory. When finding the layout of a a circuit, placement and routing are considered together by lees a 
planar embedding of a graph modeling the circuit. ‘The models rie not eaciily associate one ae 
with cach cnponent: Each component may be modeled sa a set ot of ‘nodes and edges: es a cycle. 
Interconnections among a set of terminals may be modeled asa set tof edges or a set of nodes and edges. 
For ample the branch point of a wire ¢ may be ccaiescated: asa spain under a Steiner tree representation 
of connections. Several models have been suggested. A ‘summary can be found in van Graph 
embedding techniques have the oe that placement and routing interact chaiplesely at the 
topological level. However, geometry - ~- the size and shane of the components ~ ~ is is completely gnared 
and waist Be accouated for tepacuaty: ~~ . 

In the standard approach andi in the sh embedding approsch, the onl information Sante the 
layout provided by the input is the set of nets. ln an alternate approach based on n stick diagrams (Me80) 
topological information about the desired layout is ako etovided: In stick dina. regions in various 
layers of the integrated circuit are re represented. 2s nes. For example, ‘in nMOS ec technology, a . polysificon 
line and a diffusion fine crossing represent a transistor. The lave postions of component and the 
general path of each wire are indicated. The Inyout automation omen pan expand the “aicks” into 
rectangles of the proper dimensions based on i res, and modi ‘he sciik Ft) that components and 
wires fit together with the proper spacing. Examples of such sem are STICKS [wi77] snd CABBAGE 
[Hs79]. These programs attempt to pack the layout as much as possible while satisfying the design Tele 


and maintaining the original relative positions of components and wires. This technique exploits the 
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human designer's ability to do overall layouts, ic. rough sketches. The program FLOSS [Ch77] also 

backs a rough manual hs it works fram a hand drawn mien anneal ayeccns for iia the designer 

| does the gaia layout use a symbolic reprseuation ote the layout (Gite, Pert The symbolic Hayout 
is produced by the bei seid expanded into a fully speciied layout auiomatically Other techniques 

; for detigi automation use special structures such as Programmable Lowi Arrays PLAS) with known 


layouts. Functions) spéciticadins can be automaticaly converted into. PLA implementations lay 


Special interconnection patterns can.also be exploited to ‘sutiat in design (3079) 
2.3 Summary 


. All of the above ehniques have been developed to to automate, at ast in part, the seen of 
circuit Eu The research reported in this + dimertion ar restricted to placement and routing as 
described in me ease section. Although ‘many algorithms sais been designed sad tested, little 
mathematical sda of the algorithms has heed performed. ‘The techniques of pompenily theory a are 
not regularly applied to layout problems One ceases een is the work of So, Ting, Kuh and others 
((Ku79}, {S074}, [Ti76}, [T 78h rri79}, rTs79) for pened circuit board i routing In this work, a roms 
Sei for printed circuit boards with ai vias in columns is broten up into severe! problems, ending 
ae a auniber of i instances of a routing saibhen for a row of terminal ona sats layer. om problems 
are analyzed. Necewary and sufficient conditions for a single, singe ayer routing to be optimal are 
developed. ‘ieee the model is not ewe suited to ‘itearaisl circuits ‘The erat reported in this 
dissertation also focuses: on anita subproblems of the layout problem. The problems are motivated 


by the channel routing noe of interconnections. 
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As is encen 3 in n the discussion i in the previous I hap, most 3 models of circuits for layout use 
points or rectangles to shod components. I the research eed in n Chapers4 trough 6 we use 
rectangles. The model i is described procily below. ‘In the second section nof his chap, we describe a 


graph fale in which h components are points We discus its u use in proving g bounds. on he arca a equired 


by circuits with certain interconnection paticms. 
3.1 The Geometric Model — 


We have chesen to use a recta a ‘of ani pee it nai the ponte and 
topological vases of the layout problem, Components are :ectangul i in shape and variable is size. 


Wires lie on any var several — and. are of uniform width ‘The model was aided a te design rules 


WEG & 


ciessaned in [Me80] for nMOS ectnooy. sur it is sinended to be ppb to carats technologie 


af 
es the model i is as follows. 


% 


The input for a boii problem will const tof a a set of components a and a set tofn nets. Fach 


component C,, will be a a rectangle with sivend dimensions , andy, IY; is siven locations o on the c boundary of 


a a component Ca are terminals t., j= ou oN, where n, is ie number. of terminals on component C, Fach 


seh? oF UEg 


et is a set of terminals and the nets are pairwise : disjoint Bach | net represents a callecton of terminals 
which must be eet connected: ‘The you poem es i campenens on 12 pla vlan and 
form the interconnections specified w the nets in the minimum ea | The interconnections can 

be nade in any of N ies then Ni is specified : a Prior L Thelen ar numbered 7 aug N; ayer i is 


PM 


considered to be adjacent to myers i+1 and i-I, for 21<i<N- “I. Typically, Ni is tier “The wires weed to 


re , 


interconnect terminals have a uniform width w. ‘There is a minimum spacing between wircs on a layer, 
between components, and between a wire and a component of s. 


The interconnections form paths which lie on the surface of the plane not covered by the 
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components. These paths go between nodes. A. node is a terminal.or an addigignal paint on one or more 
layers of the plane. Each path segment hebinnts Gai nodes has a designated layer. The paths induce a 
partition of the nodes into disjoint sets of connected nodes; each set contains terminals from exactly one 
net. Path eguicais in different layers-may. intersect. Path eons in the same layer cannot intersect 
unless they-are: ‘gesociated “aith ‘the same net and intersect at a node. Additional points are allowed as 
nodes so that orte wire may spi into several and: $0 that a path may change layers. A path may. -change 
from tive ito layer j, for any i andj between | anil N at any node; however, te node mu Hie on all 
layers between layers and inclusive. hye path is coer oe on eth of hee er at the node. 


dive teas tang ote layers between i 


Therefore, no path aseociated with a different net eat ee 


and j. The set of nied and pane in any one layer represen an erent a a plana Boh The area of 
ce eas ae *eSige 


the layout is rae area of the salle mtg dn te wand 


ne: ‘all components and wires. 


Figure 3.1 ihstraten 


_ The model desciibed above takes into scoot bb tageometi and toposes aapct of the 


Gaede circuit layout problem. The wires are oghety sve via The restriction that components 


be rectangular is used sll o somewhat igi te rob sially the calculation of free area for 


routing. We believe that restricting components tive rang shapes sil provides a model that is 


applicable to the layout of VLSI circuits. For thostlgia components whose layouts would fill only a 


small portion of a rectangle; e.g. “L-shaped layout maybe tei brat the conponkl into two 


be! 


or three physical components with more nearly recigula you, 


One of the advantages of the model is that. 3 o various tevels of modularization. 


Components may represent anes logic gates, or even arithmetic unit allowing shy concen devel ; 


of design. A bierachical approach to Inyout ood as that used by Preas and Gwyn r7a] can conlly be 


taken using this model. 


Additional parameters can be allowed as inputs to the layout problem in this model. Upper 
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Figure 3.1: Must ration of the geometric modvt for layout design. 


spacing 
less than gs 
. different: . 
layers 


W,X,Y,Z are nets ST aE ag ots LAE) ondary oe ee 
O indicates a node. Nodes are labeled | by the | net t being « connected. 
4S indicates layer 1; indicates dayee. 2 

There are only 2 layers. 
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bounds on the length and width or the total area used by the layout may be given. Alternatively, a bound 
on the aspect ratio, Le. the.ratio of the length to the witith of the layouk area may: be:given. This would 
restrict the shape of the laygut. Aay of these: dmitations may he: necessary to insure that the layout 
produced is, usable. For example, suppose.a.minimam aree-selution resulted in an. extremely long and 
thin chip. The shape could prohibit the fabrigation or the:pagkaging of the chip. Another possibility is 
that the circuit being laid out is.a component at a higher: level in.e design hierarchy. ‘The bounds might be 
derived from requirements at the higher level. Bounds sush oa: those. pravide @.tange of solutions within 
which the algorithm should work. Jt is alsa possible 4a dnkendece-svire Jength sequiremeats to the layout 
problem. As mentioned in Chapter 2, wise length i4.an anpartans feotor in signal- quality. Individual nets 
may be. given upper bounds ea. the. sotal .wire- length Of, theig: interconnections. (Wire length can be 
measured using the center ling through 2 wise.). Relationsaips etwoenthe wire lengthe for various nets 
may also be imposed. | 

In the following chapters, we will put restrictiens.on- the types:of layouts allawed: within: the 
above model. We will restrict ourselves to rectilinear paths.and.orthogenal orientations for companents. 
Directions “horizontal” and “vertical” perpendiculat. 42 each: ather will be chosen, ad all: components 
will be placed so that their boundaries are. parallel te. one- of, these: directions; Also, all: paths will be 
composed of horizontal and vertical segments. This restriction ia imposed: fos two reasons... First, the 
problems examined are motivated. by currently used layout algorithms.-.As discussed in the previous 
chapter, the rectilinear metric is almost always used. Second, the sestrictian limits the number of possible 
solutions and makes the problems easier 4o analyze. ‘Most.af the time, we also restrict horizontal and 
vertical segments.to be.on diffcrent leyers. Given that: only boriaentel and, vortical- wire: segrtients are. 
being used, this assumption aot omly greatly simplifies the description of allawabie paths {those that may. 
intersect at a point but.do not overlap) but is reasonable when: only two layers arc-available. In this:caag, 


two segments running in parallel ov different. layers, one-on. top. of the other,: prohibit any path from. 
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crossing from one side of these segmenes to the other. “Atty pair of terminals, one on cach side of the 
segments, which must be connected will require a path-Which ‘goes around the parallel segments. 
Therefore, one would rarcly want sich segments: ‘Also, if the scgments-are extremely Jong, they may 
cause clectrical problems. due tb capacitance. Restricting adjacent: layers to contain segments in 
perpendicular dircctions eliminates these potential problems ss 

There are scveral technical aspects of the Jayout problems which the model does not take into 
account. We discuss these below-and indicate How the model can: extended or modified to include 
them. In actual component design, a component may ‘have several terminals to which a particular 
connection can be made. . These terminals. may be cither physically equivalent, ie. they are connected: 
inside the component, or logically equivalent, i.e. you 'ehri't tell then? apart fimctionally. An example ‘of 
logically cquivalent terminals ts the set-of input terminals of an "and" gate. The problem of deciding 
which terminal to use in connecting a particular net is called the pin assignment problem for printed. 
circuit boards, and is usually solved before routing digorithi ‘isused FAR} Both physically equivalent 
and logically equivaient terminals can be modeled as sets Of terminal. A set of equivalent terminals, 
_ rather than an individual terminal, would be a member of the set defining a:net.-'A sct of phiysically 
cquivalent terminals: would appear in’only one‘net, whites St off logically cquivatent terminals would 
appear in several ncts, but no mere nets than the number of térntinaly in the’ set”. When layout ts 
complete, each set of interconnected: tehsil wvust-Gotiods at -sivest One’ terminal from each stt of 
logically equivalent terminals. ' This model of equivatent terminals ts used by’ vati Cleemput [van76]. 

‘in integrated circwits, layere are made ‘of different ‘matériais snd ‘different design rules may - 
apply. Allowing the width of wines ‘und: separation’ betweeh objects t Vary: between layers is a minor 
modification to our model, although the resulting layout problem is more difficult, However, the design 
rules for various layers can be more complicated ‘than’ the ‘model will allow.’ For example, in nMOS 


technology {Me80}., wircs in diffusion and polysilicon layers.cannot cross. Therefore, we model them as 
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one layer. However, the ntinimum distance between two wires a either polysilicon or diffusion is 
different from the minimum distance tetweon wires of different layers. This Is not provided for in Me 
model. We must wae the tirgest of tite actual required separations, Luctily, the metal layer in iMOS can 
cross both diffusion and polysitionn wires. If it coukt only tives one, then the actual layout problem 
would net be captured by our model: two layers could be tied for intetconncction, but a concept of 
“coloring” wires in cach layer accofding to the design tules would Rave to be added. Another design rule 
which we have not accounted for concerns changing ener When @-cordutting path changes layers, a 
contact cut must be made to: eomneet the tayerd tivciigh aw Rsulating Inyer: The contact cut requires a 
square area larger that the-width aie wits, ‘This is not taken into account in our area calculation. Even 
worse, some wires on @ particular layer may foe significantly wider then other. wires on the same ‘layer; 
depending on the electrical load on the wire. For example, the wires supplying power to'a large circuit 
usually look: like.the human arterial: system: there: is a Serge itaifi’ wire and a tietwork of wires of 
decreasing size-scaching every part of the circuit. We reoogztize tfiat our model overtooks many detaits of 
actual layout design, but we fect that it s «reasonable approximation Of the major tseues fn the layout of 
3.2 A Graph Embeddieg Model 

We now discuss the uses of a model in which a circult is represented asa graph. We call this 
model the graph entbediding model. Each component is represented as a node-in the circuit graph, and 
each consection to be made is represented as an edge between two nodes. “Thelayout probiem is defined 
as the problem of embedding the circuit graph in-a two-dimensional grid: graph. An embedding maps 
each node represonting a component w a node of the: grid. ‘This mapping is one-to-one. Fach edge 
representing a connection is mapped to a path in the grid graph. ‘Fite path can only contain two grid 


nodes which correspond to compenemt nodes ~ thosc that correspond w the ‘ehdpoints of the edge being 
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embedded. Paths corresponding to distinct edges are allowed to intersect at grid nodes ‘but are not 
allowed to use the same grid edge. ‘To make such an embodding possible, cach node in the circuit graph 
is restricted to having at most fowr edges adjacent to it, Le, it is restricted:40 be. of degree at most four. 
Such an embedding is an edge-disjoint homcomorphic embedding. Given an. embedding, there are two 
measures of area which we will use. The first, which we-will call the.nade area and denote A,,, is:a count 
of the number of grid nodes used as sate cme taats baal oh paths which are images of edges. The 
second, called reciangle.areg and denoted A,, is the nunaber. of nalice contained in a rectangle. whose 
boundaries lic on grid cdges and which circumscribes: alt nodes used. inthe embedding. i.e. all nodes 
counted by A... For a circuit graph, C, let A,(C) and A,(C) denote:the minimum node area. and 
rectangle arca, respectively, over all embeddings of C. Obviously, A,(C)<.A,(C). An embedding is 
shown in Figure 3.2. | 

‘Fhe model presented above can be viewed a6 describing a layout which uses only ‘i catsiaaaa 
vertical wire segments (segments in the two grid directions), each diroction on a separate layer. 
Alternatively, the nodes of the grid may be viewed aa representing unit asca squares ana plane. Each 
_ grid edge adjacent to a node represents a boundary of the corresponding square. Two panes) 
node need not actually intersect on the plane -- they may run diagonally, cutting across $ opposite comers 
of the corresponding square. (See Figure 3.3.) Thompson [Th80] uses a des ceria layout ich, 
divides the plane into such unit squares. He views connacctinns as running primarily in a-single layer of 
metal. Crossovers.are achieved by using short.runs in a second layer such as polysilicon. | 7 

The graph embedding model lumps all. terminals. of a component. into ene point... The 
point-to-point connections made by wires are pre-determined and represcnicd in the graph for the circuit. 
If Steiner tree interconnections of the terminals are desired {i.c. branching wires), these must be explicitly 
modeled in the circuit graph by adding 2 node for osch branch point and eriges from cach brancix poiot to 


components or other branch points. The model does not represent the arca required. by components or: 


Whe RRR Pr ER 2 Sehn s uphig Seis i Te PT ec aligtegeae ghee ds op aeameagy ae ae te ei 


~3l- 


Figure 3.2: Embedding a graph. 


CIRCUIT GRAPH Cy 


Ay(C) = 17 
ACC) = 20 


> | Biggie Sale Refeets Of the order of ternsinals arewad s component. 


Stes 


Figure 3.3: A node representing a square. 
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the fact that terminals have a specific order scocal the component. This order causes some routings for | 
the graph embedding model to be impossible in the geometric model (see Figure 3.4.). However, the 
model is very useful for investigating the  inplications of certain interconnection patterns, ~ 

When node area is uses, cach circuit none een mene eae: one ann ‘a ‘area. Therefore, any 
_ circuit graph with | n nodes which requires more than O{n) area: must “have an iaterconnection pattern 
which regu a tot of are wo route, ‘Lower bounds on the amount of arta required by a graph are 
proven using this meer as Note that any lower bound on de area salves bound on rectangle 
area. Thnciens [Thsoy proves a a lower bound on the arca required to exnbed a | graph asa function of the 
minimum biseetion width of a subset of the nodes of the graph. Given a graph G sia eels S, of the 
nodes of G, a set of edges in G.bisects S if the removal of these edges partitions the nodes of Giinw two 
DQ LISV2J of the nodes of § are in one. ect and £4§}/21 are in the other set; 

ii) After feitusting the edges, there are no paths between nodes in different sets. nae 


Let we be the size of the smallest set of edges bisecting S in G. Then: 


Theorem (Thompson fensop: Given a graph G with nodes of degree at most fout, 4nd a subset, S, of the 


BEecet Wag OES 


An msuperconicentrator is a graph with n designated ingut nodes and nding output nodes 
such that for any nets of k input nodes, and k output nodes, Igk<a there are k ns disjoint paths 
connecting me k input nodes to the k output nodes [Val75]. Forte .set-of input nodes, ik w, 2 bn/2J. 
Therefore, for any n-superconcentrator, G, A,(G)2> (0-116: “since for any n, there are 


n-superconcentrators with O(n) nodes, all of bounded degree [Pi77], there are graphs with n nodes which 
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require node area Q(n’). Note that altheugh the superoonceatrators in [Pi77] do not have degree at most 
four, it is straightforward to reduce the degree of each node by adding a node for each edge. The number 
' of nodes added is then bounded by the original number of edges. Figure 3. 5 iflustrates. 

Thompson also defines average and worst case infonnatéon complexities of a function. He 
derives a lower bound: ‘on the oo or worst case time ee router by: a raph t compute a function: 

iperaantors worst case) tine > 
(1/w, Maverage (respectively worst case) pimumaadibe complexity of the function) . 

where | is a special set of input nodes in the graph. Combining the results, Thompson obtains a lower 
bound on A, x(average time)? for a graph which computes an n-point discrete fourier transform of 
Ln/8.J"log’n; he derives a lower bound of Q(n7log’n) on A,x(worst case time)? for a graph which sort n 
numbers. The reader should refer to [Th80] for details. Hounds for other functions have also been 
derived by various authors using Thompson's technique, e.g. [abst0j (Sav 

Upper bounds can alo be obtained on the area to embed Various classes of graphs. Upper 
~ bounds are derived for rectangle area, A,. Any such upper icuad is apo an upper bound for node area, 


A First observe'that any graph with a nodes, each of degree at mst four, can be embedded in 


N’ 
rectangle area at most 6n*+ 3n. We sive here a modification of thé proof presented in (Val79]. This 
modification improves the bound from a (3n)x(3n) square area toa Gn)x(Qn+ 1) rectangular area. Since 
n-superconcentrators require O(n?) area, proportional te a% area is both necessary and: sufficient for 
embedding an n-superconcentrator. 

The embedding ‘achieving rectangle area of at most 6n?+3n is shown in Figure 3.6. The 
directions used below refer to the figure. The nodes are embedded in onc vertical column of the grid -- 
one node every three grid points. There is a column for each edge, cither to the left or the right of the 


_ column in which the nodes are embedded. The path representing an edge must reach the column for the 


edge from cach node representing an endpoint of the edge. Therefore, each path includes two horizontal 


Me 
_ Figure 3.5: Reducing tbe degree of nodes in 2 superconcentratar,( 
In the 
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segments between the column containing the embedded nedes aad the column for. the edge... The first. - 


(last) segment of the path is cither one of these segments or a vertical edge from the endpoint of the path 


to the grid node Just above or just below it Ifa vertical edge begins the path, the horizontal segment 
going left or right: from the second nade is eed As long as th edges of the circuit graph can be 
partitioned into two sets = those wihied columns are-to the: left ofthe nae column, and those whose 
columns are to the fight ah that at mos een iaeaees the same set, such an 
embedding exists. The odes of any graph whage node are of degree at moet Four can be colored using at 
miost six colors such that fio two edges adjacent ia ibe same node. are the kame color [ne66, p-262}.° Let 
edges with odd numbered colors have columns to > the left, ted those with even numbdred colors have 
columns to the right. At most three edges adjacent to a node onithe sathe side, as desired. Since any 
graph with n nodes, all of degree at most four, can have ‘ most. 2n edges, there are at most 2n+1 
a, ba = a poy eny ‘ : 

Other classes of graphs can be embedded ines shag 140 area. Valiant [Vat7s] and Leiserson 
[Lei80] have independently shown that sven a graph G with : nodes, each of degree at most four, if G is 
planar, then A,(G) és O(nlog’n); if G i Be tree, then Ag(O) i : On). Valiant actually shows that an O(n) 


a 


embedding for a tree can be achieved without : 10% ng a crags. It s an open question whether 
there is an n node planar graph which requires O{plog?n) arez. Leiserson proves a general result which 
relates a separator theorem for any class of Gaskets an uppers: bound. ee. the. rectangle area required to 
embed any graph of the class. | eae 

The results which we have reviewed above iflustrate that the graph embedding model is useful 
for proving bounds on layouts for particular classes of graphs and for proving time/space bounds for 
function implementation. Using it, we can identify easy and hard interconnection patterns to route. In 


the rest of this dissertation we will be interested in algorithms to actually do the layout. Therefore, we 


will use the rectangle model described in Section 3.1. 
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Chapter 4: Complexity of Layout Problems 


Regardless of the exact paiapmiss of the ayout probe we are uierened in nding an 
efficient mer which c computes an optimal payout if this i is not ot posible, we would like | an efficient 
algorithm which computes an optimal layout much of thet time ¢ and a 2 good layout me rest of the t time. 
This. algorithm sia actually sa a siieciae of goin to wove eborobins which tiatter give « a 
Layout, Again, we would like the algorithms for the  subproblems to find solutions efficent However, 
most sue aaociated with trou layout are NP complete The definition of an nes 
problem i is given below. From a practical point of view, the NP- completeness ofa problem indicates that 


it is gaa a soansied to find an efficient algorithm w which solve the problem 
4.1 NP-Completences and‘ Heuristic Algorithess 


Two major classes of probiems in bina i thoory a are ethe classes P and NP. A paren is in P 
(NP) if there is a deterministic (nondetenminisic Turing machine and a oyna p _ cg the 
Turning machine ae any instance of the problem with: an a of length nina auniber of steps ‘i 
greater than a} The eae of an nee is the length of its representation 2 as a character sales ina 
predetermined character set. We will not formally define i Turing machines here. The interested reader 
should see [Ah74]. Very often, a nondeterministic Turing machine solves a problem by “guessing” a 
uo and testing to see if its guess is actually a sohution. ideas the problem of intrest is a 
minimization (or saiualtana) problem, it is reformulated » 30 that candidate ssnutions can be tested 
independent of each other. In a new problem, a parame k is part of the input A A solution to the new 
problem i isa feasible sini to the old problem for which: the ant to be minimized (maximized) is 
less than (respectively greater than) k. A feasible solution is one which sstiafies all requirements of the | 


old problem except optimality. For ee if the orginal problem ‘ to find a minimum area layout of 


a circuit, the new Sibel. given the circuit and parameter k, is to find a Layee of the circuit of area ‘loss 


than k. Under this formulation, a featible sobision: cat, ‘be’ bested -by criteria depending only on the 

feasible solution to see if it is a sofution. When the prompt or metron of the ia is ase 
feasible sonitiobe must be compared agains cach otter to find the scsi solutions. 

‘The scipabee of steps taken by a deterministic Turing 5 machine is polynomial re related to the 
number of 0 seps taken under a model of computation that 9 coreponds to the instruction set bof a a 
combiner if the lengths of numbers operated on is taken i into 9 account tana) “Therefor any ‘problem 
which i is in P can be solved i in a ‘polyboetat number of scps sbyan ao lgoritm j ina igh evel programming 
language The number of aces executed yy an algorithm is referred to ‘as the time taken a the 
algorithm. We wind like all the problems we need to solve to abel in 1P. | 

The question of shethcs P = NP is one > of the « major open eer ee in cael) nai ‘itis is 
believed that P # NP. Probicms in NP have been found Whose thopticrateip: int P ‘implies. NP = P. Peete: 
problems are called Nasconnele propien: A problem is NP-hard if se cach saute in NE, there is a.. 
soiyiceita Panda transformation conpiitsife by a determina Toring machine in a a polynomial 
number of steps which transforms an satire of the problem i in 1 NP with an input t of length n to ose 
instance of the NP-hard sroblern with an input of Feng Am An NP-complete problem is one which is . 
NP-hard and i is in NP. There are no iwi deterministic sorts fo NP-complete problems wich 
take a number of = polynomial i in the i length of the input. The fact tha such wel sudied problems a as 
integer programming ani the  avelng salesman problem a are © NP-complete ious aongly suggeats 


tig: iy 


that NP # P. Therefore, proving a problem NPcomplet is ety sro © evidence t that any lord 
which solves the problem wil sais time consuming. ‘an . | | | 

The most common way to prove that a problem is NP-complete is to find a eduction from = 
ae NP-complete sible to the new probe whic can be executed i in 1 deterministic polynomial | 
time. Then, since the conpoudon of polyno is + still a olomia all problems in NP. can be  reduoed 


anes the known NP-complete probiem t to the new ¥ problem in deeriniie poo time, 


-%- 


When a problem has been proven NP-complete, the usual course is to try to find an algorithm 
which runs in. polynomial time and finds.a solution meet‘of the time. When the problem preven 
NP-complete is the parameterized version of an: optintization problem, ant algorithm which can find 
solutions close to optimal, if not optimal is desired. Heuristics are developed which can direct the 
algorithm. towards a good solution. ‘When the ‘algorithm is-searching through a large (exponential) 
number. of possible sohutiona,. toasting remove large: numbers of possible sulutions based on the 
likelihood that none will be optimalactutions. , 

Very often a heuristic alert is tested and validated by empirical evidence. AH of the 
algorithms discussed in Chapter 2 are heuristic singin Phesment and ieciiicg All have been 
judged by comparing the solutions they produce to manually cones ‘cnutions for the same se plese 
These algorithms are also compared te-each other to jpadge their worth. Ineontrast, inthis thesis, hcuristic 
algorithms for optimization problems are judged by the relationship of the solutions they ‘produce to 
optimal solutions. Consider a layout problem in which ans area is eed ers area aig) be the 
’ area of the layout for circuit C found bya pat algorithm. ‘Let area x) be the minimum area 
layout of C. Then define the worst case at uals ed of te algorithm, denoted we mo} as the maximum 
over all circuits of size n of aed yy(CV area, (C). Define the average case perfomance of the calor 
avg,,,(N), as the average 0 over all circuits of size n of area (CV area, ©, shee the average is taken with 
respect to a predetermined distribution of circuits. The size of a tireult-can de defined in various ways 
depending on the circuit model, e 8 the number of terminals, or the: sum n of the number of components 
and the number of nets. The size should be defined 30 that the length of the input specifying the circuit 
to the algorithm is polynomial in the size. 

Average case scvetiams is more likely to correspond to the observed performance of an 
algorithm, especially if the average is taken over “realistic” circuits. However. it is often very difficult to 


analyze. In this dissertation, the analysis is limited to worst case performance. If a lower bound on 
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afta, and. an-upper bound on area, can be derived, .an:upper bouad on wo can. be €oncluded. 
| Ideally, a bound of 1+e for gmall ¢ is desired. cs we are happy with apy constant bound. 
Hewever, as will be seen in Section 4.4, even constant ‘bouads.are not. achieved: by some common 
algorithms. 

In the next section, we wil review NP-completoncss sevults fir problems related to layout.. in 
Section 4.3, we will ‘consider : two. subproblems. of hail: witha caddie that they. are both 
NP-complete. In Section 4.4 we will ae a heuristic algorithm §oc.one: of the problems shown 


NP-compicte in Section 4,3, 
4.2 Placement and Routing: NP-complete Formulations —_ 


In this section we consider the compicxity of the problems resulting from. the decompasition of 
the layout problem into placement and routing. 


42.1 The point model: Some known results, 
First consider the placement problem used or rind i circuit boards and standard cells << 
OA Met AE 


thal components are e modeled as point and d toa wire eng * te @ ua to 0'bE minimized ov the 


layout The quadratic assignment problem i is one formulation’ maf pei 


: n : es eee ABS. AY, 
Given: Locstons 1 through m and sie natin {gl Sh.) between leatons; 
components 1 ehrovgh n and connetion mat fe} $id Sn) betwen componens, . 
ae hoe benfiok od boosie of ome 


Find: A one-one mapping. ,of components io lcaion such that, me 
COST(p) = vgpnitalid lotacmsbavalel io is minimized. 


Sahni and Gonzales [Sah76 prove that the tage rae borin sxigament problem is 


NP-complete. In fact, they prove that unless NP =P, there i is no approximation algorithm for quadratic 


~ 
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assignment for which there is an ¢ > O-such that: for all distances of the problem: 

COBH OOSTIBA) ST +e. 
The proof does.rely on instances of the preblent Witty cofnection ‘matrix values c;, greater that eft. 
Consider only instances of the problem for which thec,, are titted to a bounded range of values. The 
proof that the existence of an approximation algorithen for which dié‘ratio'of COS Tip, |) to COST(p,.) 
is no more than. 1+ ¢ impties NP = P no tolnger halls ftir ait'¥ However, the restricted problem is 
NP-complete as long as the ¢, are allowed to take on the value 0 or 1. 

The quadratic designment problem is a formutation ofthe placement problem in which all 
point-to-point camnections are apecified. ‘Fhe cost ¢y represent thie nuiiber of connections between 
components i and j. The distance dis the cstimate’ of the wite fefige needed to connect territinals at 
locations k and'h. Ifa placement problem formulation-tr whleh an esiliate Of the length of wire needed 
to connect whole nets is used, the quadratic assignment problem i is a eeeecalic case in none all nets i are of 
size two [Sah80}. Therefore, the ques mee wippiel rcs | to this formulation of the 
placement problem. It follows that this. formulation of the placement problem in NP-complete, 

Let us now turn to routing Ifa Steiner t tree interconnection Langan is desired for each h net, then 


dt PAS 


even finding the connection path for one net is NP-complete Formally: 


The Steiner Tree Problem 
Given: A set of points, P, in the plane with integer coordinates. 


Find: A set of integer points, Y, sach that the minimum fength spanning tree of PUY is minimal 
over all sets of integer points containing P. One. of twa measures of distance may be used: 
(i) The discretized Euclidean length: r ((x ray toy yyy where (x ry) and (xpy9) a are ihe 
points. The problem i is then the Discretized Euclidean Steiner ‘Tree Problem: 

». GD The rectilinear metric: fx pelt yyysh giving the Reetilincar Steiner Tree Problem. 


For either metric, the Steiner tree problem is NP-complete. If the standard Euclidean metric is 


Tiassa baie aati ii de pe SE OEE CDE ei RR er ent St a ot tr nc ae RR 


used, the problem is NP-hard, but is act known to bain NP (GarP 

The minimum spanning tree using cither metric. canbe solved in polynomial time [Ah74]. For’ 
the rectilinear metric, the. length of the minimum: spanning. tree is at most 3/2 the Jengsh of a-minimal 
Steiner tree flw78}. Therefore, any algorithm to find. minimums, length spanning tres is a heuriatic 
algorithm for finding sninimumJongth Steiner, wees with a. wors case rato af tongth over length, of 
3/2. A.discussion of heuristic. algorithms for the Rectilincas Steiner ‘Tete. Problem: can be-found: in 

The above two problems apply t0.the mode! of layout in which componests are points and 
minimum length wiring is desire,. Other problams reatad toithis.movigl.are-also NP-complete. Ting et.al. 
[1179] show that a via assignment problem encountered in their approach So touting is NP-complete. A. 
summary of a number of NP-complete problems associa with thin ode con be found.in{Sah8Qj. 


422 Te rectangle mae: seee rau. 


The sabi NP-complete sescies do not directly spy to the ‘model of yout in which 


4 hyd at 


cones are eae sad minimum area is $ desired. We sal now rove that even when no 


interconnections are needed, the placement of rectangular component ton minimize area is is NP-complete. 
Since this is a special case of the layout ken whes inercsanectloae are ss recuied, the more eee 
layout problem is NP-complete. The proof we present beiow does requixe some. exteasasumiptions in 
addition to those given in the description: of-our model im! Chapter 3.’ ° Att’ componetits and the 
circumscribing rectangle emerson ees ee re of their: sides is in the 
direction: of one of two: ac aepioaal axes. “This des pt ime deHictions on the placements 


but i is consistent with the concept of “horizontal” and "vertical" ‘being special directions which wires 
follow. We algo restrict all peints and dimensions ta be. integer valued aq: that the problem: hae discrete 


solutions. 


AbD - = £ 45 es eave = SR este? PAE ee ee ee es BAN gS MER Sigs le eae 
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Problem P1: Discrete Inyeut with no interconnections 
Given: A set of n rectangles and an integer, A. For }¢i¢n, each rectangle: r, has dimensions h, 
and w, which are positive integers. 
Question: Is there a placement of the rectangles og the plane with a cartesian coordinate system 
imposed sothat: | | 
(i) Fach boundary is perallel to one of the coordinate system axes: 
(ii) Corners of the rectangles lie on dnteger poihts in the plane: 
(ili) No two rectangles overlap; he 
(iv) The boundaries of any two rectangles are scparated by at least a unit distance; 
(v) There is'a: rectangle i in the plane ‘which circumactibes the placed rectangles, has 
boundaries parallel to the axes, and is: of avea. at: most A. The boundary of the 
circumscribing rectangle is allowed to gontain Boundaries of placed rectangles. 


Leman 4.1: Problem PI -- discrete layout with no.intercannéctions ~ is NP-complete. 


Proof: Consider only placements for which the lowest lemon corner of any rectangle i is at (0, sed All 
other placements are just tracabaciie of these. The coordinates sid a lower left corner of each fertangle 
and the orientation of the rectangle. i.e. whether the side of enh h, is in the x citecten or ie : 
direction, determine its position. Since the coordinate of = meee left corner can ony take on integer 
values between 0 and A-h,-1, there are at most A” choices for each ae A nondeterministic Turing 
machine can guess a possible placement and write it down, Given. a ee conditions (i) through (v) 
can be checked deterministically in polynomial time. This shows that iy is in NP. 

| The prot that Pli is NP-hard is accomplished o sites the Bin Packing Problem to it. ‘Since 


the Bin Packing Probicm is NP-complete [Gar79], this proves that Pl is NP-hard. 


The Bin Packing Problem: 
Given: A set of n items, each of size c; a positive integer;. also, positive integers B and.C, the 
number of bins and bin capacity, eee: 


SBe iby Ss Sie aT 2 ee 
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Question: es a ee 
. Ml heme meslansd in bik Sosnen peced C2 

Given any instance of the Bin Packing Problem, we opiate inane of Pl as follows 
There wilt be n+1 rectangics. One, called R, ‘wilt be of size by w, ‘where w = = 2B+1)C-1 and 
h = 2Bw+l1. The remaining rectangles will directly correspond to items, Rectangle #, will have 
dimensions bh, = (2B + I)<;}- and w, = 1. -‘Phe-hound on: etes:seilt be: A: = 20hk+2Bw. Note that the 
lent ofthe np ote Bin Packing Probie is ta + hg gH TA ineasons of he etn 
an A canbe caused taking no more than polynomial nen this eng 

Sint ye sin aber lip tiamladeibareng tsa! only if there 
is an assignment of the items to 1S bins Without cicdeitiig Capacity C in acy bis. “Given a bin packing, 
Figure 4.] illustrates a satisfactory placement. ‘We sow anmue week ony plaremneal etinkrang €) chpogie fe} 


eocrcsponss a: Ieee bie pases: Figure 42 wil hata Without loss of generality, let the side of ne 


rectangle R of dimension h be in the y direction Let “NeW” denote towards ower numbers in the x 


direction and right” denote towards higher numbers, let "above" denote towards higher munber in the 


4 


_y direction and “below” denote towards lower number ‘A any pin, he width ofthe layout i the x 


at * 


direction is strictly less than w-+1. I ot then area 2 Mw) = hw 2BW41> A. Therefore, aft 


rectangles, f, fy Must ie above or below R. None cnet ni beyond he gi ses of so that 
soba we fork, pejed i * 


for each rectangle, some line inthe y ection inerecs the rectangle and. Norge, Ty ry oened | 

so that is long sie is inte y direction. Onteie, te dimension of he cueing recange in the 

direction would be at last sian pt, wh te aed 1 xo for the separation betweea 
rectangle boundaries. Then: ae | 

arca > w(h+(2B+ 1x;) = Wh+2B+1D>A “selec: gakév)is set satisfied. 

We now know that’all rectangles, 1, are oritited so tht tht long sides are in the x direction 


re . eeret teh 
As 


Any placement can be modifiod sighily without increasing the area 20 that the f, give tui shale aid 
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Figure 4.1::A plucement.givens 2 bla packing, 


. 2b as es eeaeiea | 


Each row w cofitains rectangles corresponding to 
items in one bin. Adjacent rectangles are scparated 
by ote wait: Sinee the c; of the bit sum to at most 


‘= Co]: | : C, the length of the row ina most | 


(B+ De; ”) +(# of rectangle in thet row -1) 


in the Pow 
{ h= 2Bw= 1 «ah ae. ein S.QB+1C -1 


~~ area = w(h+2B) = A 


w = (2B+ IC -l 


Figure 4.2: A packing given a placement, 


| AE os see ae | 
(edi | eee ey | 


im. rowscoresponding 4 '(T] GJ! 
if there is anything here, 
J. the area is ta large; 


$0 all r within broken lines 


| | rectangle the long way is too tall 
| | resulting arca is too large 


there may be £; Sat the top and at the bottom — 


ROSES ER 1 URS RE SERRE SSS 
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below R. (Any f, and i less than two units above or Deca Dee Ro merase each anes 
least one unit to the left and right. These can be shifted ws form the fist rows abn and Dee R. The 
next rows are formed analogously above the upper boundary of the row i above R and below ae ices 
‘boundary-of the row below R. Figure 42 ihsrates.) Fach row of rectangles can be considered the 
packing of a bin. lf the rows s correspond to a egal bin packing, we are done: Suppose there a are K rows. 
Then the dinchsion of the raimacdbing reeange in te y Alvectian | must be h-+ 2K, since there must be 
a unit space separating rows from-cach other and R. if K > RB, then arca > wth + 2B) = A, contradicting 
(v). Therefore, at most B bins are used. It remains for a to eu that the sum of sizes c; of the items in 


any bin is at most C. Each item corresponds to a rectangle with h, = (2B+ 1,1. 


% ((2B-+1)c,-1) + the number of rectangles in the row - 1 


rectangles in row : ; 
S$ the width of a row of rectangles < w+1 = C(2B+-1) 
giving (2B+ 1\(sum of c, for r, in row)- 1 < ©QB+1) 
Therefore, ~ (sum of ¢, for r, in i <C + 1/QB+1)- 


Since all c, and C are positive integers, the above es 
wen in one bin side Setaeastlae WAS waekiad) gc -as desired. ... oO 


Corollary 4.1: The modification of PI removing the minimum spacing euemen iS NP-complete. 


Proof: The same proof is used. Rectangle R has dimensions h = - Bw and w= C(B+1). For each i i,t, 
has dimensions h, = c(B+1) and w, = 1. The NP-hardness part ofthis prof does not require that any 
of the dimensions be integers or that the rectangles be placed so = tate comers are, on integer points, of 


the coordinate system. g.c8ei on a : : 7 < | - 


Lemma 4.1 is presented with spacing required between rectangles to closely mirror the problem 


in circuit layout. The dimension of each componeit ‘Cai ‘be increased: by one: Unit to abeotmt for the | 


Siar sao lagt aegis ere MEE Ree oe gent She tg Cage 
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Fpctangle to be one unit too large in each 


required spacing implicitly. This will cause she cimpumy 
dimension. When spacing is not explicitly required, the. probles remains NP-complete even if the aspect 
ratio of the circumscribing rectangle is bounded. 


Probicm P2-a: Layout with haunded aspect ration end up iatarcomnections 
Given: A set of n rectangles, and a positive umber A. For ISign, each rectangle, f,, has 
dimensions h, and w 


Question: Is there a placement of the rectangics on the plane with a cartesian coordinate system 
‘imposed so that: = 
(i) Each biaiedacy' is parallel to one of the mage system axes; 
(ii) No two rectangles overlap; 
(iii) There is a Fectangle i in the plane which circumacribes the Placed Tectangles, has 
. boundaries paraltel to the axes, ‘i. an | at jee A, “and has aspect ratio (long 
~ gide)/(short side} at moet a, where er’ls @ ratiihél umber not fess than one. The 
boundary. of the circumscribing rectangle is allowed to contain boundaries of placed 


Lenuna 4.2: eter is NP-hard for any a atational number aot less than ons. 


Proof: The proof i is a reduction from bin sacking similar to ‘the proof of Lemma 41. Given G for 
1<ign, C, and B, construct R with w= oes os ve a B. Each r, is of dingastons 
h, = ac(B+1) and w, = 1. The bound on area is A = mee = aca? The aspect ratio o implies 
that the larger side of the desired circumscribing rectangle is pe most (ad) = W. ‘Therefore, assuming R 
is oriented as in the proof of Lemma 4, nae oF the an He te eo ih of “The rest of the 


proof is analogous otha fo Lemma. anise to the rede. boange tse ft be . é 


Corollary 42: If the bound on aspect ratio, a, is allgwed: a8 ap.iaput for problem P2-a, the problom 


remains NP-hard. 
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Preof: The construction in the proof of Lemma 4.2:can be compiited in time polynomial in the length of 
the representation of a; therefore, can bean input: e pa 


Lemmas 4.1 and 4.2 prove that the layout problem » we are neine is NP-complete even in the 


of components, 7 is the routne 


degenerate case when only placereina ‘4s required. Given ‘@ placeme 
Sfobkem “NP-complete? We do not have a proof of a | cecal ‘NP-completeness result for routing. 
However, in the next section, we will pee two ere results ae subproblems encountered 


in channel coming: 
4.3 NP-completeness in Chanac! Routing 


In this section, we will prove that two problems encountered i in a channel routing approach are 
NP-complete. Recall that i in. the channel routing approach the, soating aah divide, into horizontal ek 
vertical streets. Terminals tic along the sidcs-of the bweets.” Each street is niiade up of a sct of parallel 
chharinicls tn the direction Gf thé strict. Bach chiamiel br wide cough to aaa! or the width of a wire 
and the required separation between wires. First the interconnection pattérn of paths through the streets 
is ‘chosen eo routing). Then, within cach street, the =e routes owe path eieagayes ox to the 
street - nen the Cuan -- is & determined est sigma) ‘Tre cal is to minimize the overall 


layout area. see Figure 4.3 ) 


fy 


The set outing problem canbe epee at a raph problem There will be a node in the 


Es wie Nee - 


graph for cach position eras a strect at lsat thee is are terminals Lee: terminals dicty across om 


1 


each ches on a street are sreaieensiiss by one node) There why ato be a ode in the raph for ech 
Sub aa Seal ef haze aig? gk inp Bee 

intersection of streets. Each edge of the ate represents a aitioa of a street between two positions at 

which there are terminals of intersections: For cath het; we wish to fividta dubgraph which is a troe whose’ 

nodes consist of the nodes representing terminals of the nct and nodes representing strect tterscetions. 


This is a Steiner tree problem on the graph. The intersection nodes in the graph are analogous to the 


GZ 
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added nodes in the plane for a Rectilinear or Euclidean Steiner wee problem. However, we do not wish 
to find a minimum length Steiner tree in the graph for cach net. Using a minimum length Steiner tree 
for cach net does not necessarily yield a minimum area layout. ‘The eventual arca of a layout must be 
estimated when the Steiner tree for each net is being chosen. | 

The second phase -- channel miguede = cass ie actual a of ihe layout. We are 
assuming cach strect is of variable width. The number of channels used in a street dircctly corresponds to 
the strect width. Channel assignment detcrmines the actual paths in the plane realizing the 
interconnection. pattern determined by strect joan: The paths are composed of horizontal and vertical 
segments. The path segments within cach street are detcrmined independently except that the segments 
for paths which change streets must be connected -at the intersections. The horizontal segments in a 
fisrizontsl ae (and vertical segments in.a.vertical strect) lie in channels. Fach channel § the region 
between two lines parallel to the strect direction; tic lines. are. spaced so.that there: is: coom for one wire 


width and required separation between wires in any channel. Wire segments perpendicular to the strect 


direction are used to connect segments in channels to each other and to terminals. 


43.1 Channel assignment within : sre. 

The first NP-complcteness-sesukt which we present proves that, with certain restrictions, the 
routing of paths in a strect so that the number of channels used is. sainigiaed is NP-complete. The 
restrictions are that (i) all terminals to be interconnected lie in the street (i.e. there are no street 
intersections to worry about), and (ii) cach set of wires interconnecting one net uses exactly one channel. 


This problem is represented as follows. 


Problem P3: Channel Assignment within a Street 
Given: A line segment, S, containing equally spaced points numbered | through h, and a sct of 
terminals, I’. ‘The line segment represents the strect. Each terminal, t, is an ordered pair (i,b), 


~ where i is a number between | and h indicating the position of the terminal along the strect and 
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b is an element of {0,1} representing the side of the street on-which the terminal fies. For any 
two terminals (i,b,) and (,b,) with b, = b,, [ij] 2 s, where s is a positive integer chosen a priori. 
Integer s. represents: the wise: widths plas equivéd separation. Also giver are a collection of n 
disjoint sets of terminals, N, for 1<i¢n, repemmpnting nets, and a parameter k. 
Question: For each net, Di, let a, be the position of thie terminal of lowest position in N, and 7, 
be the position of the terminal of highest position. Lat [a,4} denute the set of all. points on S_ 
between a, and z, inclusive. Is there a mapping. ch, which assigns cach net a \ number sr OeENeen 1 
“and k inclusive such that for any N, and N 1 LSidsa: vee 
(i) If terminal (x,0) € N; and ae (y, 1) € N, and beyKs, aia Henin Xen; 
(ii) IFfa,z1 byzlisn is not empty, then enchiN) » * th : 


The mapping, ch, bigest the seeenienY | of a wire segrneat jot ponte a, and z, to one of k 


oa] tpiere bp 


channels for each net N, The interval fayn) i is called the ieterval of Bet, Ne. Scamenta, :perpendicular to 
the direction of the strect are assumed to go from each terminal iif: the net'N, to the segment in the 
channel. The restriction that these segments must not overlaps, represented as condition (i) above. 
Af there are no terminals satisfying the hypothesis of condition (i), i.e. condition (i) is the only 
felevant condition, then Problem P3 is the interval coloring problem... The interval coloring problem is: 
given a finite set of intervals on a.line and a positive integer, k, assign a calor {positive integer) to each 
. interval so that no two overlapping intervals have the same color and pe more than. k colors are used. 
Nets define intervals, and “channel” is just another name for “calor”. The. ingerval coloring problem can 
be solved by a polynomial time algorithm [Gav72] [Has71]. The solution to this problem uscs the same 
number of channels as the maximum over all points on §,of the number of nets, whose interval-[a,2] 
intersects the point. Therefore even if we allow wires for ong pet to use more.than one channel, the 
solution found using one channel is optimal. _ 
Without additional restrictions, the channel assignment. problem. stated as Problem P3 is 


NP-complete. 


“a tepres cera cade eth 
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Lenuma 43: Problem P3 is NP-complete. . 


Proof: ln oil is demi Ting Maan ots maping snd cio 
(i) and (ii) are satisfied. Therefore, the problem i i NP: 

We will reduce: the, circular-arc. coloring. petbiem; aa NPcampietc problem (Carl to 
Problem P3 to prose that Pi is Part Theta a ap olorrg problem simila othe interval 


coloring problem except that are ona circle rater than intervals on a Ting are use. Figurg 44 ifhusrates, 


| "The Gireuar Are Coloring Problem ; | : 
Given: Aine eta of ce and pose mer, af 


ee 
oe 


Question: I there an assignment of colors numbered 1 though & to the ares such that any two 
ares which overiap are assigned Uifferent colors? Acs which titefsect only i i crip ae 
not. considered as oveslapping. 


“we may modify any sct of arcs so 


- Since arcs which intersect at endpoints are ndt overtappin 


that no arcs ‘have eridpoints in common {see Figure 4'4); “Tie: actual lerigthi of thé arcs is irrelevant. 


"Therefore, for nates, we can number the endpoints trod 1 through 2 While traversing the circle ini some 


direction. Bach ate will be: represented as‘at ‘orderéd pal, (if), listing ‘the éndpoints of the arc as 
encountered in the traversal of the circle. . Say | 

Chiven an indtancd of the circular are coloring piubilent with fi aiGd 4nd’k colors, we will produce 
an'instance of the channel assignment problem With k-chanhels and (2h +(2K-c-+ 1)) nets, where ¢ is the 
rember of arcs which contain irc (2h,19. Bach néc will Gonitaid exaedly tivo terminals. Tatuitively, we cut 
the circle between points 21 and 1 anid tfc ic ot ‘straight: “Thid Fest in +c intervals, since c ares 


have been cut in two. The n-2c intervals which de not hav SUidissiriey’On' the Clit fine will Become the 


intervals of nets. Considering ‘unty thiese nets, any iegal assighmtnt of channels will be a legal assignment 


of colors to the corresponding arcs and vice versa. 


-=§}- 


Figure 4.4 Circular Arc coloring, 


(Fach color can 


sore heck) 


_ates (a).b}); (a2,b9); (a3.b3):(ag,b4) traversing clockwise... 


a 4.5 Construction of a channel coleman probe . 


gives.a circular age 
Given: IS oe Lgl oa 
Key: : 
; g-— ‘represents the interval of a net. 
cut Arrows point to the side of the strect containing 
the terminal. 
Cares 
Construct for k colors: .] +\ 
: {| sipe 0 | 2 4 
ee ae | Ni, 
‘ : . . i. 
3 ae . Noa . intecvals ne : 
wie ar Gee ase ee _ 
: oan Gar) 
additonal arcs | | i ' | additional 


Nip wise, jebre\ 
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There are 2c intervals with endpoints on the cut line -- two intervals Per are. if we can ingure 


that both intervals of each pair are eee the same channel, then each- channel assignment will! 


correspond to an arc coloring. We do this by extending the intervals beyond the cut line and adding nets 
which force the pairs to be amcigned to the same channel. Figure 4.5 gives ine construction. The formal 
definition is given below. abe points on which terminals lic will be numbered from -o(k-{e-1)/2)+1 to 
2n+c(k-(c-1)/2) eather than from 1 to 2n+c(2k-c# 1). - 
For each arc (a,z) which does not contain (2n,1), there will be a net {(a,0),(z,1)}. Order the arcs 
which do contain (2n,1) by increasing starting pan. For ue i" such arc, there will be a ad nets -- 
half have terminals within the negatively puitiiberkd’ poiiea: and ‘half have serena within the points 


numbered above 2n. For arc (a, z,), 1$7, <a, $2n, the i arc containing (2n,1), the nets are defined as 


follows. For 1<Sigc, let : ee 
i-l 
sum(i) = © (k-j) = i (k-4(i-1)) 
j=0 
= -sum(c) + sum(-1) 


Py = 2n + sum(c)- sutfi(i-l) 
Ng = (0) +1.0),@,D} and 
Nj = fe +j,)), ; +)+1,0)} for Isjgki; 
= {(2,1), (p +.0)} and | | 
= {07 4.0.@ T4+1)} for 1sigk-be > oe, — Dee ee 
For each i aire 1 and c, the’ application of condition (i) in the definition of problem P3 


results in two chains of k-i inequalities for the scto€ nets convespinding 6 to the i arc containing (2n,1) - 


one chain of incqualities ch(N;;) < ch(N” and pone of inequalities ch(N, f)<ch(N* | +p for 


ij+D 


O<j<k-i-l. Therefore, nets N_ and Nii must be assigned channel! 1 if no more than k channels are to 


Co tyes 2 Aare EpE Rp ES ocr Tite dg fy epi see Eg Oe oct Be are ee 2 
Senge s Re See ee © EMSS eit Ravi ae tae ak ig See ia mae 3 
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be used. Nets Ng and N, 7 can be assigned channel 1 or? without violating condition (i), but their 
intervals overlap thove for nets N,z and N, $..respeetively. “Therefore channel 2 is the only choice for 
nets Nyq and.N,%. Proceeding in this way, we ate that nets N,” and'N; t- must be assigned channel i if 
no more thar k channels are to be used. 

Given any channel assignment for the ‘complete set of nets, we color the arcs which do not 
contain (2n,1 ) by the same riumbered color as'the correspotiding channel. For the i" arc which contains 
(2n,1), we use the number of the channel assigned wo iets Ni and N,t. For afy pair of arcs which 
overlap, there is at least one.corrcsponding pair of nets‘wheest intervals overlap. Thetefore, any legal 
channel assignment corresponds to a tegal coloring using the sme number of colots as channel. 

‘ Given a coloring of the arcs using at most k cdlots, we can assign tic ricts to channels as follows. 
Permute the colors so that the 1" arc containing ‘atc (28,1) is assigned the # color. . Assign nets Njj and 
NT to channel i+j, for 1Si<c and O<j-<k+i. ‘The remiaining nets are assigned the same numbered 
channel as the color of the corresponding ‘arc. Between pemts I and 2n, intervals overlap if and only if 


ij 
However, ch(N,) = i<p+q = eh(N,), $0 nO Nj and N,, Which ‘overlap’ are assigned the same 


their corresponding arcs overlap. _ Elsewhere, no Ni, and No for Kp overlap except Ni and Nog 
channel. An analogous argument shows that the Nj are properly assigned to channcls. Therefore, for 
“each coloring, there isa legal channel assignment using the same number of channels ascoloms. = 
Since the construction uses only nets with two terminals, we have: 
Corollary 4.3: Probicm P3 restricted to nets containing exactly twe terminals is NP-complete. 
4.3.2 Channel assignment with intersections. 


The second problem which we prove NP-complete deats strictly with’ the ordering of paths in 


intersections so that the resulting street widths minimize the area. This problem is somewhat similar to 


the channel assigament problem in that if two ‘paths are-approachiag an intersection from the same 
direction in a street, and one path needs to go left.at the intersection and the other nceds to: go right, then 
they cannot share the same channel in the. new street unless they.are in the propor order whem they reach 
the intersection. The problem will be modcled using a graph to represent the stroets.. Subgraphs which 
are trees will be used to represent the interconnection pattems resulting from street routing. 

Let the strect graph, S, be some subset of a two dimessiosal grid. graph. Each: node in 5 is 
labeled with integer cpordinates (p,q). Fach edge is either horizontak:i.e, between: nodes (p,q) and — 
(p+ 1,q), or vertical, i.e. botwcen nodes (p,q). and (p,q+ 1), The:graphr §.is partitioned into streets. Each 
street is a path jn. S using only vertical or only horizantal eiges: A horizontal strcet goes between nodes 
(ik) and (j,k) for some k and i<j; a vertical street goes between nodes (k,i) and (kj) for some k and i<j. A 
node represents an intersection of two or more streets, An. edge represents the-partion of a street 
between two intersections. ‘The interconnection pase for cach act-is represented by a tree in S. Each. 
tree, T, will be called a net tree. We would like to.assign cach occurrence of an edge in a-net-tres w 8 
channel. Let ch be a mapping from each occurrence of-an exige in’ a net tree to a positive integer. The. 
integer indicates the number of the chaancl containing that cdge in the. street to.which the edge belongs. 
We require: 3 

(1) If edge e of S is in distinct trees Tyand Ty then, the eocurtence of ein. T, is 
assigned to a different channel than the occurrence of e in T.. (No overlapping wires.) 
(2) Ife, and sc niece cdg canto a eee 

S, then ch(¢,).= ch{e,). (A connection path cannot change-channels within.a strect.) 

(3) Suppose horizontal edges ¢, = ((p-1,q),(p,q)) and ¢, = (pad q)), for some 

p and q, belong to a strect, s. Furthermore, suppose ¢, and, helong to-distinct-het trees I, and. 

T,, respectively. If ch(e,) = chfe,), then; 

(i)¢, is not in T, and.e, is not in T,, (This follows from (1) and Q).above.) - 
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(ii) Suppose that vertical edges ((p.q)4p.0-+))and tpsa-1)4.0)) are in one.strest; and each of 
T, and T, contains at least one of the edges. Then chi(s,) <ch(x,), where x, represents. any 
occurrence in T, of either of the edges and x, represents ay wae in T, of either of the 
edges. This condition insures that the wire segment repeesented by horizontal edge 
((p-1,q).(p.q)) does not overlap the wire segment: sepesenad by horizontal edge 

((p.q)(p + 1.9). 
(4) Analogous to (3) bu for vertical edees (.4-14(P.9) and (Pah(P.9+ 1) in T, and 
Ty T,, respectively. If x, and x, are occurrenges of hotizonial edges ((p-1,q),(p.q)) and/or 


((p.q){p+ 1.q))i in T) and T,, respectively, then en p< chix). 


We want to find an assignment of edges in — to channels so that the resulting overall area is 
minimized. The assignment is called the intersection channel assignment since the intersections induce 
the constraints on the assignment of channels within cach street. Area will be measured as follows. For a 
given channel assignment, ch, let width(ch,q) be the sum over all vertical streets containing a node (i,q) 
for some i of the number of channels ‘ised in the street; let height(ch,p) be the sum over all horizontal 
streets containing a.node (pj) for game j, of the number of channels used in the street. Let width(ch) be 
the maximum of vec) over alt aieeot q appearing. as the secosd peorcinme of some node in the 
street graph. S. If the sxaimum is aera, then widtich) i one. bet hoi be the maximum of 
height(ch,p):over all integers p which appear as’ the first coordinate of some nol in S: If the maxiraum is 
zero, then height) is‘one. ‘Then, atca(ch) is defined a as the product of witty and height(ch). 

Using the representation presented above, we have the following probicm, illustrated in 
Figureag. 2» © | 

Problem P4: The Intersection Channel Assignment Problem 


Given: A street graph, S, partitioned into streets; a collection of net trees, T; and a positive 
integer, A. 


Toetac cafe ane: peel tee te RRS Be rie: spe aa SNE odd ots Tb Lied) Ser ey eee eT RR Se 
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Question; Is there an assignment, ch, of the ocgurrenses pf-edges in. net trees.to channels in 
streets which satisfies conditions (1) through (4) above such that area(ch) < A. 


Lemma 4.4: The intersection tanned sereemient peotiar ia NP-complete. 


Proof: Conditions (1) hoe (4) and the area of an et can be checked by a ‘nondcterministic 
Turing machine in polynomial time. Therefore, the problem i is in NPL _ 


“We show that the problemi in NP-hard by reducing -satsabilly (Gar to it: 


4 ts re ae 


Given: A boolean expression composed 0 of the conjunction of k clauses, C, for sik. Each 
ene is the disjunction of three distinct’ literals, where a literal is a boule variable or its 


complement, ie. c; = ee where Vy is x or —tx for some variable x. 


. ‘Question: Is there an » amigninent of truth Soa to hike boolean variables such that the 
© expressionis satisfied? =n 
Gira; an instance of the 3- -satisfiability Sobiew with k ae acid Vv ¥ variables, we e will construct 
an instance of the intersection channel aikigninent problem with 2k +1 horizontal streets, w+ vertical 
streets and A = 2vk(3v+3). Let § be the (2v-+ 1) by (2k) grid grape with aodes numbered from 0,0) 
_ through (2v,2k). Each path from (0,i) to (2y,i) is a horizontal street, for-any i between 0 and 2k; each path 
from (i,0) to (i,2k) is a vertical strect, for any i between 0 and 2v. Certain streets are associated with the 
clauscs and variables of the boolean expression 2s follows: 
(a): With, each: clause c., 1<i<k, associate the-hosizontal street from (,2i-1) to (2v,2i-1) and 
name it C,. The remaining horizontal strcets are yanamed. 
(b) For each variable x, 1$)<Vv, associate the: vectical, strect from (j-1,0) to G-1,2k), named X,, 
and the vertical strect from (2v-j + 1,0) to(2v-j-+ 12k), named Xp” 


(c) The one remaining, vertical street, that from (y,0) to (v,2k), is named street M. 
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We construct two net trees for cach variable. For variable X, the first net tree: i? contains the 
following edges: ee ep see 


(1) In street C,, for each i from 1 to’k: all edges onthe path from (j-1,2i1) to 
(2v-j + 1,2i-1); 

(2) In street X): the edge from (1,0) to G- . mee fe all even a peween: l and el 
inclusive, the two bias on the path from (j- Lae to G- -] 21+. Also, if ki is even, the edge 
from (j-1,2k-1) to (-128); 

(3) In street Xig: for all odd i between 1 and k-1 inclusive, the tie expt cn the path 
from Orit, 2i-1) to Qvj+1, +d. Ais, ifk is odd, ihe cage between Ovj+1, 2k-t) and 
Qv-j+1, 2k);. 

(4) In street M: if variable x3 appears in — c, uncomplemceted, then if i i even, the 
edge between (v, 2-1) and (v,2i); if i is odd, the pe saiucte (v,2i-2)-and (v,2i-1). It, res 
in ¢; complemented, en et is even, tie ode besweves (v-2ha) ana (v2ET): iti Sold the edge 


ences (v,2i-1) and (v, %. 


The second net tree for variable x;, T,, contains the following cdges: 


py 
(In street’ C,, for cach i frofi 1 tok: ‘all edges ‘onthe ‘path: from (j-1,2i+41) to 


(2v-j+ 1,2i-1), 


"from (j-1, 2-1) 40 (2, 241). Also; ite bk the edge Between G E281 a G12 
(3) In street Xj: the odige froin (2v-4-+1,8) (244 1) and forall even {between 1 
_and k-1 inclusive, the two edges bi the:path Mein (24-4228) 40 (G4 1341) Alby, if k in 


even, the edge from (2v-j+3,2k-1) (2v-}+ 1,28); 


(4) Tit strdet'M1 if vatiable 1) appears in claust'¢; d, then if i is odd, the 


edge between (v,2i-1) and (v,2i); if iis even, the edge between (v,2i-2) and (v,2i-1). If X, appears 


-61- 


in c, complemented, then if i is odd, the edge between. (v,2i-2) and.(v,2i-1).if.i is even, the edge 


between (v,2i-1) and (v,2i). 


The following observations are useful. For any coordinate p; cail~an edge. fram p,2i-1) to 
(p.2i-2) an edge down from street C; call an edge from (p,2i-1) to (p,2i) an edge up from C,, Each net tree 
T; coatains a path which begins, with an edge in, street X,down from:street Cy. The pagh goes through C, 
to street Xp goes up from Cc, to street C,, and through C, to street X;. The path continues snaking 
through the horizontal streets associated with. che. clauses, using strect. X; : go from street C, to street 
Ci, 1 When i is even, and eae street” Xig' to 26 5 Fro Gwe Can “when i is odd. Each net tree "p also 
contains a path which snakes throtgh the C, but in, she opposite dition: it begins with an edge in x, 
down from Cu uses mist Xo to yehanes from C, w c “when i is even: and uses street x; to change 
when i is: odd. For : any net tree 7 or T. pif ‘ ao ‘incomiplemented in C, thd ne in street M 
| intersecting street C; é in the same direction fom c as the edge of the tree in street xX; if x, appears 
complemented, the edge in street M is in the sae drocin from ane C, as the sae ok the tree ia atreet 
Xig (opposite to that in street X;). Figure 47 gives an example of the'constrution 
For each i between 1 and k, each net‘treq confains the edges between horizontal positions v-] 
and v+1 in street C.. Therefore, ae ed for-eacl net tree, giving a height of 
2vk for any channel assignment. Any channel assignment which gives area at most avkGv+ 3) must give 
width at most 3v +3. | . ; 
At each intersection of street M with astreet C, there are * net trees ches contain edges of M 
incident on the node for the intersection -- one ‘aie of net trees, T; and Tp for each variable x; appearing 
. (We assume that each variable occurs at most ofece in a clause. Note that ls is always 
satisfied.) Half of the cdges are up from street tG and half are down. Therefore, at last three channels 


are required for M. Three chanpnels will be used in M exactly when the boolean expression of interest is 


satisfiable. 


aie gage Rage ees tt Se eB ie a me ger ee eee 
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“Figure 4.7 Construction of the proof of Lemma 44. 
" Espression: (x1 V x2 V 1x3) & (“xy V x9 V3) 
Cy C 


Truth assignment: ky = Kp = X3 = true. 


‘street 
Cy 


street —s —_ _ — 
3 T 
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If the width ia required to be at most 3v-+3;:and street M contributes at least three, then the 
width contributed collectively by strects X, sad Xi for 1X) <v, lust be’ at most Sv. Consider the 
intersections of streets X, and Ky with reet'C, Roe iota! ‘The iige fot (24-2) to (125) in X, 
belongs to T,, and the edge from (j-1,2h1) to (1,2) in béfongs oT). Therefore, if only one channel is 
used in X,, the channel in C, containing the edges of F; misst-frive & lower number than thechannel in C, 
containing the edges of Ty (y condition (4) given for'chatide? sshiginnienty: “However, in Xjq, the edge 
from (2v-j-+ 1,2+2) to (2v-j+ LIEV) Belongs 00 'F, anid-tie -etigt Hom (244 281) to (2v-] + 1.21) betongs 
tol. If only one channet is weed iti X y the channel tn €, containing edges of F., must have the lower 
number. Therefore, only one of xX, and Xo can contain exactly one chante! If is even, the edges which 
belong to T, and TF}, afe just interchanged frorh that above; we reach the same bonchaion: Therefore, 
cach pair of streets, X/ and Xi contribute: at feast thret chaitneli to the width. ‘However, if collectively at 
most 3v channcls are contributed by stheets X; and X.,, HS}-<v, at most three channels must be used by 
each pair: Note that T; and T, "fit.together" so that-théy can‘always share‘one channel in x; or one 
channel in Xp, In the other of X; atid’ Xp, edges of Tate assigned to one ‘channel and edges of Tyg are 
assigned to a different chanael. We will associate acchanntl assignment in which X; has one channel with 
a value assignment in which x, has value “true”; we wil alsdetate a channel assignment in whictr Xig has 
one channet?-with a‘value.assignment of “false” for x). . 

Given an aasigriment of truth’ valites'to the variables x, such that the boolean expression is 
satisfied, the corresponding channel assignment giving width 3v-+3 is as follows. For each j between 1 
and v, if x; is true, one channel is used in X;; if x, is falee,-one channct is used:in Xo. Street M contains 
three channels. We first determine what trecs wif share channels at each intersection of M and some C,. 
Note that only-adjacent edges in a tree are-required to dse the sanie channel: oiiges of a ti¢e which are in 

- the same strect but not part of the same path in the street are si reQjitited40 use the same channel. Let x 


be a variable whose occurrence in ¢, is truc under the given assignment of truth values. The edges of T, 


and T, in M at C; can share achanael. This follows from the fact that if-x, 4 uncomplemented in c,, the 
edges are.in the same direction as those-in'X, and the adges in X; share a channel; if x, is complemented, 
the edges are in the same direction as thee in Xpand the edges i Xq havea channel, However tbe 
pairs of edges in M at C, will not be able to share 4 chagnel if the corresponding literal inc, is falec. 
Therefore, let the edge for x, in M sine froma G, share with the edge-up-feom C, for a.sccond variable, 
inc, Let the edge down ‘e Xp share with the edge up for, the third: variable x_. Finally, the edge: down 
for 2, shares with.the edge up fora, The new constraints am de arsignramt of chanects it C,induced:by 

this sharing are consistent with dhe conatraiis from any sharing a the intersections af C, with: XX, Xp 
Xp Xoq and Xoo (See Figure 4.8.) 

Te gine hres btwn he Cae cad ay netteee. Thesefare, 
the only condition relevant at.the intersections of M. with these strepis.is.chat a.path ina nettree which — 
goes through an intersection caanot change cheanelx, There will be @.path in.a net tree through such:an - 
intersection if the same literal appears in clauses c, ang.c, , », cousingan edgeup from: C, and down from 
C41 in one of T, and Ty The edges. arc adjacent af (v2). ‘The: agsignment to chanacls ia M'can-be 
sirct at horinantal potion 0 arbitrary amign 
cach edge in a net tree incident onthe intersaation.to.a diferent channel aranng. channels. 3: through 3. 
These edges are edges down from street C,. At C,, the.edges dawn have heen sesigned: asnign the egos 
up so they share with the edgcs down, as previously detesminad,. Now consider the: interectiag of M with 
the street between .C, and.C,.. The, channcls. for, any tree. pathestenen C,:ho-Cy hve, been. deseemnined at 
C,; the remaining edges down from.G, can be assigmed arbitrarily, e:this way..orters.clown: from each:Gy 
can be assigned channels at the intersection of M. with the hecizental steeet-before C s0 shat. paths do not 
change channels, Edges up are assigned at the, iatersoctiop.of. M-with (, ap. that.chanacls are-shared: 
properly. Thus, if there is any assignment of beplcan variables satisfying the hoolean expression, there is : 
achanncl assignment, ch, with width(ch) = 3v-+3 and arpalch) = 2vk(3v+-2)< A as desired. 


determined as follows. .At the intersection of M. with: 


eR ce a 


T T | Tye 


‘dap 
down ” represcnts the tree among T: aks sierra for variable: 
ty - ~ yequal to}, p, org. eee, a 


“pup | ee ae 


“Ty 


Ty 


U —» V indicates the channel i in C, coming edges oft net tree Ur must fave a ore number 
thap-that containing edgcs of net wee V: 


—€—9—-_ indicates that the order.could be either way, 
For any choice of direction for the two —~“——"y—. edges, no cycle is created. Therefore the . 


up ” represents the tec contwning ~ edge up. 


- 66 s 
Given a channel assignment with area(ch)< A, we must ‘stow ‘liow to construct:a boolean 


assignment which satisfies the expression. We have shown that any such chanrfel assignment induces a 


truth assignment by choosing x; true if X; contains, nly: on¢ channel and x; false. if Xp tontains only one 
“channel. It remains to show that this assignment satisfies-the expression. We know that exactly three 


~ channels are used in M. At each intersection of M with some C.. edges from two different net trees must 
: “4 q . on ‘ 


share cach channel. ‘Among the trees containing edges in N at the inversection with C, consider the tree 
whose edges in C, are in the aot numbered channel. This tree must contain ‘é cdge‘in M down from 
C;. Suppose the tree is for variable x,. The channels of €y signed to edges in, and Tip are in the 
proper order for the edges of T, and T. gp a Caner ance it 9, sore ance nee 
then the edges in M at C. for T, and T\p are in the same direction as those in X, at C, It must be street x 
which contains only one channel. Therefore x; is true, ‘and c j iS s satisfied by era x; it x appears 
complemented in in c,, then the edges in M at tC, for T 1 and Ty ar aren se same divestion atone i: Kyy. It 
must be street Xig which contains only one chaanrel. The value of x, is falar, and c, is satisfied by literal 
—x,. A fitera! satisfying each clause can'‘be idcitified by looking af the intersection, of M with, the street 


for the clause. Therefore, there is an assignment of truth values to the variables such that the boolean 


. expression is satisfied. Oo 


In the construction used to prove Lemma 4.4, height(ch) is the same for any channel assignment. 


Therefore, we conclude: 


Corollary 4.4: The modified Intersection Channel Assignment Problem in which one desires a channel 
assignment, ch, such that width(ch) < D for some given integer, D, (or such that height(ch) < D) is 


NP-complete. 


Lemma 4.4 shows that even ignoring terminals, assigning channels is a difficult problem. In the 


next section we again consider channel assignment within one street. A heuristic algorithm is analyzed. 
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An example is constructed to show thatthe. algerithin-can be snadeto do arbitrarily badly with respect to 


the optimal solution. 


Seen Mn ER RARE ge Te . wn 


In this section, we discuss heuristics for the chaneel assignment problem within a street. We 
begin with a general discussion of yoru restrictions whose removal changes the optimal channel 
assignment. Then we present a heuristic algorithm and its analysis for Problem P3, the version of the 
problem proven NP-complete. a eS 

The ordering on ch(N,), required by. condition. (0). of the, statement of Problem P3, caa be 
represented using a directed graph, which we will call the auuea (dain bias will Pee one ‘pode for 
each net. If ch(N. ) < ch(N;) i is required et ini at ssganestof Problem P3, shen there is 
an edge directed from the node for N, to the node for N, It is possitile for the sraph to. be cyclic. In this 
case, there is no channel assignment for the problem a stated (P3),! ‘If each net can use more than one 
channel -- by using wire segments in the direction Sercvesidivalal’ a0 the street’ direction ta connect 
segments in different channels -- then a channel mpepeaes my exist. The see ‘a kame to 
the street direction are called. joss. ita used for different pat like a any other ¥ wire segments in the same 
direction for ferent ne, must be spared bythe minima spacing : 

Even when jogs are allowed at any saint along eae: the channetsxgnment problem may 
not have a solution. Figure 4.9 gives an example. However, ifdogs pre allowed anywhere between points 
on a street rather than only at the points, then the channel aseignment problem is solvable in polynomial 
time [Ka79]. The model of a street used in [Ka79} differs slightly feom the formulation ‘used ia P3. For 
this discussion, we only need note feacis [Ka79}terminals on epposite sides of the street arc cither at the 
same point along the street or are at points separated by at least minimum spacing. 


Allowing jogs between points implies that an arbitrarily large number of scgments 
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Figure 4.9 Channel assignment problem with no solation:eves if jogs sreallowed. 


] 2 SIDE 0 3 . 4 


i 


f 2 = smpE1 3 <:. 


Nets: {(1,0), 3,D} = Ny (2.0), (4,D} = N3 See aaa 


{(1,1), 3.0} = Ny {(2,1), 4.0)} = Ng 


“IPNy or No inses a jog at 2 to change channels, neither Ny not'Ny cam use a jog at 3: 


Figure 4.10 Inserting jogs between points to effectively eliminate terminals across from one another. 
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perpendicular to the street direction may be inserted between any two terminals. The length of the street 
~ the dimension of the street in the direction ‘of the street -- must be variable, rather than just the street 
width being variable. The algorithm in [Ka79} minimizes only strect width, not the area of the strect. Let 
the maximum overlap of a set of nets be the maximum over ail points on the strect of the number of nets 
whose intervals overlap at the point. An optimal solution under the model in [Ka79] requires at most one 
more channel! than the maximum overiap of the nets. Recall that if there are no constraints due to 
condition (i), the channel assignment problem P3 is sofvable in polynomial time. Assume, as in [Ka79], 
that terminals are either at the same position along the street or have sufficient space between them. By 
using one channel and inserting a jog between every. pial ‘of consecutive termirial positions, we can 
effectively separate the terminals on opposite sides of the street sd that tiere are m0 constraints under 
condition (i): (See Figure 4.10.) We at most double the ‘streot length. Then the number of channels 
needed (excluding the channel we used to'create the effect of suitably spaced terminals) is the same as the 
maximum overlap of the nets. This implies that one never nedd to étigtheet the niveet by more than 
double to get a channel assignment which uses within one channel i i indie 

Since we are thinking of a street as running along the boundary of a component and we do not 
wish to think of each component as expandable," allowing streets to lengthen is not satisfactory. An 
alternative is to use jogs in a street intersection. When a jog'is needed, a segment which goes to the end of 
the street and out into the street intersection is used. ‘Ih the intersection, the path can use a segment (the 
jog) in the perpendicular street to change to another chaiinel in the original street, and reenter the region 
of the street it was previously in. This type of solution’ requires that perpendicular streets be'avaitable and 
that nets are allowed to have wire segments in two channels of a street at the same position along the 


street. Figure 4.11a illustrates. Such a routing not only affects the width of the strect containing the 


1. In fact, some components are expandable [Jo79}. 


Rohe. Sty See a ee NCS SE PRE See SAEs SER te 
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terminals being interconnected, but affects the width of the perpendicular street as well. Using this 
scheme, any collection of nets consisting of terminals in one strect can be interconnected. Figure 4.116 
shows the general connection pattern. | 

eee eee original channel assignment problem, problem P3. The proof of 
NP-completeness of problem P3 relics on the fact that jogs are.not allowed. .Figure 4,12 shows that even 
for instances of the problem derived from circular arc coloring, allowing jags. reduces the number of 
channels accded.. We do not have a proof that she channel assignment. problem with jogs in NP-hard,. 
although the author believes that this is the case. 

We now present.a.simple heuristic algorithm for problem P3 gad analyze the quality of the 
solutions it produces.! This algorithm is the basic algorithm used ia [Deu76] when jogs are not allowed. 
We will assume that the constraint graph produced by applying. condition (i) does not have any cycles.s0. | 
that a solution without jogs does exist, We first define the Jevel of a node in the constraint, graph: | 

(1) All nodes which have no edges into them are on level 1. ; 

(2) Assuming that the nodes.on levels 1 through k-1 are defined for k>1, a node is on level k if 

__ all edges entering it are from nodes on lower levels and it is not on a lowes level. 

Since the constraint graph is acyclic, the levels are well-defined. They caa be computed by starting with 
level 1 nodes and following edges. We will say that a net is at level. if its node is at level k. The nodes 
from which there are edges to a given pode are called. predecessors of the given node, and the. 
corresponding nets are called predecessorg of the given net, | 

Order the nets in-increasing order by the positina, ,of their terminals of lowest position, If 
there are no constraints duc to condition (3), then the following algorithm Snds the eptimal.solution: 


1. This analysis is part of joint research with Errol Lloyd. _ 


aT 
Figure 4.11 Probleme requiting jogs. 
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One u-shaped path per pet. Top terminals connect to the top half of the u;. bottom terminals connect 
to the bottom half © 
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Figure 4.12: Use of jogs when not necessary. 
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_ Lewto High Fill [Has71} 
Repeat for each channel, beginning we chasnel 1 and eta the channel number py one 
for each new chhone until al nets are asigned: . | 
Choose the first iaiagnad net under the above ordering. Assign it to a new channel. 
Repeat until the channel i is full: : 
Choose the first unasigned net ot whose terminal of lowest position is at a 


higher position than the corenttal of highod position in the last net to be 


_Bssigned to.the channel. Assign this net to the channel. 


Claim {Has7]}:! When there ; are no constraint due to condition es the above malas uses exactly the 


miner of channels as the 1 maximum m vel of the» nets, 


Proof: Look at.the lowest point, ay in the interval of the first. net, N, assigned to the last channel: Each 
lower numbered chaanel, aust uatalaa tek whose interval contains point a, Hf there were a channel, k, 
which did not contain such a net, then the net should haye-been placed in this channel. This follows 
because, in the net ordering, N is before the neta, if any, which were placed in channel k after point ap. 
"Therefore, all Jower numbered channels are assigned nets. whose intervals contain point ay, ‘The overlap 


at this point is equal to the number of channelsused, | oO 


Algorithm "Low to High Fill" be the basis of the algorithm to aay nets to channels when 
constraints are ‘present, ‘Channels are again ‘filled vexing with issued L At any point during 
execution of the algorithm, let the set, A, of available nets contain those nets which are ynenigned and all 


of whose predecessors are zie “Low to ee itr" is odie to choose 2 nets to be assigned only 


1. In [Has71], Hashimoto and Stevens use a different approach to prove the algorithm correct. They 
prove a differcnt but equivalent claim. “The proof given: here:is-based onthe proof uscd by Gavril 
[Gav72] to prove that his ia abi for chordal add eiihes is correct. 
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from set A (the nets are ordered as before) and to ane A after every assignment. 


Lemma 4.5: For an instance I of Problem P3, leben the maximum over al points onthe set ofthe 
meas of nets at level k whose interval overlaps the point. Let h me weil Bgnest level. Then 

(number of channels),,,( 1)/(number of channels),,,(1) : s mith, a) Ss ni 

where the ago! is the modified Low ¢ to re Fill algorithm using set A, ave(d,) i is the average of me 


d, for TSESH, and n is the numberof nes in instance L 


Proof: We know that (number of channels), (E) > maximum of the 4, and'h. 
We prove that (number of channels),,.(1) < a (4). Let ae peas the set of channels containing ; a 
net at level 1. For k>l, let C, denote te set of channels with higher numbers than the lst channel 
containing a net at level k-] and which contain nets at level k. The number of channels used by the 
algorithm is the sum-of the {C,| over all levels. Fora particular k, C, may beempty. IFC, is not empty, 
look at the first net at level k placed in the highest numbered channel in C, : denote thisnet N,. Let'a, be 
the lowest point in its interval. If any-channet in C, id not aatignied i net whose intérvat contains a,, then 
N, should have been placed in this channel. This follows because NY; must Have’ been in A when thie 
channel was assignod ~ all level k-1 nets are placts before:chatincls-in C, are assigned — and N, precedes 
in the net ordering whatever net, if any, was assigned to the channe after polit’. Therefore, IC,1<d.. 
we have: 
(number of han (DM DAoumber of fac 5 < Marlatt) . 
s maxthavad) K  minthavatd, ))/maa(h.max(d,)). $ min(havetd,) 


The numberof et at evel is at east hase 7 . 


2-14) = huaveld,) <a and * inteasata,) <0! meee = 


If the above bound onthe wor ese performance 6th ape i curt, te igo . 


can do very badly. In fact, the bound is the same as would be obtained for an algorithm which applic 
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"Low to High Fill" on nets of a fixed level, one level at a time; requiring afl nets on level k (1<k<h) to 


be assigned to lower numbered channels than those of level k+1. Nonetheless, our bound is tight. 


Lemma 4.6: There is an instance of Problem P3 for which modified Lew. to High Fill using set A 


produces a channel assignment whose ratio to the optimal selution is a(n) for avg(d,) = h = a’, 


Proof: Figure 4.13 gives the instance. There are nf” groups of nets, each of which forms an n# level 
chain. For all level numbers, k, between 1 and h inclusive, the intervals of all nets on level k overlap, 
giving d, = n*®, The algorithm, assigns only one. act to. each channel. However, there is a. channel 
assignment which does not assign nets to channels in: adew to high order. Among nets on a level,. some 
nets which are later in the net ordering are assigned. ta lower suaabered channels than nets which come 
| before them in the ordering. ‘The nets whith camé first fn the ordering céin share ‘duaiiels with nets 
which are on higher levels than they are, but whase predecessors have already been eiapied Refer to the 


figure for details. oO 


The analysis presented above illuminates the fact that an algorithm used in practice is capable of 
finding very ‘bad solutions: It forms a point of comparison for more "clever" or more complicated 
algorithms. The “low to high fill” method can be used as the basis of an exponential time algorithm to 
find optimal channel assignments by trying all possible choices for each assignment rather. than taking the 
first net in an ordering [Ker73}]. 

In the next chapter, we present an algorithm for a special case‘of channel feiithig which is not 
NP-complete. The problem is to route interconnections among two — nets whose terminals lie on the 
outside of one rectangle. This problem is reminiscent of the circular arc coloring problem, which is 


NP-complete. However, the paths around the outside of the rectangle are allowed to change “color” as 


they go around corners, i.c. the order of paths need not be the same on adjacent sides of the rectangle. . 


This order can change because horizontal and vertical wire segments are allowed to cross. Once the paths 


-16- 
Fie. Vet ec or Low High a, epee 6 tine Ee, f 
no oA ma nt inten 
; 6 ry GME wee Pe geil bal cue is an. ear 
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for Sj <n are on level k. 
ela ty Fill ensigns channels as shown above, but optimatselution is: on 
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have been routed through the four streets surrounding the rectangle, the channel assignment is four 
instances of the interval coloring problem. In the next chapter, we show how to do the strect routing 


optimally. 


Chaster 5: An Algorithm for Ronting Terminals on a: Rectangular Componcat 


5.1 Preliminaries 
We will now present an algorithm which finds an optimal solution to the following routing 


problem in polynomial time. Figure 5.1 presents an example. 


‘The One Rectangle Routing Problem 

Given: a rectangular component with terminals around its outside edge. ‘Tcrminals lie on 
positions which have at least unit spacing along the edge of the rectangle. ‘The unit spacing 
represents the width of a wire plus the minimum spacing between wires. A list of nets, cach 
containing a pair of terminals which must be connected, is given. Each terminal belongs to 


exactly one net. 


Find: An optimal routing of the wires between pairs. Paths must be composcd of line segments 


which are parallel to some side of the rectangic.. Distinct paths may cross at right angles; 


however, paralicl segments belonging to distinct paths must be separated by the minimum 
spacing. (We are assuming that there are two layers for interconnection. One layer is used for 
the line segments in one direction, and the other layer is used for line segments in the second 
direction.) All paths must lie outside the rectangular arca of the componcnt. An optimal 
routing is one which minimizes the area of the smallest rectangle which circumscribes the 


component and all routing paths. The sides of the circumscribing rectangle must be parallel to 


those of the component. 


Place the rectangular component on a cartesian coordinate system so that its sides are paralicl to 
the axes. Arbitrarily choose one axis direction to be horizontal and one vertical. Label the horizontal 
sides of the rectangle as top and bottom, the vertical sides as left and right. 

Fach pair of terminals can be connected either by a path which pepdicois from one of the 
terminals or by a path which goes counterclockwise from the terminal. ‘The directions of the connections 
determine a sct of intervals which path segments will use along each side. Once the intervals to be used 


by each path along a side arc determined, minimizing the length added out from the side is a channel 
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Figure §.1: Example of the routing problem for one roctangalar component. 


Terminals are represented as points on the rectangle boundary. .,._. 

Pairs of terminals with the same number must be connected. 

An optimal solution is shown. ous ane ; 
The height addet? to the top by this solution is 2 units: the Height added to the bottom is 3 units. 
The width added to the left side is 3 units; the width added to the right side is 2 units. 
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assignment problem. Vertical and horizontal —— bse oe too ‘one ‘path cs can 1 be extended beyond 


a comer as far as necded to meet Cech calles without snore. any sites lah Therefore, the 


channels in the street along one side of the component can ‘be nigel independent of the assignrnents- 


within other streets. Within each street, there are no constraints de sto terminals across from one another 
(condition (i): of problem P3 in Chapter-d), We have’ four instance of the interval coloring problem ~ 
one for cach side. ‘Therefore, the length added out from each side is equal to the maxiriuan ovettapiof the 
intervals on that side, and it suffices to use paths which only change direction to round a coraer of the 
componcnt or to connect toa terminal. | pH 

We must determine. the ditoctian of each connecting path so that the resulting, area will be 
minimized. We consider two types of connections, Pairs of terminals on the same side or adjacent sides 
of the component are called local connections. _For these, it suffices ws chose the direction which goes 
around the fewest sides. A path which ais the lng w way adds. at keast 0 one e unit in cach dimension to the 
circumscribing rectangic. A path which got the short Way Canifiat add more fhan this. Therefore it is 
never better to go the long way around. 


The second type of connestion contains those ‘hich ao (ram the 8 sew the ght side or the 
top to the bottom. The choice of direction for top botidan: connections is (beeen er/ of the direction of 
left-right sbanecione and vice versa, since regardless of the direction used fora top-bottom connection, 
one unit is added to the horizontal dimension of the circumscribing rectangle. Therefore, we have two 


instances of the following problem: 


Top-Bottom Routing Problem 
Given: Terminals on a top and a bottom, 4 2 of top to botion connections, and some local 
connections on these sides. 


Find: A direction -- left or right ~ for cach top to bottom connection so that the resulting total 
vertical dimension is minimized. Each connection is made by going around to the left or the 
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right as indicated by the direction. 

In this description, connections which originally went from the top or bottom to an adjacent side are 
considered to be local connections which go to the very edge of the top or bottom. (The portion used to 
turn the corner is not of interest here.) 

To solve the above problem, we first reducc it to a problem of assigning 0/1 values to clements 
of two vectors so that a matching on the vectors is maximized. The vectors represent the top and bottom 
terminals. ‘The representation removes unnecessary information about local connections. The value of 
"0" or "1" represents the direction “left. or right-- of the connection at a terminal. The matching 
identifies which segments will share the samc channel in the final routing, 


The major phases of the algorithm to solve the new assignment problem are as follows: 


1. Partition each vector into regions within which matching can be localized. After assigning 
values within a region, the regions are recalculated. The algorithm iteratively assigns valucs within these 


regions until there is only one unassigned region for cach of the top and bottom vectors. 


2. For the remaining top region and bottom region, the algorithm, assigns values from left to 
right along the top region, maintaining certain properties of the number of “Q's and "1"s in portions of 
the vectors. These properties picid a maximum matching for the vectors. If may happen that not all 
Properties can be satisfied simultaneously when same element is assigned a vale, At this ree say a 
failure has occurred. When a failure occurs, the algorithm may still be able to determine an assignment 
which suffices to guarantee a maximym matching. However, it may be necessary to try. both values for 


the assignment. 


3. When the algorithm must try both choices for an assignment, it docs not treat both choices 


equivalently. For the choice of value "1", the algorithm is ‘applied recursively to complete the 
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assignment. For the choice of value "0", the algorithm is applied with modification. We assume that the 
choice of "0" leads to a better solution, not just as good a solution, as the choice of "I". Therefore, as 
soon as there is evidence that the "0" choice will lead to no better a solution, the algorithm may stop 
pursuing this choice. In particular, if the situation in which the algorithn must try both choices rcoccurs, 
it will turn out that the second "0" choice leads to no better solution that the first "1" choice. This second 
"0" choice can be climinated. Only the "1" choice is used, and the search done by the algorithm is thus 


bounded. 


The technique of bounded search is crucial to the algorithm. Without the ability to bound the 
search, the algorithm would have ‘exponential ‘running time. The following sections described in more 
detail the algorithm outlined above. During this description, a large amount of notation will be 


introduced. The appendix to this chapter summarizes that notation. _ 
5.2 Reduction to Maximizing Matchings 


The Top-Bottom Routing Problem is reduced to the problem of assigning valucs within two 
vectors to ‘maximize a matching. The vectors, T(1,m) and BU1,n), represent the terminals at the top and 
Bottom, ftespectively, numbered from eft to right. (We will drop the T and B when it is clear from 
context whether we are referring to a top or bottom terminal} A value of "0" or "1" is used to represent 
the direction of the connection at each terminal.” Define the value fisnction VT: T(i,m) —> {0,1,2} as 
follows: | “3 | | 

VIG) = 0 ifthe direction of the path from terminal i to its pair isto the left 

1 —_sif it is to the right 
? if it is undetermined 
Value function VB: B(1,n) > {0,1,?} is defined the same way. 


Let p:TU.,m)UB(1,n) - TU,m)UB(,n)U{*} be the pairing function. ‘Tcrminal p(x) should 
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be connected to terminal x.. When p(x) = * (don’t care), the terminal x is connected to.a terminal on an 
adjacent side. Initially, the pairing function is known, and the. value functions have 0.or J values only for” 
the local connections, We would like to find: VT and, VB for. which all.directions are determined and are 
consistent with the initial values given. We will say, that.a value function V' is.consistent with V if 
V‘(i)= V(i) whenever Vn#2. We also say that V' extends V. For any pair of value functions we require 
that V1(i) = V.BQj) if p(T) = BG). 

We would like to balance paths to the left with those to the right on the top and bottom 
simultancously. [et my, be a function matching terminals in Tm) of value"0" to terminals with 
higher index in T(1,m) and of value "1", ic. matching paths to the left and right. ‘The function is formally 
defined as a snelaame partial function from THgn) to. ig such that if m,,(i)=j, then. VT(i)=0, 
VIG)=1, and j >i. peline My, similarly. For a lanl eal VT and VB, . My y be the maximum over all 
matching functions my, of Irangetm,, ph ic, My, is the nae numniber of matches among T(1,m) for 
a given VT. Parameter Mvp is defined similarly. 

We will reduce our oudhe problem toa problem of assigning "O's and “I's to maximize 


My7+Myp. This is justified by the following lemma: és ai . . 


Lemma 5.1: Given a Top-Bottom Routing Problem with local connections represented by value 
functions VT, and VB,, the problem of finding. VT, and VB,under which all directions are defined and 
for which the vertical dimension of the Creunnenbiae lr cad is minnie’ is equivalent to the 

probicm of finding VT, and VB, suc that My, + Mya, is prayed over all VIT._and VB which map T 


and B to {0,1} and are consistent with VT, and VBy 


Proof: ‘The matching function my, corresponds to determining which horizontal segments above the top 
~ side will share the same channcl. If my, (i= j then in some channel, the segment ending (goi ng left to 


right) at terminal i is followed directly by the segment which starts at terminal j. Each channel begins 
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with a segment running from the top left edge (a segment for whicti Mr {i)=0 dnd’ ptiERE,n) or pi=*) 


or with a segment which beging at a I-valued tert . 
segment in each channel is a'scgmient funiving to the Bp rit bags CVT Ayet and p(h€Bt,h) or pi)=") 
ied ander’ yr. ‘Figure $.2 Mlustrates 


or a segment ending-at a Ovahtied terminal which’ is’ uinma 
sbbive thE Wide: thercfore, we have fou ea 


possible configurations. ‘The channéts tah be th iaiy-ort 
assignment of segments to channels for cach Myr The argument cari ‘be ‘hevciied ‘so that we find a 


matching for each’ assignment to channith, “ Thett“is'd oite-t8*ne coltespondenc’ ‘betwech channel 


assignments snd -enatchings giving: 


number of channels used = ‘number of segments going to ie ae + umber of dameictied “Is 


This gives: total vertical diviension i edie to matching functions i ae, | 
© ht of reetangie'(denoted Hy + # chvasini toy 4 came tito a 
h + # top scgments going to left edge + # unmatphed “1 'sapiop:. 
+ # bottom segments going to right edge + # vomaiched “0"s at botiom 
ht # top segments going to left edge + #-V sawp: Iangettiyy | | 
+# bottom segments going to right edge + # “Os at He : (im, 


Note that: # "1"s.at top. = _# top samen going Hah ode + dap segments going to neither edge 


and therefore: “A top segments pois weft + "1s at tip 


2 op teach 


“gett F 


This gives: ‘minimum al vera dimension oer ll matching functions fr VF and VB, 
ah + # top segment i W beeen selirocld pee vi Mya) whe #C- m™,, + Mv) 


Where C = (# local connections top and bottom) + as ee a ica 


instance of the routing probes. Therefore, minienizing the, heighy, is equivwentto-masimizing 
My; + Myp- TREN yes 2 ae ae “sis aaa Pa Eth a gee hy ee B.. 


| We now present the various phases of the algorithm which finds value functions VI, and VB, 
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Figure $.2: Possible configurations when filling channels. 


pli) = + pj) = Bk) 
Penner channel 
VT(i) = 0 VI@ = 1 
myfi) = j 
pti) = j oo pea o 
_ Channel 
VIqi) = 1 v1@) = Q VI(k) = 1 VI(q) = 0 | 
unmatched — ~ my) = k ' unmatched 
pti) = Btu) | pq) = k 
— channel 
VI(i) = 0 VT(@j) = 1 ‘VI(k) = 0 VI(q) =1 


my() = j My (k) = q 


maximizing (Myr + Mya). 


5.3 Defining Regions 


The following description will be in terms of T. An analogous development is assumed for B. 
Let T(x,y) denote the terminals from x to y inclusive. In any tet interval TC.i, for any i<m, we would 
like to maintain the property that no more than half the seoninals are I-valued. This property must hold 
if every !-valued terminal in (1,i) is to be matched toa O-valued terminal with smalier index. If more that 
half the terminals in a left interval are 1-valued; sas of these 1-valued terminals cannot be matched 
| under any m,,,. If at least half the terminals are 1-valucd, any ?-valued terminals in the interval should 
become 0-valued if we are to maximize the number of matches. Sienifarly, we want no more than half 
0-valued terminals in any right interval to insure that these O-valued terminals can be matched. Because 
of the local connections, it is not always poaibie to fiabutady, these propertics en any left or right interval. 
Instead, we define regions within which these bounds hold. Within these regions matching can be 
localized. | 

For a given vatue function, VT, ket ZEROS,,(S) = {i€S[VI(i)}=0}; ONES, (S) = 
{i€S|VT(@)=1}; UNDET,{S) = {i€S|VT(i)=?} where SCT. Define the property iy 
(similarly OK-0) wien) is true if and pel if JONES y (0 s Morr + Da, (We. us the notation 
"ONES,,,(x,y)" rather than "ONES, -(Pxy))") Also define Put (VT, 1y) (similarly Full-0) which is is 
truc if and only if IONES,, (x.y) 2 Ta(y-x+1)1. Property Full- "BVT | is truce when at least half 
the terminals in ‘T(x,y) have value "b" under VI. Note that Full-1(VT,x,y) and OK- “U(VT.x,y) can be 
simultancously true only if IONFS,,..y) = “A(y-x+)), an integer. Full-1(VT,1.i) indicates that the 
terminals in UN DET.) should become 0-valued if we are to maximize Pe number of matches, 

We will now define the regions of T within which matching can be localized. iaseue under 


VT are formed scanning T from 1 tom. A new region begins when the previous region has at least half 
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"]'"s. Similarly, right-regions under VT are intervals of T scanning from m to ] and counting “0"s. © 
First define functions jan ry, mapping T to itself: 
hy) = 1 
hi) = i if Full l(VT by p(-Dai-1) 


ly,Ai-l) otherwise, fori > 1 


ry-{m) =m 
ty ,G) = j if Full-O(VT.I+ 11+ )) 
refi} otherwise, for Can 
The function 4,, induces an equivalence:relation on-T(1,m) under which two terminals i and j 
are equivalent if and only if hy) = lyf). The resulting equivalence etaases are intervals of T(1,m); the 
lowest element of each is an clement of the range of 4... These equivalenite classes arc called /eft-regions 
under VT. A new region begins when the previous region has at least half "1"s. Simitarly, r/,. induces 


right-regions under VT which are intervats' whose highest element is anelement of the range of r,/,- 


Lemma 5.2: (a)OK-1(VT,/,-(i)i) unless. /y,(i)=i and VIQi)=1 


(bOK-O(VT ir, (i) unless ry =iand VIQ)=0. 


Proof: We will prove (a). The proof for (b) is analogous. 

If 1) =+ and VI(@*l, thea [ONES (2 0=L 4(i-i+ 1) and OK-1(VT,ii). EF A Ki, 
then not (Full-1(VT,/,,(i-1),i-D))._ In this case \)=h i) and 
[ONES,, (4,41. $ ONES, Ay_(i-D.DI +14 Fated A +1 = L&e-h, fi) +)I+1 


and OK-1(VT,i,,{i),i). | oO 


On all the left-regions except possibly the last, at least half the terminals are 1-valued. Such 


regions are called full left-regions. We define a delimiter for the full regions. Let L,,;=m if 
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Full-I(VT,/,(m),m) and 4,,(m)-botherwise. Let Ry, Tif Ful-0(V1; Ln) (3) and 6,{1)+ E otherwise. 


Lemma 5.3 shows that the full left-regions can overlap the full: right-regions‘only when all terminals in the 


interval of overlap are 0/1-valued. 
Lemma 5.3: [f Ly, 2 Ryp then UNeee vitRypl a= 0, ie. al ietions are known on (Ry Ly). 


Proof: Assume L,, > Ry, The proof counts the number of “0"s and “1"s necessarily in (Ry, yp) 
using the definitions of Ry, and Ly 


Case 1: JONES, (Ry ply 2 Enos vy Suppose Ly is in rangetty,) Then 


(Ry ,-Ly,7) is a set of right-regions of T. Since Full is true for any  fight-r0gion of terminals greater than | 


Ry>p it must be that FulH (VT,R,,.L,7). However, we have assumed that there:are at feast as many "1"s 
as"0"sin (Ry pLy;). This gives | | 

JONES, <RyL HW = = [ZEROS yp Rypphy pt = = Mb Ry t I) 
and no elements of Rypl yy have value "7" under VT. . 

We now suppose .L,,- is-not in the range of r,,. Thea 

EROS Alyy + Lt lyst DW CT Mey yt lye + D+D1 and 
typly pt Dar fLy,) Interval Rypiyfl yp is a set of rahcreaions and Full0 holds. This, gives: 

ZEROS) (Ryply pM = ZEROS, (Ryfy fly H- PROS (Ly, Liyg(lypN 
IZEROSYRyply pl > FC y yp Ryp DIEM yy Ly At D -LUtLy Ry HDs, Le 
Full-O(VT.Ry 1-7) again implying no-eloments af (Ry. ij) hae vatoo "Funder VT. 


Case 2: JONESY (Ry plyp < PROS (Ry Lye 1 Thi case is s proven in the: same ‘manner as case 1 
: ; Neadbhy eee | Aha 

with the roles of ie and left-regions interchanged with nt sles of Ry, and raheem, ‘a 

The matching within the full le-regions (I lyr) or the fat right regions (Ry ,m) can be 


maximized meee: of the rest 1 of T. Our grt uses this. propety to break up the enor 
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Theorem 5.1 formalty states the independence of the regions. 


nearen 5.1: Let VT, bea given initial value funtion 5 dening Ry, and lyr chet every VT consistent 
with VT, and every matching function » Myp meres isa paerine function m’ vr agrecing with Myp On 
Q-valued terminals in (yy, +1,m). but. saghine sats 0-valued ieninal in a, ivr) to a 1-valued 
terminal in (lily }and such that jrangefm’ > frangi( ‘Analogously, there is a matching 
function m’,,, agrecing with m,, on I-valucd terminate im (LR yz -1) but matching each I-valued 
_ terminal i (Ry. #) to 2 O-valued terminal in ie ead range ls at feast as targe. 


te 


Proof: We will only prove ei existence of m'yr- The proot of the existence of m"\, is similar. The 


; proof’ is by induction on the number of ote-eion in a Leyty ) 


Basis: (yy, ) contains at most one region. tee 7 
If Pa then (yy) is empty. If tyr 1, then VT,(D=1 a (Lys) does not contain 
any "0"s. In either case, the theorem is trivially true. . 

If Lyp, > 1, then ber flvn 1. Thetefore, OK-HVTgLiy;,). We also have 
Fult-l(VTy1.Lyp,) and can conclude JONEB 5 (Ly = ‘Ly ‘For any value function, a 
maximum matchiag:can be found by scanning T from 1-to m. Each time a "I" is encountered, it is 
matched to the lowest numbered yet unmatched "6". As tong as there arc no more "1"s than "0"s in an 
interval (1,i), there will be an yet unmatched “0” in (1,1-1) to-match a 1-valued terminal i. In the case we 
are considcring, OK-1(VT,,1,i) holds for every iSly;.. Therefore every 1-valued terminat can be 
matchcd to a 0-valued terminal if all ?-valued terminals in Cleyp, )-beeomne O-valued. Since there are 
’*Lyy ds -valued terminals in. (1, Lvt } under VT), such a matching function would match all 0-valued 
terminals in (1, Lv, ). Fer cach 0-valued or ?- valued terminal under VT, in (LE yy 1) we can associate the 
1-valucd terminal to which it would match under the above matching. Given any value function VT 


consistent with VT, and any matehing function, m,,, the desired tatching function m’,, can be 
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created by matching cach 0-valued terminal inl yy }Me the. the. T-valued terminal it js associated with 
and using the meetin defined by ™vr for 0- valued terminals i in (vty am). since m’ vr marci at jeast, 


as many "0"s as Myp lrange(m’, vot is at least as anes as angen (Real that matching functions 


sone trone so ha te domain of machin fusion he same se as age) 


Induction: Cicada cain By the argument used in the: hae, the matching of the "O"sin 
this region can be restricted to the *Ys in. sion. “Fherefore, we may someve the first. region and 
consider the remaining tcrminals.as a new prabem. : ‘The peneen the fermaining. terrhinals. are 
unchanged and we ee the aeye assumption The desired m’ Myr! uses the palate of baie in the 
first Ieft-region to “1's in the first left-region and the matching of “O's in sth remaining f full et regan 


to "Is in the remaining full left-regions stained t by the inductive ssmumption, oO 


5.4 Assignment within Regions 


Given Theorem 5.1, we can show that it suffices to consider functions VT consistent with ¥T, 


which assign all ?-valued: terminals ia, Lvs, ) the value and all t-vabied terminals i in (Rvs, zm) the 


value "1". This also holds for VB, Theorem 52 states thls formally. ‘The algorithm for aesigning vaiwes | 
begins by assigning these values to the originally ?-valyed terminals. and their pairs.“ When a conitict 
arises, i.e. THEM. L yy) and Bj) = PATADE(R yy) or TAEER yam abd. BG) =. MROKECL yy 
cither choice can be made. We choose to.define the algerithes so. thet aesignments ace mada in the 


1. Make all?-valued terminala in(LLy, )and their botiom pair O-valucd. 
2. Make all valued terminals in (Ry, anand their bottom pales I-valued. 
3. Make all ?-valued terminals in (I,Ly }and their top pairs O-valued, 
4.Make all ?-valued terminals in (Ry, 0) and aheir top pairs I-valuod. 


This order gives a preference to assignments dictated: by the top segions. Executing each of 1 
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through 4 defines new functions VT. and. VB which extend VT, and VB, and which may be extended to 
value functions achieving the maximum, matching (Miyg + Mya) These new. functions may induce 
new Ieft-regions and right-regions, The algorithm repeatedly computes new L, Ry Lyp, and Ryp, 


and applies onc of | through 4, cach in turn, until there are no full regions containing ?-valued terminals. 


Theorem §.2: Let VT and VB be any ate of value functions consistent with VI, ; VBp and the sine 
function, p. If there are ?- valued terminals under VT, in (il, L, V1 ) then there are functions VT" and VB’ 
consistent with VT, VB, and p, such that all ?-valued terminals i in (1, oe ) under VT, are 0-valued 
under VT’, and My." + Myy 2 My, -+ Mype Similarly, there : are valuc dinctions assigning ?-valued 
terminals in (yr, m), (1, Ly) or vay sf) the values oo “Q", oF ey" ‘ewaclivel, without decreasing 


the sum of the matchings on T and B. 


Proof: We will only prove the statement for fixed ‘anita Ly): Proofs for the other intervals are 
analogous. Suppose there are ?-valued terminals in (Lily }under VT If these terminals arc 0-valued 
under YT, we arc done. Suppose there are such jerminals which.are not @-valued. under VT. Define 
value function VT" to agree with VT on (Ly + 1), but assign each ?-valued terminal-under VT, in 
(LLyp) the value "0". Let my, be a maximum matching function for VT which-matches each O-valued 
terminal in (L.Lyp,) under VT to a I-valued terminal. in (ILyy,)» By Theorem 5.1, such a.m, must 
exist. Let my, be a matching function for VT’ which agrees with Myp on (Ly_ + Ln) and matches 
each O-valued terminal in (IjLyp,) under VI" to a I-valued terminal in. Gibyp)- Again, myy is 
guaranteed to exist by Theorem 5.1, Then, 
range(m | - rango(m,. = [ZEROS yA Lyy. y- ZEROS, (LLy; Yad 

But, My, 2 lrange(m,,) and: My = [range(m,,}. Teorales, Magy: vr 2.4. The corresponding 
VB’ differs from VB by changing at most d"?"s or "1"s to "O's This can destroy at most d range values 


of any Myp giving Myp- My,‘ sd. Therefore Myr’ + Myp' 2 Myr + Myy . 
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“It remains to assign 0/1 vaines to any talied terminals of © TlLyy +1.Ryy “1) and 
BiLyp + LRyp -D. where VP, and VB, are the last extensions. of VT, and VB, Obtained by the above 
‘assignments, Intervals T(Lyp +LRyp “land BiLyp +1.Ryg -1) may be empty. Given Theorem 5.1, 
we can maximize the sum of riatches within T(Ly, +1Ry_ “l) and Bil yp + LRyp “1 independent of 
the full regions. We remove the full regions and define a new T and B containing only the remaining 
alia We get a problem on vectors T(1,m’) and mK 1,n') with new pairing function, p’, and new initial 
value functions, VT", and VB‘y, such that 7 | | | 

| m' = Rye lyr -Land n' = Rint, “1, 
Vi‘ g(i) = VEL yy, +i) and VB g() = vB yeti 

The pairing function, p’, corresponds to the the old pairing function when both terminals of a pair are 
within (ly + LRyy and (yp + LRyp 1). If a terminal within one of these intervals was originally 
paired with a terminal outside the intervals, then the new function pairs this terminal with * (don’t care). 
Any terminal originally paired with a termninal outside these intervats is 0-valued or I-valued under VT, 
or VB,. Since the pairing function is only needed for ?-valued terminals, changing the pairing of such a 
terminal to * does not change the problem. The new value functions, V1", and VB',, induce exactly one 
left-region and one right-region on each of T(1,m’) and B(t,n’). ‘None of these regions is full. 

The assignment to the new vectors is made using procedure SCAN-ASSIGN. <The procedure is 
called with inputs (TG,m’),B(L:n')p'.VT",,VB’,). "Thts’ procedure oe top: ie ‘reassigning — 
top-bottom pairs of ?-valued terminals the vatue "O" whenever this assigtiment would not create a 
bottom right interval, B(s.n), with more that half O-valued terminiits. When such an intérval would result, 
SCAN-ASSIGN ties to reassign the terminals the’ value "I". IP this would create a top left interval, 
TU1,q), with more than half T-valued terminals, the [procedure stops and retarns "FAH". The procedure 
which called SCAN-ASSIGN must handle the failure. If no faitute oceurs, SCAN-ASSIGN continiies 


assigning until there is a full top right-region, Le. a right interval with at least half “Os. Any zemaining 
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?-valued terminals are assigned the value "1". The procedure SCAN-ASSIGN is defined in Figure $.3. 


The following two lemmas. and theorem prove the correctness of 


SCAN-ASSIGN(1(1.m),B(1.n),p,VTyVB,) when it succeeds in assigning "0s or "1"s to all ?-valued 


terminals in T(1,m) and B(1.n). 


Lemma 5.4: Lct V1; and VB, be the value functions assumed by variables VT and VB after considering 
terminal ‘T(i) in either the first or second FOR ~ loop during the exccution of 
SCAN-ASSIGN(T(1 mm), 14 ILm)p.V" fy VB) Tf VT, and VB, do not define any full regions on T or B, 
then for any i, 1<igm, VT; and VB, natty the following properties: 

For all k, 1¢k<m, OK-1(VT,1) = 

For all k, 1$k-$n, OK-O(VBjk.n) : 
Furthermore, if SCAN-ASSIGN(T(1.m),B(1,m),p, VT, VB,) enters the second FOR loop, then for all 


value functions VI; and VB, where TG) is considered in this loop, PERS 1m)|J=F'%4m1. 


Proof: We first prove the shpeities for functions defined in the first FOR loop by induction on i. 

Basis: The properties are true for VTy and VB, by hypothesis. 

Induction: By inductive assumption, the propertics are true after considcring terminal T(i-1), ie. 
after the (i-1)" execution of the loop. Function VT, differs from VTi, at most for T(i). For this to affect 
OK-1(VT;,1,k) for some k, the following must hold: JONES), (LK) =L ak, k>i, VT,_,@)=?, and 
VI(i)=1. However, if JONES; (= Lk, then VT,(i) is not assigned "1". Therefore, for all k, 
OK-1(VT,,1,k) holds. For OK-~X(VB,.k,n) not to hold for some k, it must be that: ZEROS), (k.n)l= 
L&(n-k + 1)J, k Spi), VB, (p(i))=?, and VB(p(i))=0. However, if ZEROS), | (k.n)f= L'a(n-k +1), 
then VT;(i) and VB (p(i)) are not assigned "0". Therefore, OK-0(VB.,.k,n) must hold for all k. 

If SCAN-ASSIGN enters the second FOR loop, Iet h be the first value of j considered in this 


loop. Execution of the first FOR loop is completed because Full-O(VT,_,.¢,m) holds for some q. We 
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Figure 5.3: Definition of SCAN-ASSIGN 


SCAN-ASSIGN(T( 1, m), BL, n)p.VI .¥B,) 


/ Assigns 0/1 to ?-valued terminals in T and B where Ly, = =0, Ry, =m+, L om" Ryp, +1 and 
no region is full/ ics 


VP=VTy and VB: = =VB, 
FOR i = STEP 1 UNTIL. (3q)Full-O(VT.q,m); 
IF V1{i)=? THEN 
IF (3s) pli) Ris.n) and {ZEROS yal = =Li(n-s+ lJ THEN 
_ IF (39}i5q5m and JONES, yal ai= Liked’ THEN 
RETURN (FAIL i, vr Sag 
ELSE VIG): = Land VBQpQ): = | 
ELSE VI(i):=0 and VB(p{i)): =0 
END ag. Same 
FOR j = iSTEP 1 THROUGH m 
IF VTQ@)=? THEN VTQ):= 1 and VB@p):=1 
END 
RETURN (SUCCESS, VT, VB). 
END SCAN- -ASSIGN 


Cie: war hs sree Sieg ang: 
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know that h>2 and qh since VT, _, agrees-with VT, an ‘Fh,m) and for no k does Full-O(VT,,.k,m) hold. 
If q#l, we know OK-(VT,_.1.@D holds -by the: first part of this ‘proof. ‘This implies that 
ZEROS, (141K 2 F%4(q-1) TV since no terminals in T{I,h-1) ave ?-valued under VT,_,. Then 

ZEROS), (1 my > F'4(q-1 Y1 + FMem-g +1 2 Pe . 
However, PERO. if my < rim otherwise, the fint laop would not have been repeated for h-1. 
Therefore, 7 ied Ae m)| = [%m1. the second FOR i does not assign any "0"s. Theres: 
for all VT;, h<j<m, VFEROSy, m)| = Thm. Alo, since OK- VB ,- vk m) holds for any k, — 
OK-Q(VB sk, m) holds for any k. 

It remains to show that Of (VT, 1 sk) holds sd al k. ‘Sine any VT; age with VI, 
T(,b-1), OK-1(VT},1,k) must hold for k<h. Suppose that for some ih. and some ejzh, OK- “IVT, 21k) 
does not hold. Since the number of I-valued terminal increases on cach iteration, OX -1(VT,,,1,k) must 
also not hold. Since no terminal is ?- valued under VT. EROS (, m)| = ['%m1 implies 
ONE iE (1, ee Lusmd; therefore os 


For any i>h, [ZEROS yy (iim) = HEROR (iim)| <7 sie aa 


This implies JONES, vi, (e+ Lm> Lm(mtk+ D+ Dd. 
Then JONES, (1b) < Limd- - L¥e(m-k)4d, 


giving ONES (1,k)| < r4k, contradicting our r assumption that OK-1(VT, rl *) does not hold. QO 
Lemma 5.5: For any call to SCAN-ASSIGN as ia Lemma 5.4, the first FOR loop is completed with i<m. 


Proof: If the lemma does not hold, then after i=m nin ine first FOR loop: 
ZEROS, “(my <P. 


However, by Lemma 5.4, JONES, (1m) § Lam, implying ZEROS, (1,m) > T'4m. QO 
m m 
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Theorem 5.3: Let VT,.and VB, be initial value functions which induce no full regions oo T and B. If 
SCAN-ASSIGN(T(I.m), B(La) p, Vig, VB,) returns (SUCCESS,VT_VB_}: thes VI_ and VB, 


maximize My. + My, over all VT and VB consistent with VT, and VB, ne 


Proof: Since OK- WOT, Ly) holds for Isksm and ONES a sn = Um Myr = L%m, the 
maximum possible Since OK-0VB,, ik as holds for Istsn cach Oralued bottom terminal can bd 
matched, and Mvp = VEROS yy al ay. If Mvp = Land then no lrges matching is , possible and 
we are done. em Myp <LYnt. Let Myp =LY¥4at-A, An. lea ues and vB be =i other ae 
fenctions consistent with p, VI, and VB). mee must show that 
| Myr + My SMyy_ + Myg = = LMA +Lumd. 
Since the initially ?-valucd terminals are assigned a values i in top-bottom pir we have: 
IZEROS,, (mn IZFROS vagClnd = FEROS my PEEROS yy ay 
= rama +d and 
IZEROS, ,(1.m) - [ZEROSy, (Lm) = FZEROS (0 (ZEROS nah 
Let [ZEROS,,(1.m} = rem+z, where 2 any saat Thea: 
ZEROS, _(1.n) = (ZEROS, (104 (EROS, EROS, (L0N) 
= f Ham 1+ 240 em Yo A)= Lead+rA. 
Case 1: 2<0. Then My; + Myp < ZEROS, {1.m)| + ZEROS 1.0) 
<< Pam t+2+b beds A Cie Leb A, 
Case 2: 0$2z. Then My, + M va SNES (Lm + VEROC slot i 
< LimbrtLinlted = = Limit Lind. | ‘Oo 


It remazns: for us to deal with a retara: of (FAIL pele Before describing how a failure is 
handicd, we prove that for certain value functions, the value of a terminal can be switched from one of 


"0" or "I" to the other without decreasing My,+My,. ‘These value functions are consistent with the 
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initial value functions used as input to SCAN- ASSICIN and assign only "O's and "1"s. ave assign an 


imbalance of "0"s and "1's in some left or right interval. We will use this ability to change the values of 


terminals to prove that a value functions returned by SCAN-ASSIGN do not assign too many or too 


few "0"s and "1"s to sahleve: an optimal matching. 


Lemma 5.6: Let VT) and V1}, be initial value functions which define no full regions on T and B. Let VT 
and VB be any function conor! with Vly Me and P wader which no pies is SVR: 

A) For any x, Istsm, if JONES, he [ZEROS vill x» > 1, then there is a terminal i in 
T(x) such that VT,()=7, VIC) = 1, and the functions vr and VB’ obtained by changing the values of i 
and p(i) to "0" are such that Myy'+ Myy'2My.+Myp. | 7 

B) For any y, 1SySn, if IZEROS, (y.n){ - JONES, ,(y.n)I > 1, then there is a terminal i in 
By.) si that VB,(i)=?, ViKi)=0, and the functions VT’ and VB’ obtained by changing the values of i 


and p(i) to "1" are such that My-'-+ Mya) My7+Myp. 


Proof: We wilt only prove A. The satiot for B ts analogous. Since een ees 
there are at least two ie in T(1,x) unmatched under any matching function for VT. 

Let there be a matching function achiesing Myr for which the terminal of lowor index among 
two unmatched "1"s is ?-valued under VT,. We ‘will dow by contradiction that such a matching function 
must exist. Given such a matching function, change the valuc of the terminal of lower index to "0", 
changing the value of its bottom pair as well. This defines VT’ and VB‘. The new "0" can match the 
second unmatched "1", giving My~=My,+1. In B(1,n), at most one range value of any matching 
function has been destroyed; therefore, Myy'2Myp-l.. It follows that My7/+ Mya 2 My, +Myp. 

We now show that the matching function described above must cxist. Figure 5.4 illustrates. 
Suppose that for any matching function achieving My, at most one 1-valued terminal in T(x) which 


was initially ?-valued (i.e. under VT) is unmatched, and that if there is such a terminal, it is the terminal 


Figure 5.4: The proof of Lenuna 5.6. 
Claim exists: a 
ae en es See cea en 
“yj ; “y" Pane tne iz Mee edt. 8 
unmatched unmatched 
7-valued under VT9 
Otherwise: 


any 1-valued terminal under VTp is matched to a 0-valued terminal here 


- 1 ; ‘i ee ae 


all 1-valued terminals under VT are also 1-valued under VIg . 
a many "1"s under VT 


all 0-valued terminals are matched to l- valued terminals here 
Figure 5.5: Definition of C. 
‘set C includes all "1"s ” 


assigned by SCAN-ASSIGN 


of highest index among all unmatched "IMs in T(x). Then, for any matching function achieving My,, 
there is at least one 1-valued irene under VT, which is ened: Consider all maching functions 
achieving My, for which ; a maximum euaiker of valued sortie in Ta) are © matched. Among these, 
consider =“ those aivctiid for which a maximum number of originally ‘valued, now re valued 
terminals in Ta, x) are sineebed (either zero or one mh eins. For each of ie note athe indices 
of any terminals in T(1 x) which 2 are + uerneched and are 1- Srahued under Vt, Cae the matching 
function for which the spice index, say U, a obtained. No nly -valued terminal in tl, u-1) which 
is 1-valucd under vT i; unmatched. bisa O-valued terminal in Td, u- in must F amatch an yall 1-valued 
terminal in TL u- », Otherwise, the bachine function could be modified so > that the fending O-valued 
terminal matched u. If the offending valued terminal had bien unmatched, this would pinidiice a leer 
matching; if it had matched a terminal in T+), ‘this wonild produc a matching under which mere 
1- valued terminals in T(1,x) are matched. If the offending temminal had matched an nally i valtied 
terminal in T(1,x), matching it to u would produce a matching with | a bitect vitae of unmatched 
* L-valued terminate which are ?-vahued under VIy; if it had etched a terminal ir in Tu+l x) which is | 
1-valued under VTy, this would produce a matching with an caniached termina which i is 1-valued under | 
VT, and of higher index thain U. Any of these possiblies coateadict 6 our chinice of matching fusctioh, 
There can be no o matched 1-valued terminals in TO, u-t) which a are initial 2-valued nor any y unmatched 
terminals of this type. All 1-valued terminals in TU eae are = hid 1- ake! Since all O-valued 
terminals in T(1,u-1) must match 1-valued terminals in ee we have: 
[ZEROS yr wl S JONES y7(1 u- 1) = JONES, (l, wl 
Since no terminals are ?-valued under VT, this implies ONES, Q, u- I ar ‘wD, aes the 


hypothesis that there are no full regions sndec VT>. | 0 


ps ad gee: = * 2 = one oS gaat oat a ve Bk ai 
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5.5 Handling a “failure” 


Let us suppose tha SCAN- “ASSIGN one (PAIL . vr, 7.¥B} ). We w wil now describe what is 
ksows about T, B, and their valve functions afer SCAN- ASSIGN returns We wi wil use Se this information 
to handle the failure Procedure SCAN-ASSIGN iscniind s the vale, k of bop variahie i | when 
SCAN- ASSIGN encountered ied lure, Therefore we : know that Tips yee. vB, 1, are the value 
functions returned. There i isa ak such that JONES, a ay = Lad We wil lise the smallest such, 
q- | There is an s satisfying te lowing three properties: o mins “ EROS VB. say = 
Lens Da and (ii) cach terminal in m0, k: 1) which has Raed eet the oo a by 
SCAN- ASSIGN is paired with, a terminal in n Mm), ie. ONES, sal EDONES a EN OK6AD. 
This last Property: follows fom the fact that scan-assion only assigns ee! wen, a te a 
(AsxptV€ Bsn) iad ZEROS aco = on Je “Tere we can —-* an 8, with each 
?-valued 70) made r vaed: Once EROS (0 oi = = Lh Ds forsome VB, the ‘paar of ec oft 
in the interval dors not change in inter t nections: Therefore, chootng the galls ofan eae 8, and the . 


a7 


4 sociated withthe faire, we have an ind, which satis all hoe properties, In fac, we we the 


ae Tate. a ie! 
— s which stir all three properties, ‘Not iat VT, = =V8, 000)=2. \ We also now hat cach rs 
Say in Bian) which is signed a = or a ae value by SCAN-ASSION rust be pied with a 


Se +e 


terminal in TA, k-1), since no Iniaty aed termina in a Thm) or thei amie ne have been 


Lb te oy ad emi hn pi a Me) (Fae. 
For an Spine semenent ey sl renaining ated ermine nT.) ould bre“ bas For, 
an optimal signment a the boom, al emsning Paved minal in, An) shuld boone "I" 
Obviously, on C these are conflicting goals. We wish to find 0/1. asignnens for the remaining ?-valued. 
terminals in C (including k) for which some extension will achieve the maximum sum of matches. To do 


this, we first prove that we need only consider valuc functions which assign a number of "0"sor "1"sinC 


Sie estes te Ae SR oh Rio ge RTS atk epee Ns Re GRRE Ulta Sree Meee oe, Saag e i ee 
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within a certain range. Let: 
oF JONES yy A, 
= EEROB ys (CH. 
y= [UND RT, (OF 


To simplify the statement of thie ierresaa to follow. | ‘et * “assume the standard state after 
SCAN-ASSIGN returns ‘FAIL mean: Assuming VT, and VB, are initial ‘afte fuinctioiie which define 
no full regions on ‘F and B, ket SCAN-ASSIGN( Tm), BO,n), Pe La ove return (FAIL k, VT,. 


vee Pao rite ty bom we hav Ged em tone” 


Lemma 5.7: Asqume the standard state after SCAN-ASSIGN. rowrns: "FAIL". Let. VT and VB be 
arbitrary value functions consistent with VT y VBy and p uodgr-which there are no 7-valucd terminals. If 
c,>1, then for any x, 1$x $c,-1, there are valuc functions VI" and VI’ also consistent with VT, Vi, aad 
p under which there are no ?-valued terminals such that; — 

| Myr + Mya 2 My; + Myg and JONESY (OL = c+ x, 
Ifc,= 1, there are such functions VT" and VB’ for which c, < JONES,,4O) = ¢, +1. 


Proof: Ifs > 1, tet D be the set of initially ?-valued terminals in Tq) whose pairs are in B(1,8-1). 
Then: JONES,,,(1.q¥ = IONES,,{C)| + JONES, (D) + JONES, (1.0) | | 
Note that JONES; (1.9) = JONES, (Lay + JONES, (C) a Lua 

Case 1 JONES,,(C) > c, + x, where ifc,=1, x= 1. | 
We use induction on JONES, (C+ JONES, (DJ. 
Basis: JONES, (C)| + [ONES,,(DW = c, + x. Then JONES, (C)| = c, + x as desired. 
Induction: Assume that the lemma holds ift . os 

c, +x < JONES, (Cj + JONES mNSetaeL A i>0. 


When lONES,,,{C)| + JONES, (D)f = c,+x+i, then: 
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JONES, (LQ) = ¢, +x +1 + JONESy, (Lqy= Liq +x+iD Lig +2 
IFJONES,,(C)| = c, + x, we are done. Otherwise, we use Lemma 5.6-A. 

ONES, ,,(1,q - ZEROS, a> Lqd + 2-(F Mql-2)>3 
Therefore, _ Lemma 5.6-A, there are VI" and YB” such that My + Mya > Myy + Myp and 
| JONES, (CY + JONES, (D) = JONES),,(C) + JONES) (D) - 1. 
By inductive aapeae there are VT" and VB" such that Myr" + Myg” 2 Myr’ + Mygy and 
JONES y7-4CH = Esk | | 

Case 2 JONES, {C) < Sc, +4, where goa x=0. 

Then IZEROS AC) ZC tc, + & (c, + x)) = = Cy + Cy-X. iam, let E be the set of initially 
?-valued terminals on Bés,n) whose pairs‘are in T(q+ Lm). Let p(C) denote the bottom pairsofC. 

IZEROS (sn) = ZEROS, (PC) + ZEROSy_(B + IZEROSyq (6.0) 
Weknow [ZEROS\p (tm = ZEROS, (H(CO1+ mOsm = Lia(n-s+1)4. 
We use induction on anes vat PCC) + ZEROS, (PA. 
Basis: |ZEROS,,(p(C))| + [ZEROS (EW = cy +¢5-x. Then |ZEROS\,(p(C))| = cy +c,-x, implying 
IONES, ,(p(C))| = c, +x as desired. : sa 
Induction: Assume the lemma holds for: 

Coteyx S [ZEROS oC + ZEROS fF Sptsya ti fori i (> 0. 
When [ZEROS yy(P(C)) + (ZEROS, g(E) = cytes ati: 
IZEROS vps) = Gyteyx+i+|ZEROS, 2600 = Lu(nst paegaeiZ Meat DLA 
If |ZEROS,,(p(C) = cy+c,-x, then eo = te ag desired. Otherwise, we use 
Lemma 5.6-B. _— os 

- (ZEROS, (sn) - JONES, (6.0) > Li(ns+ NI+240 W(ns+1)1-2) 23. 

By Iemma 5.6-B, there are VT" and VB" such that Mygy + Mya! 2 Mya + Myg and 


[ZEROS yy(H{C)| + ZEROS, (FN = IZEROSya((C)Y + IZEROS, (EM 1 
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By inductive assumption, there are VT" and VB" ‘such that My.” + My" > Myr + Mya and 


JONES, AC) = C) +X. . oOo 


We now know that there are valuc' functions achieving the maximum sum of matches which 
assign at least 28 ‘many "I's in C as SCAN-ASBIGN' has: siuigned. However, we do not know if the 
particular choices for I-vatued terminats wif fead to ani optimal'solution. We will now show that except 


for possibly one torminal.which-we can caslly find, SCAM-ASSIGN has made’ good choles. Note that 


the distinction made between c, +1 and'c,>F as in Lomna $:7 6 necomary: Wheat c>}, regardicss of the: 


extensions of VT, and VR, being considered, there are enough necessarily unmatched "1" in T(1,q) of 
unmatched "0"s in B(s,n) so that we can trade off top matches-againgt Bottom matches, When c,=1, we 
know that an optimal assignment wilthave at least one wArtatéNed:"}’* on T(1,q) or one unmatched "0" 
on B(an).: However, we Catiiot teathe these against eactbother because when a "0" is unmatched in B(s,n), 


the extra:"0" in T(1,.q) may match a1" in Feq-++1,m). Bimtlady;-an extra “1” in B(sin)-may match a “0” 


in B(1,s-1). The algorithm must uy both ienoies - 7 eves 9 one ueeuenee a in TQ a) or one 


unmatched "0" in Bn). 


The termidal which may have been incorrect asigned is ae smallest numbered a valued 


terminal in B(s,a) which is valued undet VB, Ken this is termina i "Given t two ermine: u and y, 


with u< v in T and p(u) < pv) in B, assigning 3 to u and "1" to v is always preferable to assigning “" 


to u and "O" to v. ‘This is because any terniinal whieh ver pt¥) can satets-when they are-0-vatued, WoOf 
p(u), respectively, can match when they are 0-valued; any terminal which can: be matched to u of pfu) 
when they are }-valued, can be matched to v-er pth iepactiiely, when they are l-vatded. Therefore, 
assigning "0" te u and p(u) and "1" to v:and-péy) gives at icast as large a-matching on each of T and B as 
the opposite assignment. ‘Therefore, if terminal tand-its pair: are-smalter: numbered than p(k) and , 


respectively, then it is better to: have terminal i O-vahved and torminal k I-valued. than vice versa. 


However, the proferred assignment is not consistent-with VB,, and VB, j, We define new value - 
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functions which maintain the important properties of VT, , and VB,.,, but:allew. the preferred 


assignment. 


Definition(sec Figure 5.6): Assume the standard state after SCAN-ASSIGN returns "FAIL". Let i be the 


terminal of smallest index in 1(s,n) which is ?-valued-under VB, and I-valued under VB, ,. Certainly, 


pti<k. If p(k), let h = i; otherwise, let h = p(k), By definition of h, all Jevalued. terminals ia C are. 


paired with terminals in 1h,a), Define VT, and VB, as follows. If h. p(k), then VT, = VT; and 
VB,, = VB,.,. Ith # p(k), then VI;, and VB, agiee with VT, , and VB, , except ath, k, ph) and p(id, 
where: 

Vig =1 igo =? 

VB, (ptk)) = 1 VB, (h): =? 


Note that OK-1(VT;, ,1.x) holds for all x, 1¢$x<¢q, and. OK-O(VB,, y.0) holds for.all y, sSygn. The 


“standard state after SCAN-ASSIGN returns FAIL” will now iacludeh, VT, and VB, as just defined. 


We will now prove that VT, and VB,, are satisfactory ¢ extensions of VT) and ee i.e. they will 
lead to a pair of value functions which achieve the maximum maiching My; Myp ate 5.8 ives 
properties of extensions of VT, and VB, which will be pease to prove that it is nucle to consider 


only these extensions. 


Lemma 5.8: Assume the standard state after SCAN-ASSIGN returns. "FAIL". -Let VT and VB be > 


extensions of VT and VB, 


A. If there is a ?-valued terminal in T(p(h),q) under VT,, which is }-valucd under VT, then, for - 


any l-valued terminal, i, in T(1,q) under VT, there is.@ mefching: function:m,,, achieving My which: 


matches cach 0-valued terminal in T(L,q) to a 1-valued terminal in T(1,@) aod leaves i unmatched. 


B. If there is a ?-valued terminal in Bfs,h) under VB, which is O-valued under VR, then, for any 


0-valued terminal, j, in B(s,n) under VB, there isa. matching function Mg achieving My, which matches 
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Figure 5.6: Definition of h. 


Case l: h =i 
T: 
1 pti) m 
: p(k) on 
B: 
all"1"s in here after SCAN-ASSIGN are "1"s 
under VBo 
Case 2: h = p(k) 
T: 
1 pCi) k q m 
| emer s p(k) n 
B: 


all "1's in here after SCAN-ASSIGN are "1"s 
under VBo 
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eid gE py lite oR a et chop SSN DC RN RR a RE RE eee ea 


each ]-valued terminal in B(s,n) to a 0-valued terminal in B(s,n) and leaves j uihmatched: © 


Proof of A: If suffices to show that [ZEROS (2.9) < (ONES) {x.q}{ for all'x in T(1,q). If this is true, . 
then we can produce a matching on T(1.q) which ee each “O" in TO a) toa "1" in T(1,q) and does 
not match a given 1-valued terminal i We define bt eaching function by scanning T(1,q) from q to 1 
and matching cach "0" emmgiantered ois yet unmatched "1" of higher indcx other than i. ‘Fhe fact that 
ZEROS, ,(x.0) < |ONES,,,(2,q){ guarantecs that we will find such an unmatched "I" for each "0". This 
matching can be extended to a staiching which achieves M, by letting the matching agree. with any 
matching function achieving Myy on "0"s in T(q+ 1,m). 
For any x in T(1,q), we will show that ZEROS, &x.q)< JONES ;(x,9)1-by showing that more 
than half the terminals in T(x,q) are 1-valued under VT. Suppose | < x < p(h). Then, by hypothesis: 
JONES (x.0 2 JONES, (xa +L 
Bee SE NONES yi C100 Hib Toda? dvs Keay Se ORV Egham toes 
HONPS,;, (ual Bbq] LMG-DJ > L%4(q-x+1d and : 
JONES, (x.q > Lwaxt Ds + 1, as desired. | 
Suppose pth) <x <q. We claim that JONES, (La-D<L&(s-1)4. If true, this gives: 
JONES y;(x.q 2 JONESY; (xa > LYegd- Lii(a-I)4 > Lisa + 1) as desired. 
If x > k, then JONESy,_ (1.x-1)] = JONESy, (a-1 CL4G-1)4. Otherwise, we would have chosen 
x-1 as. q, butx < a. Ifx < k, then since x > p(h), p(h)stk. In this case, 
ONES), (Ls-1) = JONES, xt “1 
since VT, (p(h))=? and VT 1-1(P(h)) = 1 and the value functions agree Srsepahets else on T(1,k-1). Since 


OK-1(VT,,, 1-1) holds, JONES yy (1.x-1 ¢ LM(-1). 


Proof of B: It suffices to show that JONES, ,(3y)1<|ZEROS,,(s,y)| for all y in B(s.n) so that all 


1-valued terminals in B(s,n) can be matched to 0-vatued terminals in B(s,a) without using j.- Produce the 


_ tt: 


pealehing by seonnbay tiney:s 0 8--AS te. the. praot oF Awe, peste: Miat secre: ean halt she terminelsin 
B(s,y) are 0-valued. For y > h, by hypothesis: _ are . 
IZEROS, (ay) 2 (ZEROS yy (sy) ew” 
We know EROS yy (and = L(g IJ and. ify Ca, OK-O(V Rpg 1m) holds, Therefore: 
ZEROS, sy 2 Lignat D- LMindy+ 1+ Dd Zh M(y-s+ 1d and 


AEROS (yl & L'*O8+ Dd +1, asd 


PROP EME ors pomarygt bs 


For y Ch, My + 1.n) contains p(k) and al] I-valued terminals in C under VB... Therefore,. if 
ZEROS yy, (7 + 1a) =.LMOry}4, we would have chagon y+ des 8 but ya 8 Therefore: 
2ROS ARF i44 ; te Ast S & AU we 
IZEROSyy (y+ L9M + EROS ya. (y+ L aS LY(Oryhd and 


|ZEROS, (sy) 2 ZEROS), (ay) >, b(nstDd - L&(n-y)d 2440+). O 


Now that we have Lemma 5.8, we can prove a two part theorem whieh, when combined with 
Lemma 5.7, proves that in ‘ouf séarch for’ an ‘optimal pair of value Ainctions, it is sufficient to consider 
only value functions consistent with extensions of VT,, and VBy. If c,>1, the esiensions give all 
-Aernnteale ca Tha) gee. Riga 0 A: values, if cy = 1,,the terminals h and ph) are the only terminals 


in these intervals whose values are not fixed. 


Theorem 5.4-A: Assume the standard state after SCAN-ASSIGN returns "FAIL". Let VT and VB be 
value functions condistent with VT,, VB, and’p for which no terminal is"?-valued and such that 
|ONES,,,(C)! = c,-+ey-1. ‘There are extonaioné, VT. and VB, ‘of VT, and VB, consistent with p and 
(i) no terminals are ?-valued 
(i) all ?-valued terminals in ‘T1.q) undet VT, with pairs in (1,s-1) are O-valued under VT, 


(ii) all ?-vatued terminals in B(@,n) under Vi, with pairs in (q+ I,m) are I-valued under VB, 
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(iv) if Cy2 I, all. 9- valued terminals in-C.under VT;,, except p(h> are 1-valicd-under vi. (cyl of 
them), and VT,,(p(h)=0 


()Myp_ + Myp 2 Myr + Myp, 
Proof: Properties (ii) through (iv) define all valued of VT, and VB,_ of ‘T(1.q) and B(s,n) except those 


of p(h) and h when c, = 1. Ife, = 1, we will begin by letting VT (P(h)) = VB,(h) = 0. ‘This may be 


changed. The valucs of terminals in 'T(q+ 1m) and B( $1) whose pais are in 14s,n) and T(1,q) are also 


Figure 5.7 illustrates. We stifl must specify the values under VT. and VB. of terminal pairs with one 
terminal of cach pair in T(q+ 1,m) and thé other in Kis-1). For 9-valued termtinals under VT, and VB, 
of this type, VT, and VB,, agree with VT and VB. 


We prove Myr + Myr = Myr + M,, by comparing the maximum matchings in various 


intervals and combining. Let My, be a matching function achieving My, and my, be a matching. 


function achicving My 


1. Consider B(1,S-1). Let S, be the set of ?-valued términals in B(1,s-1) under VB, with pairs in T(1,q) 
which are 1-valued under VB. ‘The only terminals in B(1,s-1) which are 1-valued under VB, and 1-valuod 
under VB,, are those which are ?-valued under VB, with pairs in T(q +1,m). This follows from the fact 
that any ?-valued terminal under VBy which is 1-valued under VBy is in B(s,n), and any ?-valued 
terminal under VB,, with a pair in T(1,q) is O-valued under VB... Functions VB,, and VB agree on all 
?-valued terminals in B(1,s-1) under VB, with pairs in T(q+1,m). Therefore, every 1-valued terminal in 
B(.,s-1) under VB... is I-valued under VB. Each of these terminals which is matched under Myp can be 
matched to the same 0-valued terminal under a matching function for, VB ey Tet Myp_ be a matching 
function for VB, which matches "1"s in B(1,s-1) as My, does. The matching under Hyp. Pt 13,0) will 


be defined later. We have: 


1 


Figure 5.7: Assignments under Theorem A4-A. — ae Ae 
“[')_ snamecnesempsreanert i TERA INSEAD TES AE REST RTE TAD) 


exce@th | . 
"2"s under VT, 
LE" $a MOE-W Eogy 


"2"s under VBg, 
"1"s under VBgy 


"2s under Vig 
"0"s under VV gy 
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jmatched "1"s in B(1,s-1) under Tp, | = 


matched “1"si in BQ), $1) uid ya |matched 1 sin 8, finde mg 


jmatched "1s in B(1,s-1) under mm, | > Imatched “In fiir) under my [S,I. 
cet. 4 : Pes 


ll. Consider T@q+1.m). Let Sy be be ~ set of Sioaia terminal under VIg in T(q+ 1m) with pairs in 
R(s,n) which are 0-valued ‘under VT. Each @- valued terminal ia Ta+l, m) under VT, is O-valued under 
VT, or is ?-valucd under VTg and VT, and is pairod with a terminal’ .") Bl. s-1). Therefore, every 
O-vatued terminal in T(q+1,m) under VT,, is 0-valued under VT. Every such "Q” which is matched 
under m,,, can be matched to the same “1” undér a imag Function for VT... “We let Myy_ be such a 
matching function on T(q+1,m). | 
| |matched "0"s in T(q+ 1,m) under my be 
jmatched "0"s in Tiq+ 1m) under my] - [matched "0"s in S, under m,,J 


_ [matched "0"s in T(q+1,m) underm,, | > [matched "0"s in T(q+ I,m) under my - [Sol- 
sf 


Hl. Consider B(s,n). Since VB,,(h) = 0, Lemma 5.8-B holds. Under VB_., each L-valued terminal in 
B(s,n) can be matched to a O-valued terminal in Bisa) while leaving any particular 0-valued terminal 
unmatched. Complete the matching Mp defined on B(1,s-1) in I by such a matching on B(s,n) which 
leaves h unmatched. We know that: - 
JONES, (C) = JONES, = ¢,+¢,1 and. 
[1-valued terminals in B(s,n) under VB,, which are ?-valued under VB, and have pairs in T(q+ 1m) 
= |terminals in B(sn) which are ?-valued under VB, and have pairs in Tiq+ Lm) 
= |1-valued terminals in B(s,n) under VB which arc ?-valued under VB, and have pairs in T(q-+ 1,m)} 
+ B,J 
Then, |matched "1"s in B(s,n) under mys. | = JONES yp (an) = JONES,,(s,ny + [S,] 
2 |matched "1"s in B(s,n) under my,| + [So]. 
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IV. Consider T(1,q), We sujet consider c,*= | and o>! separately. . 
For c,21, there are.other terminals: in: UNDEF, (C):beside hy These terminals must lie in 
-T(k+ 1,q), since no terminal except h in T1,4)-8 ?-valuad uaderVT,,. Since ke-p(h) and these terminals 
are 1-valued under VT a Lemma 5.8-A applies. Under VT... vory 140" :ie Teka) ‘¢an:be-matched toa 
“1” in TCL@). Complete the matching: defiaed: on Thqa-t nin Uf with such a matching on T(1.q). 
We know that: 
2PROS ACH = HERO, (OH ty #2 
}O-valued terminals in T(1,9) under. VI. which are ?-valued.under VT, and have paitsia B(1.s-1)} 
=.|?-valued terminals in T(1,q) wader VIE, with pairnin B(Le 1) 
_ = |@-valued terminals in-T(1,q) under VIE which. are Y-valued uader Vy and have pairs in BLS 
+ [8,1 
_ Then matched “O"sin T(1.q) under my, | = IZEROSy, Ahad + {ZEROS,,,(1.q) + [S,| 
> Imatched “O"s in Tha )ainder my]. +:16 
Combining I through IV gives: | | 
Myr, + Myp,, & lrangeliny, M+ lrangetnys Im leangelmy gh + bengetmy,) = Myy+ Myp 
For C= 1, since all ?-valued terminals in. T(t,q) under Vi, are @rvalued under VT, 
ONES AE) = ONES), (10h; Since OK: (VT, 1.x) holds for all x in. Pg), all L’%4qJ “1"s in 
T(1,q) can be matched. to 0's in: T(1.q) under VI... The definition of Myr on T(1,q) is such a 
matching. . 
jmatched "0"s on T(1,q) vader my | =L’qi = IZFROSy, (a) > (F 4q Tt. 4q4) 
2 (EROS (Lg + 1S) G HL MQ) 
> |matched "0"s.in TU,q) under my7}-+ 18,4 (F 4q Fk 4qJ). » 


(i) If q is even, the above implies, with | through Ill, that My, + My, 2 My, + My, as 
x x 


desired. 


hated RETIN PEGS 


Ni ee 


“In. 


(ii) If any of the inequalitics deduced in | theotigh FV. ts titan at 
Myr_+ Myp Myr + yy HO hase ing 
Mer +My 2 My Mig as Gebieed: © 
i Assume (i) ard (ii) do‘dot hold. IF|8,60; thew >> Kees 3 
«Legs = = fenaichedt 0's ia TH eee 


terminal in’ 1(q¢-1,m) uader VT" is Ievaluod: ander ¥F ‘Also, “atty Tevahded teri @#'Lim) 


matched under Myr is matched to a "0". Id F(q4+1,4)“and' is entered to’ the sme "0" by Why 
Therefore, there must: be an unmatched: I-vatoed ‘termina: unider’ mii de Tq 1m). We can extend 
Myr so that the unmatched "0" in T(1,q) matches this "1". Then, 

Imatched "0"s in TU betes! t= mabe “O° sbagurta and 

Myr + Myp Myr + Myy - : 

We now suppose that |S |>6. We know that 
‘(matched "I's in BQH under yy. [ = batched "864 me emt ‘Mya - IS i 
(ZEROS yy (16-1): {ZEROS (18-1) + $6, Therefore: 
(ZEROB yy (ie th - jmiaathed:"1"s in: BGLs-19'tnidet- yy: f= 
_ AZEROB fhe 1 + |S)  jatohed "rsin Beet) under mg + B22 

There must be an unmatched “0” in B(1,s-1) under Mp since no"0" in aed) matches:a "tin 
B(s.n) under this matching function. “When my was defined on 4a) bi ff, it was'teft unmatched. . 
Modify VT,, and VB. so that VT (P(h)): = VB, th)'= I: ‘Ou: AEROS aim, myp is 
unchanged. On T(1,q), Lemma 5.8-A now applicd and’ we van cee aati sl valued a , 
on T(1,9) arc matched to 1-valued termisialé on TU,@). Pee 


Imatched "O"s in T(1,q) under my, | = : Lkqd as before. 
mn 
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On ZEROS yy (Leb. Miyp. is catended s0 that a prévioudy aanimatched "0". is matched to h. On 


ZEROS yy (80) Miyy isunchaaged. Therefore 
Myr +: Mvp. Pa Wrangelmyy M+: Wraripede yg Yee feange(m, + {range(my,)| = My, + Myp- 
Note that we have changed VT. and VB, so that p(h) and h are 1-valued. . oO 


Now that we have Theorem S4A, we naow how to “handle “aires hed op By 
Lemma 5: 1, there. are vale function vt, and VB, ‘consent “aid | - VI, and a V0, for which 
JONES (Cy = =¢,+¢,1 and which ahioe the seombealie os My, $ Muy pe all functions consistent 
with p, V1, and VB,. Given these functions, we apply bemirca: 5 AS to eae that wel are extensions 
of VT, and VB, consistent with p and satisfying (i) through (iv) of Theorem 5.4- -A which achieve 
My, + My, _ Therefore, for alf tetttials which are Valued of 1-Valued under VT, and VB... ¥ we can 
fix ie valves to be those under VT, ndVB, for RI terminall Whose’ values “are dictated by (ii 
through (iv) of Theorem 5.4-A, we can fix their values as dictated. We now have new initial value 
functions under which no terminal.in 'T(1,q) or B(s,n)} is ?-valued. We are guaranteed that there are value 
functions which achieve Mvr, + Mvp, ‘a which are consistent with these new initial value functions, 
Some terminals within Tq + 1,m) and B(1.s-1) may have acquired 0/1 values under the new initial value 
functions We en recursively apply ihe algorithm, beginning with, computation of, Icft-regions and 
right-regions, to find extensions of the new initial value functions which maximize the sum of matches on 
T(Lm) and B(1.n). The number of ?-valued terminals, bas, decreased, since at least h and p(h) have 
changed from ?-valued to O-valued. : 7 
When c,=1, we first need a modified version of Theorem SAA, ‘When ieee Lemma 5.7 only 
guarantees that there are functions cuisine the sais tai matching, Mv, + Mya, with IONES ,; (Ch 
equal to c, orc, +1. Theorem S4-A can only be appiied if JONES; (CH = Therefore, we prove a 


modified version, Theorem 5.4-B, for the case that (ONES, (CH =¢+1 


TEE RAIN SS 


IM: 


Theorem 54-B: Assume the standasd state after SCAN-ASSIGN retamns "FAIL". Let ¢=1. Let VT 
and VB be valuc functions consistent with VTy, VBy and p fer which-no ‘terminal is valued and soch 
that (ONESy,(C)| =.c/+1. There are extensions, VI, and VB,, of VI, and VB,, consistent with-p 
and such that: 

(i) no terminals are ?-valued, 

(ii) al 2-valued terminals | in T al 3 under VT, with pairs in rt #1) are 0-valued under Whey 

(ii) all ?-valued terminals in Tis) u under VB, with pairs in 1 Tat m) a are I-valued une 


VB, 


[- 4 


(iv)My. +Mo, >M.. + Moo. 
(iv) 4 a VT VB 


Proof: We nced to modify the proof to Theorem 5.4-A so that the roles of intervals T(1,q) and B(s,n) are 
interchanged. Functions VT,, and. VB,, are defined as. before except that we begin with VT, (o(h)) = 


VB,,(h) = 


1. In B(1,s-1), as for Theorem 5.4-A. 

If. In T(q+1.m), as for Theorem 5.4-A. 

Il. In B(sn). Lemma 58-B does not hold, but FEROS, (s,n) = ZEROS yy, ae “Therefore, 

OK-0(VB,,.y.n) holds for aly in Bis.) ‘Al L4(0-8+ 1d “0"s in B(.n) can be matched to "I's in sa) 

under VB,,. Asin the proof of Theorem $4-A: | ee Ne oo 
TONES, g(60) + Sy) = JONES, (ny) and 

matched "1"s on B(s,n) under ™\a,. |= JONES), (sn) - (F (ns 1) FL 4(n-s +1)4) 


> [matched "1"s on B(s,n) under "a + Bi ru ‘(n-s+ VL Kast 14) 


IV. Lemma 5.8-A does hold and we can define my, _s0 that p(h) is left unmatched and 


PETRME oe 6 + Sagde i thet Late See bce ide we es eet 


WS - 
{matched "0"s on T(1,q) under my | = IZEROS,, (1a) Z 
= ZEROS) (1.9) + [S,1 2, Jmatched "0's in Té.q) under my + 1S 
We have Myr + Myp P34 Myr + My, - (n-8+1)1-L A(n-s+ 1)4). If n-s+1 is even or 
any of the inequalities in I-chrough IV are strict, we have the desirediresult. Assume not. If {Sl =6, then 
L'A(n-s+ Dd +1 =4matched "1"s on Assn) under Mypl- 
There must be a 0-valued terminal in B(1,s-1) which matches-4 1-vahved terminal in R(s,n) under. Myp- 
This 0-valued: terminal under VB ig-also 0-valued under VB. and unmatched under Mvp We can 
modify Mp, that the unmatghed."1" in. B(s.n) matches this. "0", giving ha 
| Myp + Myp & Myr + Myp 
If {Sy} > 0, we modify VT., and VB... We have: 
JONES); (q+ 1,m)| - [matched "0"s in T(q+ 1,m) under Myr | 
= JONES,,(q+ 1m) + IS: (matched "0"s in T(q+.1,m)-under my] - So) 2 2 
We conclude that there is an. unmatched 1-valued terminal in T(q-+ 1,m) under m re Let VT. (p(h)) = 
VB,.(h) = 0. Now p(h) can match the unmatched I-valued terminal. in ‘1(q+1,m) and franglo, Wi 
increased by 1. In B(s,n), Lemma 5.8-B now applies and we still have 
|matched "1"s in B(s,n) under Myp 4 =L&(n-s+))J - 


Weconclude that My, + Myg | 2 Myr. Myg- a 


Given Theorem 5.4-B in addition to Theorem 5.4-A, we can dediies that there are value 
frit comiaeat wih VT, and VB, and satisfying (ii) and (iii) of Theorem $.4-A/B which achieve 
the maximum of Myr + Myp over all functions consistent with p. VT, and VB, However, we do not 
know what the value of p(h) and h should be. If the terminal pair is assigned "0", there is an extra "O” at 
the top which may be matched to some terminal in T(q+1,m), but one O-valued scrritial in B(s,n) will be 
unmatched. If the pair is assigned "1", the opposite is true. Therefore, we do not know which to choose 


until we have completed assignments within T(q+1,m) and R(1,s-1). Since p(h) and h may be the only 
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?-valued terminals in T(1,q) and 1(s,n) under. VT, and 'VB,, we cannot usea recursive application of the 
algorithm. Therefore, the algorithm must try both ‘possibilities for h andip(h). If we were to apply the 
algorithm recursively for each choice, the running time could be exponential in the number of terminals. 
Therefore, we use a modification of the procedures we have described se far: The main saeune: which 
given VT, and VB, finds VT-and VB which thaximize My. +'M,,,, is talled MAX-MATCH. The 
modified version is cated BETTER-MATCH. » 

When we are trying value “\" for h and pth), we will recursivefy apply MAX-MATCH, 
beginning with the computation of left-regions and right-fégions: on ‘FeLm) and: B(1.n) for initial vafue 
functions which extend VT, and VB, as dictated by. properties.(if) and Gi) of Theorem 5.4-A/B and 
under which p(h) and h are 1-valued. When we are trying value "0" for h and pth), we use initial valued 
functions which are the same as for the first choice except that pth) und h are O-vafued: However, in this 
case, BEFTER-MATCH is uscd. Procedure BETTER‘MATCH does not look for a pair of value 
functions consistent with the new initial value functions and: which ‘maximizes Myr + My, over all 
consistent functiou pairs. Rather, the procedisre idotis for a pair of consistent functions which maximizes 
M,, + Myp over all consistent function pairs and achieves’ a’ better inattthing than any function pair 
consistent with the initial value functions when pfh).afid hh are‘1-valued: ‘Consider all function pairs 
which maximize the sum of matchings on T and B over all flinctions consistent with p, VT, and VB... ‘If 
under one of alee h and p(h) are pve: Gen. a RIRcOS pair achieving a better matching with p(h) 
and h O-valued does not exist. popking only or beter machings allows uw o bound the search pace of 
the algorithm so that exponential rumsiag 6 time is avoided, i beeee that we could equally, well have chosen 
to look for maximizing pehines capa the ~ choice ane Lead matchings under the "1" choice. 


i ceena'S'9 ia vecd ts bound the sean: 
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Lemma 5.9 (see Figure 5.8): Assume the standard state after SCAN-ASSIGN returns “FAIL”. Let c,=1. 
Let VTyo atid VByo be value functions consistent with VT,,. VB and (i) and Gi of 
Theorem 5.4-A/B under which pth) and hi are O-valued. i under some matching function an which 
achieves Mya there is an unmatched "0" j in Bs. then there are valued functions VT,,; and VB,., 
consistent with all of the above and under which p(h) and h are I-vatued such that 


Proof: Let VT, and VB, mice with — and vi, everywhere cxcept at pth) and h, where 
VT,,;(@(h)) = VB, ,(h) = 1. Define mr to assign exactly those matches which ase assigned by 
mT je cone involve p{h), where MT — My; VT 9 Ben | 
My i= Jrange(eny | J > kangen Mri =Myy 
Since VB, ,, is an extension of VB, and VB, _.(h)=0, Leia 5.8-B applies. Let mB 0 be a 
matching function for VB.0 which agrees with MB for each range value in Bla }) but matches each 
"1" in BGs,n) to a "0" in Bis,n) while leaving h unmatehed— Any vO" in (Let) which ‘was unmatched 
under my, a is unmatched under m'yp ee cis 
|range(m’‘ VB.6 )}| = |matched "1's on B(1,s-1) under my ol + JONES yp on 
= lrange(myy 1= Mya 
Let ME be a matching function for VB 11 Which assigns all matches that mB assigns, and, in 
addition, matches h to an unmatched 0-valued terminal in BLS). Then 
Mvp > lrange(my, = lrango('yy +1 = Mya ot! and 


.M +M M +M oO 
VT ent VB ex} 2 VT x0 VB a0 


Lemma 5.9 allows us to modify the assignment procedures described so far while pursuing the 
0-valued choice for p(h) and h. While pursuing this choice, suppose that MAX-MATCH would make an 


assignment giving a new pair of value functions, cither after computing Icft-regions and right-regions or 
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Figure 5.8: Configuration in proof of Lemma 5.9. 


If: 
matched to something 


VT. .o(p(h)) = ee 
Tj “ye 


p(h) ¢ q m 


o_o 


B: 
cy VB,,o(h) = 0 
unmatched "0" RK unmatched 
Then becomes: 
unmatched unmatched 

Vea) = TQ ny ‘ oC aa 
T: 2 2 

l "9" S h n 

eo _?—______0>-"-—_____- 


B: 
So VB,, (hb) = 1 


matches 


Sienna a 
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after processing a return of "FAIL" from acall of SCAN-ASSIGN. Suppose we can deduce that there is 


a pair of value functions, VT, and VB consistent with pairing function p and the new value functions, © 


which maximizes My-+My, over all such consistent functions but for which a matching function 
achieving Mya, leaves a "0" in B(S-1) unmatched, Then BETTER-MATCH seein sia siiy 
pair of functions consistent with the new value functions, Nape will igus better matchings than the 
funetions found when pursuing the I-valued chile fir p(t) and h, This true because Lemma 5.9 can be 
applicd to VT, and VBj, to give functions vielding ‘ geod’a‘sum of matchings under which h and p(h) 
are:1-valued. 


If the pair of new value functions which:-has been rejected is one:of two choices after a call on 


-SCAN-ASSIGN which refuraed "FAIL", then this choice is hae BEYTERMA TCH never needs 


to pursue two choices. If the pair of new vive functions is defined by. assignment to left-regions and 
Tight-regions, then by Theorem 5.2, any pair of value functions. maximizing MyptMyp over value 
functions consistent with p and the new value functions also ain My;+ My, over value functions 
consistent with p and the old value antics Since shel mich a function under which a "0" in K(1,s-1) 
can be left unmatched, no function consistent with tie oe valued fancton will produce a sum of 
matings better re that siieed under the. I-valued choice. for h and p(h). Procedure 
BETTER-MATCH need not pursue functions consistent, with the old: valve functions either. Therefore, 
BETTER-MATCH returns "(null,null)", indicating that’ the O-valued choice for h jad pth) will not yield 
a larger sum of matchings than the 1-valued choice. For the same reason, BETTER-MATCH returns 
"(null,null)” apt when SCAN-ASSIGN returns *FAIL*, © 

The main procedure, saint “MA‘ TCH, is presente ie: Figure §.9. Our original routing problem 
is solved by calling MAX-MA‘ TCHCTEL.m) BA, mhp.v FVD) where WT, and VB, represent the local 


conncctions. The modified pee, BETTER-MATCH; used whist processing a 0-valued choice for 


p(h) and h, is ‘presented in Figure 5.10. Suppose BETTER-MA’ PCH(TC,m),B(1,n),p, VT, VBp) returns 
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- = 5.9: The main procedure. 


MAX- MENCRD (dl, m), BL, n)p.VT oVB,) 
/Finds extensions of VT and VB, which maximize. te Myy+Myg ore all such extenstions/ 


VE: =VT, and ap = VB, 
Compute 1 vp Lyp Ryp andRy, 
DO WHILE there are say e-valued terminals in T( Nes ), TRypmsE 1(1,L-y,)..and BeRygeh) 
FOR cach ?-valued terminal x in Td lyk DO 
; VI(@): =O and VBipiay=0 
END 
Compute R,,,, and Ry, /assigning “0"s only affects a aa cael 
FOR cach o-valued ‘nels xin Typ DO 
Vi(x):= 1 and VBpG=1 
END . . 
Compute Ly, andy Jassigning' ee affects feft-region/ 
FOR cach ?-valued termina y in BQ, Lyp)DO 
Vity): =O: and Vite =0 . 
END 
Compute Ryp and Ryp 
FOR cach ?-valued terminal yin MRyyah po. 
ViKy): =1 and Vipodi= 


— Compute Ly, and Lvs 
END 
IF there are any ?-valued terminals in Thy RyyDa and met wt 1 yy DT HEN DY 
. om: = Ryrlyrl. 
n':= Rvp “Lyp -] 


VT is such thet VFI) = VT@,,**) 
VB’ is such that VB’ = VB(Ly_ +x) 3 : 
~ p’ maintains pairs cofiesponding | to those under | p whien Leek oii of the p pair are 
in 'T(1.m') and BAI n'y. For other terminals in TEl.m’), 201), p’ assigns ‘"*" 
(STAT Us, vt p ¥B):= = SCAN-A ia mov 
IF STATUS = FAIL k" THEN BO... - 
VT:= vi, and VB’: = VB, 
Caleulate.q, &:h, andicy. eet tk 
IF p'(h)#k THEN VTi): =1 aad veg’ (oy. 21 
FOR ally ia '1(1,g) other that p{h).and DO" - 
IF owe ? THEN 
_ JF p(x) is ie A081) THEN: ~ 
VT"'(@): =6 and VB(p Wh: =O 
ELSE VI'@): =I and VB(p'):=1 . 
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FOR all y in B(s,n’} BO j 
IF VB(y)=? and p‘(y) is in oe m .) THEN 
VB(y): &F and VP¢p'(g)): = 1 
END 
IF c,> 1 THEN DO ar 
VB‘h): = and VE{p'fh)): <0 
(VT, VB): = ent alti Balan )p.VT.VB)) 
END se 
IFc, = 1 THEN DO 
VBh): = 1 and VT'(p(h)): = 1 
(VT,,,VB,):= Seeing Bln dp VT.VB) 
VB‘th): =0 and VI'(p! 
(VT 9 VB.) = BEET eat TCHCTC a 1, Ln’ ‘sp! V1I",VB)) 
IF VI,=nult OR Myr *Mva,, < Mir, i*Mya,, THEN 


(Vt..VB):=(VT,, VB.) 
ELSE(VT,, VB,):= =(TeVBQ) 
END 
END eee 
FOR each x in T(Ly,-+ LRy_-1) DO 
VQ)= VE Orb ys) 
_ END 
FOR each y in KL. vpt |. Ry 1) DO 
we VB (x-Lyp) 
END 
END 
RETURN(VT,VB) 
END MAX-MATCH 
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Figure 5.10: The modified pecciare: for O-valued — 
BETTER-MATCH(T(1,m),B(1,0}4p.VTy¥B,) 


/If can guarantee that a pair of extensions of VT, and’VB, which maximize My,+M,, over all such 
pairs of extensions allow an unmatched @vahied totuhinatst in: B(],d) under a pe matching, then 
returns (null,nuil}. . Otherwise, returhs a. pair Of extensions of VT, and VB, If d<n, it is assumed that 
B(d+1,n) has no ?-valued terminals under VB, and that matchings an aia and B(1,d) can ‘be 
maximized independent of each other./ 


VT: = VT, and VB: = VB, ° 
Compute Ly, Lyp, Ryp and Rog: fo , 
IF there is a right-region in BU d) containing can one terminal THEN 
RETURN(null nel) 
DO WHILE there are ‘any ?-vahued terminals j in Td mn) TRY), B(LLyg), and BR yn) 
FOR each ?-valued terminal i in TU Loyy) E DO, 
IF p(x) is in BR, d) THEN RETURN ie 
ELSE V1(x): =0 and seats =0 
END 
Compute Ry, and Ry, 
IF there is a ahiresion in 1(1.d) containing ati one : terminal THEN 
REFURN(aoull,null).__ 
FOR each ?-valued terminal x in Te BO 
VI(x): = 1 and VA(p(x)): = 
END . 
Compute Ly, and Ly, 
FOR each ?-valued terminal y in BULL yp) DO 
Vy): =0 and VI(pfy)): = 0 
END 
Compute Ry, and Ry, 
IF there is a right-region in 1(1,d) containing exactly one terminal THEN 
RETURN(null null) 
FOR cach ?-valued terminal y in BER, ,.R) DO 
Vy): = 1 and VI(ply)): =1 


END 
Compute L,, and Ly, 
END 
IF there are any ?-valued terminals in TMhyy+ LRy; 1) and iL yt Rya 1) THEN DO 
m:= RyrLypl 
n= Ry p eva -l 


VY’ is ich that VT'(x) = VT(Lyp+x) 

VB’ is such that VB'(x) = VKKLY +x) 

p’ maintains pairs corresponding to those under p when both terminals of the pair are 
in T(1,m’) and Kn’). For other terminals in T(1.m’) and B(1,n’), p’ assigns "*" 

(STATUS, VI, VB,):= SCAN-ASSIGN(T(Lm ),B(0'),p’, VI", VB) 
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. IF STATUS. =. "PAIL-4" THEN DO. 
VT’: = VT, and VB’: = bs 
Caleulate g, 4h, anda, 
IF c, > 1 THEN RETURN (ut 
‘Fo, = THEN DO ~ 
FOR all x in DO 
et oH g.00 
IF p‘(x) is in B(1,s- )) THEN 
VT"(x): =0 and VB(p' (x)): =0 
“BLSEVT'G):=1 and'VB(p'(x): =1 
__. __ fOnly & and p'(k) are. assigned by. the ELSE statement/ 
FOR all Yin Ban ‘I Da. 
“TF VBYyy=4 and p'(y) is in 'T(q+1,m') THEN 
a ; VAGhe] aed yt ‘(p'@)): =1 
END 


(VT,.VB,): = (BETTER-MATCH(T(,m’),B(,n'),n'p,VT.VB) 
IF VT. = nuff THEN: RETURN ( (null null) 
_. END. 
END 
_ FOR cach x in TL t+ hLRy DDO. 
bas = VEG Ty 
END 
FOR each y in (Ly, + LRyp- 1) DO 
VE) = VBCa-L yg) 
_ END 
END 
RETURN(VT,VB) 
END BETTER-MATCH 


A. 


VT, and VB, which are not "null". Functions vi; and eis pad not snaiooith Myr+ Myp ‘ over all value 
functions consistent with p, VT, and VBy. However, if My +Myp, beter that the maximum sum of 
matches for the 1-valued choice, then M vt Mv, does manimire Myy-+tMyp Lemma 5.10 states the 


property of Myr,+ Mvp, which allows us to make this conclusion. 


Lemma 5.10: Let VB, and d be such that if d¢<n, no terminal in B(d+1,n) is 1-valued, 
ZEROS, , (4+ 1.n} = L'%4(n-d)J+1, and for any O-valued tenninal j in B(d+1,n), each 1-valued 
terminal in B(d-+ 1,n) can be matched to a 0-valueid terminal in B(d+1,n) while leaving j unmatched. If 
BETTER-MATCH(T(L.m).3(1.n).4.p,VTVB,) returns "(nullnull)", then there is a pair of value 
functions consistent with P. VT, and VB, maximizing Myr+Myy over all such ieee functions and 
for which there is a matching function achieving the maximum matchinig’oil B which does not match all 
"O"s in B(1,d). If BETTER MATCHET(L 0) LW VTQ VEY peta sort? with VT, and VB, 
not “null”, then there are no ?-valued terminals under VE, and vB, and if = are VT and VB 

consistent with p, VT) and VB, for which MyrtMya> Myr, +My, then those: are VI, and VB, | 
_ consistent with p, VT, and VB, for which My, + Mya, 2 My;+Myag caer itader which diisdmaian 


function achieving M,, leaves a “O" in B(1,d) unmatched. 
8 


Proof: By induction on the number of recursive calls to BETTER-MATCH before the call being 
considered returns. | 

| Basis: no recursive calls. By Theorem 5.2, any extensions of VT and VB, defined in the 
WHILE loop, which assigns within full regions, can be extended to functions maximizing My; +Myp 
over all functions consistent with p. VT, and VB, We know tat id <a, d+ must be the end of a full 
fiiredioa: Otherwise, we would have: 


[ZEROS (d+ 1,n(d+ I <P e(r(d + 1)-d)1, which implies 


vay 


ONES, , (d +1,0(d+ 1) > L4(r(d + 1)-d).. 
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In this case, not alt "1"s in Bd +1,16d+1)) could be matched to “Os in B(d-+ 1,n), contradicting part of 
the hypothesis. It follows that € > Ryp,- Lat VB ané-VT be equal to VD, and VT, oF extensions of VB, 
and VV defined inthe WHILE loop. If under VB there is a right-region containing exactly one terminal, 
Bay j, then: 
ZFROS, (dM = ZEROS, J+ 1.d+1 2 PAC) 141 = Lad + DI+1 

For any extension of VB; there will be a matching function maximizing the matching on B'which matches 
all "1"s in D(d+1,s) to "0"s in B(d+1,n) and whict leaves a'“O" in BG.d) unmatched. Therefore, 
"(null,null)" is seturned correctly when: there is a right-tegion containing cxactly one terminal. 

If under some VT defined as above there: is a@ -valued terminal ia T(1,L,,,) with its pair in 
B(Ryp,4), then by Thoorem 5.2, there: are extensions VT, arid VB,,, of VB'and VT which assign this 
terminat the value "0" and maximize the sum of matchings on T and B over all extensions of VT and VB. 
Then: 

IZEROSyy Ryy@l AO MARy_ HIF 
There is a matching function achieving My, which matches all "Is ia B(d+1,n) to “0"s in. B(d-+ Ln) 
and leaves a "0" in B(Ryp.d) unmatched. Therefore, "(null.null)" is returned correctly. 

Let VT, and VB, be the final extensions of VF, and VB, upon-exiting the WHILE loop. Ifno 
terminals in T(Lyp +1.Ryr -1) and B(lyg +1.Ryp -1) are t-valuod, then VT, and VB, maximize 
M,‘+Myp, over all functions consistent with p, VT, and VB, There are no functions consistent with p, 
VT, and VB, for which My7-+My,>My; +Myp .and(VT,,.VB,) is correctly returned. If there are 
?-valued. terminals let T(1,m’) and B(L.n') be intervals TELy, +1Ryz -1) and: By, +1Ryp -) 
when renumbered. Let VT! and VB" be the restrictions of VT, did VB, to these intervals. Let VT, 
and VB_ be a pair of valued functions which maximizes the sum of matchings on T(1,m’) and B(1,n’) over: 
all functions consistent with p’, VI) and VB_. Theorem 5.1 guarantees that the function pair which 


agrees with VT, and VB, on Tat ARyy “l) sac Was LR yy and agrees with Vf, and VB, 
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elsewhere maximizes the sum of matchings og T and Raver all functions consistent.with p, VT, and 
VB,- Whea SCAN-ASSIGN, retures (SUCCESS,VT,VB,), VI,: and. VB,_ maximize: the sum of 
matchings on T(1m') and. B(Ln') and. do aot: teave: any. terminals“ ?-valued. ... Therefore, 
BETTER-MATCH correctly returns functions which maximize M,,+My, over all functions consistent 
with p, VIyand VBy Rr ee eet es eee | 
If SCAN-ASSIGN  reuums . (FAILANT,.VB,):-and c,d, thea’ -by Lomma’7. and 
Theorem 5.4-A, there is a terrainal s and » pair of functions; VT_ and VB,, which maximizes the sum of 
matchings on T(1,mn’) and.B(1,n') oves all Sunctions consiatent with pT and VB, sn: Sor whlch 
IZEROB,y (3.0) = LYla's+ db: : 
| Any ruatching function for VB, jeaves.a "0" in'B(sn').unmatched. The finctions.on T(1.m).and B(tn) 
which agree with VT, and VB, on T(yy +1Ryy. “lend MLyy +1.Ryy -D and'agree with VI, and 
VB,, elsewhere maximize M,,, +My, over all functions consistent with p, VT, and VB, Given these | 
functions, there is a maximum matching on: whic, ttches tach 5 in Rig .n) to a "0" in 
B(Ryp A) (invoking ‘Theory 5.1) and:teaves 2:°0" bn B+ Lm Bi 1G BUA) unmatched, ti 
case, "(null,null)” is returned corrostly. de 

- Induction Let ie mad NB, defined when c,=1 
upon the reuumn of SCAN:ASBIGNCT(L ey). B(La)4.N7, SVR cneptae pi) and h where VI,fo't)) 
= VA(h) = 2, Let VT and VB. bethe picnic Vand Vi, aking) nd th @-vahded anit 
VT, and VB,, be the extensions malting themt:drvalued, The yee fysttions woot Sor the recurve cal 
to BETFER-MATCH ate VT, aad VB,y. By Lemwna 5.7-and-Theoteans 6-8 and'5.4-B, any: fuaction 
pair maximizing My +Myy Ove all functions consistent with VE, end VB; macishiees M,;+ My, 
over all functions consistest with p’, VT and VEX, For any. function pair comsitoat with VT, and | 
VBo there are Lis(n'-s+1U+1 “0's in Bian’), where: sis defined under the standard state after 
SCAN-ASSIGN returns “FAIL”. One such “0” munt be unmatched underiatly matching: 
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Ef the recursive call-to: BETTER-MATCH using T(1,m’), B(Ln’), n’, p’, VT,, and VB_, returns 
“(null.aull)", then by inductive assumption, there is'a paix of valuod functions which maximizes the sum 
of matchings on Tm) ahd B(L.n’) over aff funeeiond consistent witt p’, VT_; and VB,, and for which a 
‘maximum: matching of BELn} feives 3 "0" in BELA’) wnivaitetied. Either this pair of functions also 
‘traximizos-My +My, over all fractionpaite consiatehit with p’ VT, and VB. or a function pair for 
which:p'th) and hare O-valued maximizes My+My, over these functions. In-ofther case, there is a 


value function pair which maximizes M,..+ My, over all functions consistent with p’, VI, and VB) and 


for which a "0" in Bf6.n') is unmaietied undet some maktinuit‘maiching for 1X(1,n'). As for the case when 
¢y> 1.in the proof of the basis above, "(null,nuth* we Uktabity-résurtiod by the original call to 
‘Suppose the recursive call to. BETTER-MATCH tetims VT_.and VB. which are not “null”. 
These functions are consistent with VT, and VB! wince’ VF," and ‘VB, afe consistent with VTi, and 
VB,. Given any pair of value factions on T(1,m’) and B(1,n’), let the expansions of these functions 
denote those functions on T(1.m) end: BEl.n).-wiich ‘agree’ with the given . functions on 
Typ +LRyp “land BiLyp +LRy, 1) and agree with VF, and VB, elsewhere. The original call 
to BETTER-MATCH returns VT, and WB, whieh are the expansions of VI. and VB.. Since no terminal 
is T-valued under VT, and VB., none is ?-valued under VT, and VB, Suppose there is a pair of functions, 
VT and VB, concioent with P, VT ‘iss VB, for whi Mut Myp bess He Mya, By Theorem $.2, 
there are functions VT, and VB, consistent with p, VT. and VB, for which Myr, +Myp, Pd 
My, +My,. Since pairs VT, VB, and'VT,VB, are both contistent with VT, and VB,, we can use 
Theorem 5.1 to deduce that the restrictions of these ‘functions to Tey +LRyr -D and 
BiLyg +1.Ryp -1) mast be such that ‘— 7 
Myr + Mya’ > Myr +Myp- 


where VT; and VB; are the restrictions of VT, and VB,. By Lemma 5.7 and Theorems 5.4-A and 5.4-B, 
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there are functions, VT, and VB., consistent with p’, VT, and VB, and assigning no "7"s such that: 
M t+M. f >M ‘tM a? r,s m8 be ee BG 
vi} Mya) & Mypi + Myg’? Myp My 

Case 1: Suppose VT, and VB, are consistent with VT,,.and VB,,..-For VB,, there. is an 
unmatched "0" under any matching function on HL.n’). Let VT, and VB, be the expansions of VT, and 
VB). Then Mor >t Mya, = My, +Myq.. and there isa, matnhing function achieving. My, which 
leaves a “0” in B(lyg +1.Ryg -D) GL) unmatched... Functions. VT and VBy afe'the desired VT, 
and VB. 

Case 2: Suppose VT, .and VB, arc consistent with VB, , _.. By inductive assumption, 
there are functions VT, and VB, consistent with. p' » VB, aided Myr tMypy > 
Myr) +Myp, and under which some matching function achieving My,’ Ieaves:a "O": it B(Ln) 

co. 
unmatched. Since VB, and VB, , are extensions of VT, and VB.,-VT; and VB, are consinent with VT! 
and VB,. Let VT, and VB, be the expansions of VT, and.VB,. They.are consistent with VIy and:VBy. 
We have: | oer dake = 
Mvp 2 My +M Mee oe: 
Myr + va, = Mv, va, 2 Myr + Myg - a , 
Under VB,, there is a matching function achicving Myg. which heaves “Of in. Mkyp +4Ryp “Ih a 


Lemma 5.10 ‘completes the technical erlopment needed ‘to ven the correctness of iad 


sigoridiim. Theorean'S’S states the cotiecinesd of the main procs MAX-MATCH. . 


Theorem $5: Let VT and VBy be ee ee the 
only ?-valued terminals. are members of top-bottom, pairs... ..Then MAX-MA'ECH(TELm)BiLada, 
VT» VB,) returns functions VT, and VB, consistent with p, VI, and:V By which maximize Mit Myp 


over all such consistent functions. 


Proof: The correctness follows from the development presented in this chapter. The formal argument ig - 
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by induction on the aumber of recursive calls; it is similar.te that used.in proving Lemma 5.10 except that 
verifying the correct return of "(null,nutl)” is not necessary. We only present the argument used when 
C)=1 after a call to SCAN-ASSIGN. it is ander ~~ — BETTER-MATCH is called by 
| MAX-MATCH and.Lemma 5. 10's nended i in the pent. 

Assume the sriodderi. sine after the call: to BCANSASSION in MAX-MATCH returns "FAIL". 
Let VT, and VB, be the extensions of the functions. weturmed xy SCAN-ASSIGN under which aff 
terminals in T(1,q) and Bés,n') except h and pi{h).are O-valued or I-valued. Let VF, and VB, extend 
VT, and VB, so.that p'(h) and:h. age O-valued, and VT. Tj AONB, ¢ extend them so that p'(h) and h are 
1-valued,. By inductive assumption, the functions Wry ead-¥B., eeturned by: MAX-MATCH(T(L,m’), 
B(1,n'),p', VT, ,,VB,,) maximize My, + M yy over all: fanctions consistant with p’, VT,, and VB,,. 

Suppose BETTER-MATCH(T(1.m’),B(1,0')s-1,p',VT, VB, 4) Fomurns “(oullnul)". Then 
there is a pair of functions VT 9 and VB..9 consistent with: VT, 9 and VB,, me maximizing Mae Mvp over 
all functions consistent with p’, Ving and: VByp for which there i. a:matching. function achieving Myy 
which does not match all "0"s:in B{is-1). By Lemma:S9, there are value functions VT; and VBL, 
consistent with p’, VT, , and VB, ; such that Myr +Myp 2 sass FO Mvp Then: 

Myr +Mvp z My +My 2 Myr tM ye 

and VT, and VB, maximize My,+My,: over. all functions consistent. with p', VT, and VB,. The 
expansions of VT, and VB_, maximize M, + M,,, over all flanctions consistent with p, VT, and VB,. 

Suppose BETTER-MATCH(T(1,m'),B(1.n'),s1:p',VT, 9 VB,,) returns VT, and VB, which are 
not “null”. Suppose VT,, and VB,» do not maximize M,--+My,, over-all fanctions consistent with p’, 
VT,» and VB.,. Let VIi.g and VB, ,.be value functions consistent with p’, VT,, and VB,, which do 
maximize My,+My, over all such functions. Then, Myr +My > Mvr +Mya By 
Lemma 5.10, there are functions VT and VB, consistent with p VT p and YB., for which Myr +My, 


> Myr ao 5 and for which a matching function. achieving. Myo, leaves a"O0" in Béls-1) 


VB... 
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unmatched. By Lemma 539, there are-value functions VT, and VB,, consistent with p’, VT,, and VB,, 


Myr Mya, >My +My = Myr +Myp “Therefore, 

Myr +Mya, 12 Myr Mv 2 Myr, +Myp My + Mya x 
Functions. VT,, and VB,, manimize My#M jy She a Ratio: cotalnte wieh p’ vq. and VB, and 
are correctly expanded and returned by MAX-MATCH. | 

_IEVT, and VB. 9 do maximize My + Mi, over al fonctions consistent with p’, VT, and VB. , 
‘then the maximum of Myr tMiyp ad Myp +My: machinists Mi, #M,,, over afl functions 
consistent with p’, VT, and NB.. pene of functions achieving the? maak iertunti (with ties broken in 


favor of VT, and VB,,) are correctly expanded and: retursied by MAX:>MATCH. ns 
5.6 Running Time of the Algorithm 


We will find an upper bound on the running time ofthe algerithm without making amumpttons 
about.the details of the: data structures-used to implement the alyoriditn. | We will assuine'that each of the 
following takes one 'step: addition to.a counter, testing the yatuecof & seales vatiable, atid nesignment wa 
scalar variable. hiaiguaiar toa funcGon variable ts asqurbed wo be aieignment to one scalar variable for 
cach domain value of the fuanction;.this allows a:Rinotion 4o-be.rebreseied a titray. “THe number of 


number of ?-valued. terminal paics unter he india alu etl bea i the — iter 


| Procedure SCAN-ASSION does not use afty' ae ee 


1 ‘This follows the unit cost model fot random acetim machines 6f Abo! fof, 
The logarithmic cost model, which sino ont te manbe oie nga a eps amb, 
would multiply any-runing time by tog,fei Fn fue 1. hist G47 Is _ 
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The procedure considers each terminal in T at most once, For each ?-valued terminal it considers, it must 

lest the sizes of ats of O-valued terminals and J-valued terminals and make an.apsignment when possible. 
Ifthe ses are recomputed for each ?-valued terminal conssesed, the. nynmber of steps used to process one 
t-valued terminal is af mast k(m+-n) for some constant k. The sopping test of.the first. FOR loop must 
also be computed, but only when: an assignment has been made. This test also requires computing. the 
size of a set of. terminals, Let SCmaua). denote the. maximum number of steps taken by 


SCAN-ASSIGN 02 any m (op terminals and.n bottom terminals with u.?-yalued terminal pairs when a, 


?-valued terminals pairs are assigned by SCAN-ASSIGN, Then 

SC(m,nu,a) $ k, (a+ 1m +n), for some constant k,. 
If a,¢u, SCAN-ASSIGN, returns "FAIL" and. (@,1)-8 She nyraber of ?-valued terminal pairs 
considered. Let SC(m,n.u) denote the maximum nuasber.of steps taken. by SCAN-ASSIGN on any m 
top terminals and n bottom terminals with w -valued secs pal. if 
SC(m,n,u) & maximum over all a,Su.of$Cim,n.ua,) $ k,fut m+ n).=. Olutm+n)). 

. _ Procedure BETTER-MATCH may use snicener SGAN-ASSIGN and a recussive call to itself. 
First BETTER-MATCH computes. the full JeA-regions and. cight-regions and. assigns to ?-valugd 
terminals within them. This is done in the WHILE loop. After assigning within aay of (1,L,,), (Ry pm), 
(1,L yp), and (Ryg.A), some regions. must be recomputed,, This takes, Ofm+ 8) steps. Once regions are 
computed, iesting,for the conditions under which "(null,qull)” is returned takes only a constant number 
of steps. Suppose a, initially ?-valued terminals are assigned.Q/} values in the WHILE loop. Then its 
execution, including the initial computation before gntering it, takes at most k,(a,+ 1m +n) stops, for 
soe constantly | eo ae 

After the WHILE loop is completed, if sy Sas pairs remain, SCAN-ASSIGN is called. 
Preparation for calling SCAN-ASSIGN takes O(m’+n’) steps to assign to VI’, VB’, and p’. If 


SCAN-ASSIGN returns "SUCCESS", then merging the functions that SCAN-ASSIGN returns with the 
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functions on the full regions to create the functions returned by BETTFR-MATCH takés O(n’ +n’) 
steps. If SCAN-ASSIGN returns "FAIL", it takes ‘O(m' +n’) steps to ‘assign to VT" and VB’ and calculate 
“4. h, and c. “If c)=1, BETTER-MATCH #s called tecurtively, retuming functions of “nutl*s. If the 
recursive call returns functions, then the processing required'to retuiin' their expansions takes O(m'+n') 
Let B(m;n,u) be the maximum number of stepstaten by any call w'BETTER-MATCH with m 
top terminats, n bottom terminals, and u initially t-vatied paits. Such a tats BETTPR-MATCH under 
which a, initially ?-valued pairs are. assigned in the WHILE loop “aed a, piairs are assigned: by 
SCAN-ASSIGN takes at most.’ : 
k,(a,+ 1Xm +n) + SCtm'n' waa) + Bom’,n'a-a, -a.<F) + kim‘tn’) steps: 
for some constant k,, where B(m,n,-1)=0 by definition for any'ni and n. Using the bound for sc gives: 
< ka, + m+n) + (4,41 en’) + Kom" n'y + Bin! ua 1) steps. 
We know 1<m'<m and l<n'<n. Fork, =k, +k,+«;, thie mimber OF steps is bourited above by 
(a, +0,+1Xm-+n) + Rim’ ni'u-a, aD). re ae : a f “expression 1 
Then Bm.n.u) < (maximum over at mS, n'Sn, and Sa, +4, Su ofeiprestion ny 
Order all tripfes (miniuh for m> i, n> to ahd w>0; leicicdgrapitcaty: “We prove by induction 
on this ordering that for afl audh tripks, Blm.n:u) k.{u+ Tain} for'somie constant k,. “The basis, 
(11,0), is trivial. For any triple, (21.1,0), B(m,n:0) & (vith), 2, bate Lt the Gropcetton be 
true for any triple fexicographicatty smaller that (mit,) where WO, We'tndw chat 1. <or-and 


I-¢n' $f whenever BETTER-MATCH ‘ts: afldd. Whei'a, +4, Cu, Osia,-a-t<ul! and 
we can use the inductive assumption. When atau, way-arl = “1 but k{ua,-a i489 = O® 
sit an upper Bound on the contribution of thé recive call ta BETTER-MATCH, ance no recursive 
call is made. Therefore: eC ee ae 


expression t < k,(a,+a.+ 1\m +n) + kfu-a, a Xm'+n) : 


aoe R gb a RCTS AE NR SR REISS RE ONS ME RGR 2 SOIL BRR TEE TF 


expression 1 ¢ k (a, +a,+1m+n) + k,(u-a,-a.(m +n) 
using ky > bande! +a! > m+n | 
$ a+ Dint+ny 
and D(m,n,u) <$ k(u+iim+n) = Om +n)). 

We now turn to MAX- MATCH. The bound on the number of steps used in the WHILE loop is 
of the same order as that for BETTER- maT by the same argument. Procedure MAX- “MATCH does 
not do the testing to return "(null,null)” which BETTER-MATCH does, but this only adds a constant per 
assignment to the number of steps used: by BETTER“MATCH, and onty affects the constant for 
MAX-MATCH'’s: bound. The pre-processing. for -8CAN+ASSIGN : and: the post-processing when 
SCAN-ASSIGN. returns "SUCCESS" have the same:bound on the mimber of steps used as that for 
BETTER-MATCH. However, if SCAN-ASSIGN: returns "FAIL" a récursive call to MAX-MATCH 
and possibly a call to BETTER-MATCH are used. ‘The pré-processing and post-processing for these calls 
take O(m'+n’') steps. Let M(m,n,u) be the:maximum number éf steps used by any call to MAX-MATCH 
with m top terminals; a er oee eae and u inisially ?-valued pairs. ‘Then fer constants k, and k, 

M(m,0,u) < maxiewe over all m'<ai,n'<n,and OSa, +a.<u of | 
k (a+ 1m +n) +S‘ n'u-a_a,) + B(nt,nura,-a-1) 
+ M(m'n’,0-8 -a+1) + tar ‘+ n’) expression 2 

We now prove by induction on the ordering of wiples that . | 

M(m,n,u) < k,(u+1)"(m-+n) for some constant ky. 
The basis and cases where u = 0 are straightforward if kgok, +k. i ae the proposition is true for all 
triples lexicographically less than (m,n,u) where u> 0. Then 
expression 2 < k (a, +1Xm+n) + kj(a.+1)m'+n') + k.(u-a,-a(m'+n’) 
+ k,(u-a,-a)*(m'+ n’) + k{m'+n') . 


<k,(u+1Xm+n) + k(u-a,-a (m+n) forky=k,+k, +k, +k, 
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Letting ky > ky, we have 
expression 2 <. k,(u? +ut+1{m+n) § £,a+d%em +n) 
and M(m,nu) < k,(u+1)(m+n) = O(u%\(m-+n)) =-Oeea+ap) ~ 
We have shown that our algorithm runs in a number of steps sae Saaeadl to the cube 


of the number of terminals. 


5.7 Summary 


In this chapter, we have presented. an algorithm which. routes the connections betwoen pairs of . 
terminals located on the outside of a robtenptlar cbiigibaeele ‘The rauting assumes horizontal and vertical 
wires are on separate laycrs and uses minimum. area. The problem jé reduces to an-sssignment problem 
on vectors. Procedure MAX-MATCH. solves’ the. probletr on.. vectors using- procedures 
BETTER-MATCH and SCAN-ASSIGN, - Procedurcs SCAN-ASSIGN and’ BETTER-MATCH have 
O((n-+n)’) running time and-MAX-MATCH has Ol(en+-n)") senning time.:. Procedure MAX-MATCH 
is actually used twice in routing terminals on a rectangie--once (cine bition cascactions and once for 
left-right connections. Before using MAX-MATCH, Jocal connections :ate: assigned directions in a 
predetermined manner. Assigning the lacal connections sases O(0) computation steps, where t is the 
number of terminals on the rectangle. ‘Each call to MAX-MATCH uses. O(t’) steps. Therefore, given a 
rectangular component with t terminals around: its boundary, the problem: presented in Section:$.1 can be 


solved in polynomial time or’). 


i Es Dera egret aa a a eee ee et EE Ce ER gE ae Se 
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Appendix to Chapter S: Notation | 
The notation used in Chapter 5 is listed here, with definitions, in order of appearance. : = 
T(1,m) denotes top terminals 1 through m. | 
B(1,n) denotes bottom terminals | through n. 


VT: T(1,m) —> {0,1,2} and VB: R(L.n) -» {0,1,7} indicae the directioas.of connections from top and 
bottom. terminals, respectively, as follows: - 


VT / VB @) = 0 if the direction of the path, from. terminal i to its pair is to. the left 
1 if it isto the right. 
? if it a undetermined 


p:TU.mUR(1n) TC,m)UBU,n)U{*} is the pairing function indicating what terminals should be 
connected: ; : vey ee Daeg aye aise tt dgtit Py Real suehetiay i ; 
For the following ere there are —" definitions for B and VB: 


Paeeroen 


my,,: T > T and is the function sathing terminals in TCL) of value "0" to sega bs with higher index 
in T(1,m) and of value "1", given value function VT. If ngs j, then Os 0, VTG) = 1, and i<j j: 


M,, is the maximum over all saiahing | hibclioes: oc for a als function VT, of aah 
ee = {i€S|VT(i)=0} for a value fuiseiin vt ‘ad 5 cr | 

ONES,,,(S) = {i€S|VT(i)= 1} for a value function vr sad s 1. 

UNDET (5) = {i€S|VT@)=} te value function VT and SCT. _ 

ONES, (x.y) is simplified notation for ONES,,(T(x,y)), otc. Pat BEE MSE TS 


OK-1(VT,x,y) (similarly OK-0) is true if and only if (ONES y4 (x.y) SLU(y-4-+1)J. 


Full-1(VT,x,y) (similarly Full-0) is true if and only if JONES,,(x,y)| > FA(y-x + 1) ae 
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lyy:T -> T and 1,;:T — T are used to define left-regions and right-regitiak ox T undet a ‘given’ vafue 
function VT: 
1,1) = 1 
Li) = i if Full-l(VT,1,,G-1)i-1) 
hyzli-1) otherwise, for i> 1 


Ty,{m) =m 
ry, 4) = j -if Ful-O(VT,i+ Lry i+ DY 
Ty {i+ 1) otherwise, for j <m 


Lefi-regions are defincd as the equivalence classes’ iftduccd by the equivalence relation on T(1,m) under 
which two terminals i and j are equivalent if and only if fit iso eh 


Right-regions are the equivalence classes induced ve the ecuilnce ras using vr 


Lyy i is defined to ey m if Full-1(VT, hyde), “i and to > equal wae otherwise. (1, ae are the-full. 
left-regions of T under VI. 


Ry is defined to equal 1 if Full-Q(VT,1 seftbond to — a otherwise. ym) are the full 
tight: epee eae va see 


Ne cities definitions are made after SCAN-ASSIGN Fetus = AL k, VT i 1 1 u 


ue 


q is the smallest i>k such that pide oa Li 


s is the smallest j satisfying the following Murer eee: 
OpwH2j : 
GD (ZEROS, Gm =L¥@it)s ees. 9 
(iii) each terminal in T(1,k-1) which has ‘been assigned the value “" e SCANASSION 1 is 
paired with a terminal in BG.n). wae og eT FD 2° 


Cis the set of initially ?-valued torminats in T(1,q) whose paiss arcim Bfga). 
¢ = IONESy, (Ch 


CO = ZEROS, (Oo. 
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h is defined as follows: Let i be the terminal of smallest index in B(s,n) which is ?-valued under VBp and 


L-valued under VB, _,. Certainly, p(i)<k. If i<p(k), then h = i; otherwise, h = p(k). 


VT, and VB, are defined as follows: Ifh = p(k), then VT, = VT;.; and VB, = VB, If h # p(k), 
then VI, and VB, agree with VT.) and VB. except at h, k, p(h), and p(k), where: 

VT,.(k) = 1 VT,(p(h)) = ? 

VB, (p(k)) = 1 VB,(h) = 2. 
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"Chapter 6 Discussion of the Algorithen 


&.1 Removing Ascomptions 


_ In this chapter we make some observations about the algorithm we have just presented. The 
algorithm is of the channel ‘nine variety. The streets to be'wsed by ach connection path: ure chosen: by 
the algorithm. Only one choice is considered for local connections For top-bottom and left-right 
connections, the procedure MAX- MATCH makes the decision. "Three assumptions are crucial to the | 
working of the algorithm: (1) only pairs of terminals need to be connected: (2) there are only four routing | 
areas -- one parallel to cach side of the rectangle; (3) the segments from each terminal to the routing area 
(perpendicular to the routing direction in this arca) cannot conflict regardless of their tengths. 

We do not necessarily need to have terminals around the outside of one rectangle as long as the 
above assumptions are satisfied. In fact, terminals may also lie along a rectangular boundary. 
circumscribing the rectangular component, i.c. on the opposite side of the routing area, as long as these 
terminals can be projected on the rectangle to obtain an instance of the original problem (sec Figure 6.1.). 
This configuration might be found when connecting terminals of a functional component or set of 
functional components to bonding pads. ‘The bonding pads He along the outside edges of the chip, and 
the routing regions lie berwten these sa and the rest of the Sei circuit. Mintoizing the area used 
for the interconnections minimizes the size of the chip. | 

| When the third assumption above is removed, we have alrcady. secn in Chapter 4 that the 
resulting routing problem is NP-complete even for one routing channel. with terminals ‘along its side. 
When either of assumptions (1) or (2) is sean our division into local connections and opposite side 
~ conncctions no longer limits the choice of directions for each connection. Consider the problem for one 
rectangular component when it may be necessary to interconnect three or more terminal When iiss 


terminals are involved, there are three choices instead of two for the type path ured to interconnect the 
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Figure 6.1: Alternate configuration of terminals for our algorithm. 


cannot be 
terminal here —>z|-— 


cannot be 
terminal here 


cannot be terminal here 


indicates a terminal 
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terminals. This can be seen by dividing the boundary ofthe rectangle into thre pecs ~ each pice pee 
ficken daisesat teraclaity wd ee Gomaday wu tavcsied: ‘Any path which flows two ajaceat piss of 
the boundary can be used to connect the tenninal. When the terminals are on thee different sides, two 
of the possible path typed use more sides that the third. However, we cagnot eliminate these possibilities 
ft is not true that the third choice is always as good as the other twa Figure 6.2-A gives a 
counterexampic. Comnepting these terminals involves three sides of the rectangle TRETOENs, we have 
atte independence ofthe two dimensions which we had when rong -tap tou ani fight 
connections. I tre terminals are oly on wo es, ten the pth ye wing ll ur ds cannot be 
climinated. This & shown in Figure 62-8, Its a. open problem whether thé one component routing 
problem can be solved in plynomic time when sc of ted or more rial ead 10 be 
interconnected. _ ats 

~ Consider retaining the restction that only paics of terminal ae interconnected but allowing 


terminals to ie on one of two rectangular components which are placed side by side Ignore for the 


moment any connection which involves terminals which Wie on the boundatis of the street between the a, 


two rectangles. The resuling rowing problem is mar to the one componént rating problem excep 

at pala can tae “shor cut” trough the set beoween the wo rectngis However, this option 
| dds many choles of paths o be wed to inerconnet terminal The top and bottom sides are no longer 
independent of the left and right sides. Figure 6.3 hares he dificubs presented bythe short cut 
The complexity of this routing problem is also open. 


62.4 Spacal Com for Pracetars MAX-MATCH 


The procedure MAX-MATCH, with its bub procedure, solves the ietsiiatiod probtom 
defined in Section 5.2. We now prove that wen all erminas ate T-rlued under the intial valu 
functions there is always an amignment of U's and "I" ch ta he se of maxi matchings 


et st SR ir a ee 
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Figure 6.2: The Ove Component Rotiag Problem with mar thant terminal net 
Terminals with the same number should be connected. — 
A: Three terminals on three sides -- three choices. 


1 
BEST OF 


THE THREE 


path uses 4 sides - ath isa 4 sides : 
routing addes 2 unit to height routing adds 3 units to height 
2 units to width 3 unit to width 


path uses 3 sides 


routing adds 3 units to height 
al ~ 2 units to width = 


B: Three terminals on two sides ~ three choices. 


1 
preferred whea kh < wl 


2} 


(3 + hXl + w) = 3+3w+h+hw = area 


(2 + hX2 + w) = 4+2w+2h+hw = area 


a IEA SERRA orig Teg ER oP 


er pie trot Baars on weal Sarg hein eet aR ee Sid gee aes fant ee ine hee ae Te 


-14@- 


Figure 6.3: Extending the routing problem of Chapter 5 to two components. a 


area = (2+w\h+3) i ee eae Beene 
=6+2h+3w+hw ae hee, 
area = (¥w2h) 
= 6+3h+2w+hw 
but: 


area = (24+w\(2+h) 


area = + wX3-+h) 
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swithin-one of the upper bound obtained when the epecific prising of top sermiaals to bottom terminals is 

disregarded. We always sumene: that. ail Tvalued.serminals.are- members: of toptbettom pairs. under p. 
When alt terminals are ‘thyalued wader the initial: value functions, the sumberef top terminals, m, must 
- qual the number of borwum terminals, n, 


Lemma 6.1 Let VT, and VB, assign "7"s to all terminals in T and B. | Then 


MAX-MATCH(T(1,m),8(1,m),p.VT,,VB,) returns VT, and VB, such tht Myr +My, 2 mL 


Reid 


_ Proof; Initially, there are no. full: regions and SCAN-ASSIGN is called. If SCAN-ASSIGN returns 
(SUCCESS, VT,,.VB), then Fm 7 -tenminale are O-valued. in each.of T-and-B.. We know OK-0(V 13,1) 
holds. Therefore, m is ever. ‘The matchings give Migr Mvp, = 4m +'am = m: 

_» Suppose .SCAN-ASSION -retoras (FAM RWT, VB). Let.g,.8 shy and:c, be as defined in 


Chapter 5. All terminals in T(k,m) are 7-valued under VI_, ‘Therefore, q. = kc, = 1, and kis odd. We 


‘know: 

JONES)» (1,k-1))'= JONES yy (s,m = b Mek = 4(k-1) 

f : m4 r 5 
1 - = ROS... a ina (Aiko Ly = A(k- zx f:} 
and | ages Das ZEROS yy (5.0 S: BEROS,, (ik = “a(-1) (a) 
Also, L'a(ras+ 1)J + Akl) = (ZEROS, am)}.+ JONES, (s,m) < {m-s+ 1)-1 
t ne: 9 

implying Mek l) SP Me(eas+ 17-1 (b). 


Combining (a) and (b) gives: 
L’A(mes+1)4 = Va(mrs): = 4(k-1) and (ZEROS), (.m)t = ZEROS, (1.k-D)} 
Hither VBj(6) = 1, or s = pth). Otherwise, JONESya (6+. 1.mik. © JONES, (san)h = ‘A(en-s)_ and 
B(s+ 1m) contains p(k) and.the pairs of all "1"s in T(1,q). ‘FRerefore, we would have chosen.s+1 rather 
than s. In cither of the two possible cases, s = h. 
Let VT,, and VB, , be the valuc functions defincd when the "I" values tried for h. We know — 


[ZEROSyq, (s+ 1m) = (mrs) and OK-O{VB, ,x.m} holds forall x in Beam). Therefore, s+ Lm) i 
x 


eae 


a set of full right-regions with no t-vahied serminsts.. All ewes) "0S in 1X34 1m) can be matched. 
Uniesss = lors = 2, intewat.B(J.s) contains one teft-region wid one sight-rcgion, neither of which are 
“full. Since {ONESy,. (1, ky = ("KT and OK-HVT,;.1 1:8) holds for afl’ in TULA-D), TUK) is a set of 
full right-regions in which %4(k-1) "O"s can be matched to “I‘s in“T(1:k), When's =2, in which ease 
k = m, we have | . | | 
| Myr +My_ = Mi) + Mom) = wl | | 
and the desired a Prey When s> 1, Tek+-Lm) c contains one ieft-region and one ae 
neither of which are full. If 5 = 2, B(1,2) 4s a full left-region andA€1) is valued under VB, |. ‘Tormitials 
- B(1) and p(BE1)) = ‘Hm) are assigned the value:”0” by the recursive call to MAX-MATCH. “Terminal 
BL) can match BQ), but ‘T(m)° is sniiaicheg.: “The: vesutting “inaximum matchings: sunt: to 
A(m-2)+ 14 Yetm-2) = me] as desired. {fs > 2, let Tm’) be the renumbering of Aik + Lm). Procedure 
SCAN-ASSIGN is called for T(m’) and B(1,s) using the appropriate restrictions of VT: and VB... 
If the second call to SCAN-ASSIGN returns (SUCCESS,VT.,.VB.,), then 
VEROSyq (13) = [ZEROS (1m) = PibenT and Mi! +My LAO Hem" = a 
Combining this with the matchings on die fuft regidns gives 
My +My = m’ + Vefars) + M(R1) = mek + bE = eb. 
where VT,, and VB, are the expansions of VT", and VB which agrce'wih VI), and VB, , on TX1,k) ated 
B(s+ 1m). | 
_ If the second calf to SCAN-ASSIGN returas{FAIL V1 VB), let q’, 8, and h’ be the special 
terminals for this call. ‘Under the initial value functions for this call ta SCAN-ASSION, ait porrminale in 
“Vim’yand all terminals in B(L,s-1) are valued and si I-valied: Thereford; oie’! is odd, anid 
; JONES, ,.48'8-1)| = ONES, A1.t-D} Mek or 
IZEROS,,,45'.8)| = Ltele-s'+1)4--¢ (ZEROS) ALES, Yeck!-2). 
~ Also, JONBS 43’, + ZEROS, (8. ae wee) +L FL Mees Sev +ty1. 
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Figure 6.4: The One Component Routing Problem with no local ee es 
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corey 
We 


HOST 0 rei OF? i giogd wiswenge grew vd obit ig Init gee sty band + aan © 


' 
iE nee . ae 08 aE! 3 Prod: rahe paltue ubqetie ytov 5 tot seta vowton qias98 BE 


hare. enol 388 fiw sie al ib ‘gilts MSY ianeS) WAM 


: vertical directions, we will 
ating dao tag balsas nog ti 
me ang wh pte re tt hg ave ch i 


roe 


od ano aprini ob tang Mit sa pation Bivens soothe ein od We. shay odds ia howasixe od ABs 


| = 148- | 
‘start later. The exchange causes the two originally connéctod segments to overlap but no other segments 
to overlap. (Soe Figure 6:5). tn this way,-we can modify the rotiting 96 that:there is oily ene horizontal 
segment for cach connection. “This sepmont must te above te terminator terminal to which it needs to 
be connected since some segment originally did, and segments have not been shortened, just merged. 
The vertical segments whiely conncet to servetitats ‘can be added of the vertical layer; Also, after the 
modification has beer made to each side, ‘segments which itist connect af the corners cin be lengthened 


or shortened to do'so without causing conflicts, © a Me i a | | 


this upper jan = the amount of epee added by using seek yer ‘te the two directions 
is actually achievable for a very simple routing problem shown in Figure 6.6. ee 

The argument used to prove Lenuma 6.2 tenis fora region above a'component side as‘long as 
ail wires ure cither perpendicelar: or parallel tothe side, segments of arbitrary feigth enterding om | 
different terminals perpendicular to the side will not intersect, and soginems whictr rit paralict to the side ; 
can be extended at the ends ae the side without oe pimatintend Big pamtaeke of eoniows © can be 


interconnected: not ony seein 


| Lt exo consider the‘ One Conipontnt Roun Problem wen tere are mare than tw ayers 


_ Vertical’ wire segments, ‘some: puraliet segments ieniat ton on i ater: "The coidctioof | 
eegments dn. different layers tno tonger trivial: ‘We opto ec ayery mi comets, no 


‘segment Intcrnceting: the connection point eat eon sy taper ec ee te coveted 


40 put sogments which conncet on adjatent ayers: aod abtiiseg ada shia soited a 
Sipe at wo wi na ep snipe br are nn 
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Figure 6.5: Construction for Lemenn 6.2. - 


A. Splitting positions. 


. de eit er. : - _ one channel, two layers 


becomes 
SPS EE LE TRINA RA ENT 


b. Merging segments for one connection. 
_ for layer connection. ees (pee " F 
a . . a for vertical connection 


after splitting: 


ightoft =< - -rightofl 
~ right of 2 ‘ 
after merging: 
~ right of 2 right of2 - - 
rightofl - - — “ 


— — rightofl 
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Figure 6.6: Instance of One Component Routing Problem realizing upper bound of Lemma 6.2. 
When one layer for horizontal segments and one layer for vertical segments; 


n channels used 


~{- nf2. channels used 
neven 


BP ting tse vanes see oo iets tompititte Dapeng he ein ge BE CRE ot ick abs 


Se mane: i ie PBs tas Sei sere eo 


than two layers. We will let oll dt smonnlpared! mee for herineatal: segments and all shen numbered 
Let “track” denote a RE ee as & ‘vas layer for each direction, the 
number of channels wed for routing along. she of te compat i the saci ax tie autbiet of tracks 
used. When there are H hpi forsegmen in tach diction, the number of channels required is as little 
as 1/H of the munber of traci required. More than ‘this number of channels may be needed to allow 
segments to consect at corner preety. ‘The connections at corners place constraints on the layers used. 

All layers in a channel may aot te pablo. It is important i) observe that there is no longer guaranteed to 
be an optimal routing in whieh each path vse at me one segment sunning ace sak side of the 
rectangular component. Permitting a path to use more than one ‘segment along a a side allows the path to 
change layer. In this way, layer may be beter wld, Figure 6.7 itustraces, / 

Let us suppose we- a iwiah te- meron ee routing tne segment ane each sida per path. We 

can choose which segments wit share the same track: eg. sing the shitching function of Chapter 5. 

Given any such choice, we must azvign the resulting scts of segments to channels anid-layers in the 
chanaels so that the connections at corkers are properly made. A legal assignment minkmizing the area of 
the layout is desired. We have bebe up the"muigmnen ipa Sears to tracks info two phases: 
| determining what segments will share the same track, arid- determining. what track they will share. “When 
there is only one layer for each direction, the second part of this problem is trivial -- any assignment will 
do. | 

If we require that all pairs of eoneecting nent by on macent layers, then we can model the 

constraints imposed by ~ connections hetivoga segments using a directed graph. The graph contains — 
one node for each set ot aan sharing a track. Let a node represonting a set of segments above the 
top of the component be called a tp node. Similarly, the terms bottom node, left node, and right node 


denote nodes representing scts of segments on the Indieated. sides. ‘There is an edge between two nodes 
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Figure 6,7; Layer assignment for the Onc Component Routing Problem. 


It is no longer truc that one channet per strect for each net is always optimal. 


Given 4 layers, no assignment of layers uses onty one channel on‘ cach side. 


One attempt is shown, a ian ed 


Graph representation of track connections for first track assignment 


o 2-right 
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when the sets of segments repmesemed by the nodes contain segments which must be connected at a 


comer. Edges sre directed as fellows: frum.a sop node to:# toft node, representing a connection in the 
top fet comer: fram 2 icf node toa battont note, sepresetting a connection in the bottom left corner; 
from a bottom node to a right mode, represénting  vonnection-in the bottom right corner; from a right 
node to. a top node, se Geiioes a-connection 1 the top tight curner. The ooh consists of disjoint 
chains and cycles. . Figure 6.7:shows the graph for the segments of that example. Since connections are 
always between. adjacent layers, i suffices to determine:the layer for cach set of segments. The actual 
channel used by each set of segments can be arbitrarily chosen, just as they could for only-two layers. The 
only. restriction is that there be only-one set of segmentsion each layer in any channel. ‘Therefore, we wish 
‘to assign a layer. number to each node of the graphy wo that left and right nodes have odd numbers and top 
and bottom nodes have even numbers. Adjacent nodes in the graph must be numbered with adjacent 
| numbers. For any given side, each oscurvence of a number will be in-a different charinel. We wish to 
find a numbering which minimizes'the area of the cortesponding layout | 
| - Let us now consider graphs consisting only of disjoint cycles. We restrict our attention to these 
graphs to illustrate how cyctic constraints can be handies without worrying ‘about graphs with different 
aumbers of tracks on different sides, which can occur if chains are present. A method of assigning layers 
w nodes in a saphiceasabig deb) of cycles is us follows. Choose some cycle: begin with some top node 
.of this cycle. Assign layer 1 to this node. Moving is beth directions away from this node, assign the 
nodes along the cycle -- one node in each direction at each seep:-- im the following pattern. Going in the 
direction of the edges (clockwisc), beginning with i=1, assign in consecutive steps: i+1 to the next left 
node; ito the next bottont node; i+ | to the next right node! i+ 2:to the next top node. Increase i by two 
and repeat the pattern. Going apposite to the dircction of the edges (counterclockwise), beginning with 
i= 1, in consecutive steps assign: 12-1 to the next right node; i to the next bottom node; i+ 1 to the next 


Jef node; +2 to the next top node. Increase # by. two.and fepeat the pattern. ‘Fhe first node reached by 


aM: 
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- both directions will cither be a top or a bottom node. Since these. are being ansigned the same mumber at 
"the same step, the cycle is cosapleted properly. if-the step assigning 2H to a right node clockwise and a 
left node counterclockwise is reached before the cycle is Gniaed. a pattern. of descending numbers is 
used. For clockwise, beginning with i=: 231-1, assign: i to the next top notie; i-] to the next left node;.i to 
“the next bottom node; i+] to the next right nade, Dectesne i by two arid sepast: For counterclockwise, 
beginning with i= 2H-1, assign:-i'to the next top, mode; +].t0 the neatiright node: i to the noxt bottom 
node; i] to the nest left node, Diecienes 3 by two: and repeat.’ Again, # tee sans alg Gs & right 
node clockwise and a Ieft node counterclockwise is-reached before the cycle is finished, the original 
‘pattern of ascending numbers is repeated. By. repeating. the: ‘twox patiesns through increasing and 
decroasing numbers alternately, a cycle of any siat canbe labeled, (Sor Figure 6.8}: ABcr the fir.eycle i 
finished, a second-cycle can be labeled-beginaing withia the next step.ofithe pattern which assigns to a top 
or bottom node. If the next step which assigns to: one of the two sides. assigns to atop node, then an 
arbitrary top node of the new cycle is chosea as the starting point. If the step assigns to:a bottom node, 
then an arbitrary bottom node is chosen .as-the: starting. point: Jir thin; way, cach cycke enn be assigned in 
turn. Each pass through inceeasing er decreasing labels, aces one channel oft each side.’:At-most one layer 
is. unused in any channel. ‘One pass through incressing; sumbes aasigns nustbers 10 1 groups of four 
nodes. Bach group consists of consecutive top, lei, bottom; auld tight nodesoe: a cycle. One pass through 
decreasing numbers assigns numbers to: Hr} much graupe.: *oesehoe she number DF channels uncom say 
tide under this aesigneeces mothot:ia:st mnt Sikes $} des . 
al Ranke) + en Seip heeud cae equation 1 
where nfcycle) is the sumber af groups of four neden im.the cycie! Note-dhat when a gsaph consiats of 
caly cycles, it containg the same number of andes for cach side. "Paeeofare equation tia cquatio. = 
» bal in, + the mumbes obcpeleey Bele 
"where a, isthe number of nodes for each side, The number nf chanecla-vand om any side at lea 17 
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~ the assignment technique to the the number of channels used on a side by an optimal assignment is at 
most 4H/2H-1 + 4H/n.. For2 < H< 1/12n,, this ratio is less a three. It is reasonable to ee n, to 
be much larger that H, since we do expect to have’ mang more interconnections, and therefore, segments, 
than layers for interconnect. eg or ee | | | 

The algorithen prescated in Chapter 5 minimizes the number of tracks used its routing ardund a 
rectangic. ‘The algorithm’ produces the sets of sckitienia which will share tracks, called 5 eins of the 
segmenis. However; she packing may not Iead to a minimum nee seit An ataienats packing of the 
segments may allow sabe layers to be used in some channels, resulting in bess channels being used. Even 
the distribution of tracks to the top and bottom or left and right may not be optimal.. ‘The algorithms for 
inp-Licktoon fliling ‘Goes ik tiy eb opliinice the Giewibudion of top tracks and bottotn tracks used. It only 
minimizes their sum. The “wrong” distribution might result in half-utilized cheapest In addition, our 


original choice of paths for local connections is no longer sufficient to find an ofitimal routing. Figure 6.9 


=e 


shows that the optimal routing may require that a local connestiga usa path ardund all four sides of the 


6.4 Summary 


The algorithm presented in Chapter 5 finds an optimal routing when certain restrictions are 
placed on the problem and on the allowed routing paths. In this chapter, we have discussed the 
repercussions of removing some of these restrictions. When we allow solutions which use only two layers 
but allow horizontal and vestical segments on both layers or which use more than two layers, the 
. sanction wo longer always fadi’en oplimel adie. In these case, the algorithm may be considered a 
heuristic algorithm, The algorithm limits the range of solutions it considers but finds an optimal solution 


within this limited set. When two layers are used for both horizontal and vertical scgments, we have 
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Figure 6.9: The One Compenent Routing Problem with more thin two layers. 


Connection for 5 goes around 4sides. 


ae 


2 py. adds 2 units to horizontal dimension 
a adds 2 units to vertical dimension 


Connection for 5 goes around 2 sides. 


adds 3 units to horizontal dimension 
| adds 3 units to vertical dimension 


——— horizontal layer 1 
+--+ horizontal layer 2 


| vertical layer 1 
; vertical layer 2 
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shown in Section 6.3 that the abaavoe Jengih added ta each dimension by the routing produced by the 
algorithm is never more that twice that added by the optimal routing. When separate layers are used for 
horizontal and vertical segments but there. , more, than one layer for each direction, the algorithm 
minimizes the number of tracks used in each direction. However, a new problem is encountered. Bach 
set of segments sharing a track mist be astigned to a channel and layer so that connections at corners can 
be indie sine sai eae of the resulting layout in minimized. Given sefs af segments to. be 
assigned, we do not know an algorithm to find an- optimal assignment, Atso, given a collection of 
segments, we do not know how to pack them into feos so that the jeecliiag area, rather than the 
resulting number of tracks, is minimized. Finally, given a One, Component Routing Problem, we do not 
know how to choose directions for connection paths so that the resulting afea, rather than the resulting 


number of tracks used, is minimized. 


In this thesis, we have examined problems encouetered in circuit layout from the perspective of 
complexity theory. Our goal has been to pain better understanding of the problems and discover better 


techniques for their sohation, We have presented ‘ome prof iW: ‘NP-completeness -- one for an 


orthogonal layout problem with rectangular component of arbicary size, one o for the channel assignment 
problem in a street, ‘and one ‘for the inncieection channel txignment siblen ie. channel assignment 


over all streets and intersections We also nave sbalyiad a common heuristic algorithm: for Nani 


i aa as This analy shows that the ratio of the scenber af Guuinels ined ti the | 


solution eoduiced by the algorithm to the optimal imnleet is not bounded by any condiaat: This 


algorithm and its analysis provide a reference point sasinat which other soda tec can be compared. The 


developmen of beuter sigorchms for channel csigament remains a ope for research. 
‘The proof of NP-completenem of the channel ssignment problen holds only when each net is 


required to use at most one chacinel, ie. channel ssignment ‘otthout joan The complexity of the channel 


signet problem with jogs remain pen The ator conjectures tat the problem is NP-complete. 

Tn practice, jogs are allowed within a strect. “Therefore, the ‘saniyels ‘of known channe! assignment 
lgoridis which uae Joga andthe development of beer algnitums are useful areas of research. The 
lower bounce on optimal channel aagnments without opt do nt apply when jogs are allowed, Analysis 
of sludons when jogs are slowed neemay to determined wef! lower bounds, 

In Chapter 5, we have presented an algorithm which finds the optimal channel routing for a 
special case of the layout problem in ‘ebanaied time. “Among other requirement, the layout protien 
must involve only one rectangular component, and all nets vast contain cxcly two terminals We have 
shown in Chapter 6 that when any restrictions on the problem are fenoved. i algorithm is no longer 
Giariainct to find the optimal solution. ‘Assumptions about satisfactory routing paths no longer hold. 


Wisse. wets abd ccwwesd’ta Wave eaiibe: than ic Getvainale bia al ollie reabicinek of Gages S hold, the 


REAR SAREE YES aI og 
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complexity of the layout problem in epee: A: modification <if the algoritfitir presented in Chapter 5 or a 
completly different algorithm may find an optimal solution in polynomial ime. Tie aibor ee been: 
unable to find such an algoritim or produce a prot of NP-hardness ‘An alternate relaxation ofthe 
restrislone aires ty exeponence whit he ine E7 A, but does not allow terminals on the adjacent 
sides. The rerction to two terminal nets i retained. ‘The abiity of paths wo run between the two 
rectangles is added This verion ofthe problem has so defied anata More generally, the achniqus 
used inthe algorithm may be ued as heuratisfor he general routing problem Further research is 
necessary to determine the quality ofthese techniques as heuriatics. | | 

hn itreting combinatorial problem i gpd by the rowing pole 0 of Chap 
Circular are coloring and routing around a rectangle resemble each other except tha, forthe rectangle, 
pats can go one of wo way around the rectangle and paths can change "colr” a ey net. ‘Suppers. 
Be eee A A Coens oe arcs” can go around he cen either direton More precinct, 
we modify he circular are coloring problem fellows Pair of point on acre ae given Tee peobiewt 
to choose one of he two ars determined by each pir and color the chosen act thatthe numberof 
colors used is minimized. Since crear arc coloring is NP-complete, we sumpect that this modification f 
Wehave chosen towel wich repens component a recangles ad wie  bjets 


wth uniform width and air interwire pacing. ‘The design rules for ategrated circuit design (M80) 


4ah4 YS FE ERS ON at 


are more complicated. Width and spacing ay fo ayer o ae. Under the model we me, the wide 


ofthe required widths and spacings must be we. ‘Ths sglifeion removes the completion of 
particular techndlagies fom the topological and geomen: specs of the layout problem. Probloma 


FEE tH TEs 


which vary sical fom our model, ach as over nd pound wie, equi pant siti, 
We assume that no topological information about the desired lout is given with the pritiom 
specification. This i in conan with the sick diagram approach reviewed in Chapter 2, which has 
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;. topological information as part of eon Both techniques.have becn found useful in practice, but lead 
to, differant probleme. The study of algorikhgrs fr, papblems,cogountered.in-alternate approaches to 
layout from that. in this thesis is an area ofcomplomentary rescatch 

Ta Chapter 3, ia addition to presenting the medel we have.used, we discussed a graph theorctic 
model for layout. This graph theosetic model, has proven very useful.in-obtaining bounds-on. the arca 
required by interconnection patterns represented by various classes of graphs... Thore are many open 
. problems. in this. arca of research. ‘Lwa,apen problemns.of great interest are to-determine tight upper and 
- lower bounds on the area nesded.to embed an arbitcary, planar graph and.on the.area needed to embed a 
shuffle exchange graph [Lei8Q], [Th80} An aljernate cost. measupe Jo-these prasented in- Chapter 3 has 
recently been propascd by Storer {StR0].- He, also proves the NP-bardness.of a embedding a 


‘gtaph in. the grid both under the equ CABUTES: BCH 


In the problems we. have. considered, placement and:couting are restricted.to be orthogonal. 
Resgarch is needed on the cost of this assumption. How much space is lost by allowing only horizontal 
and vertical wires? How much is gained by adding .a. third direction for. wires, e.g. 45 degree wires. 
Tompa [To80] has shown that far onc rset Comtaining sereainals: an.both sides suck, that all connections 
are of the form. "connect the if* temninal on one side to. the i!" terminal on the secend side", the optigal 
routing uses arcs of circies as wire paths. We have alsa gasumnad sat horizontal and-vertical. wires He en 
different layers. We have indicated how much this assumption costs for the problem of Chapter 5. How 
much this assumption costs in general is an opcn problem... When there are.only two: layers for 

interconnect, running two wires in. parallel, one on top of the. other,.on- separate laycrs across a chip 
separates the two sides of the chip. ‘This suggests that restricting each layer to wires in one direction is 
desirable. | 

At present, most integrated circuit designs are dose with only. two layers for interconnect. 


However, in the futurc, more layers may be available. When amore than two layers are. available 


restricting routing to wo orthogoual: directioins'inay the ‘ery coltly: in’ séitlie'o of atea of layouts: Should 


more directions be introduced and one direction: be used for each layer?’ Should the “different layer, 

different direction” technique be abandoned completely? For-more disi'two'layers, a ‘major area of 
research is to determine the restrictions which siiplify the routing problem sufficient! allow 
development of good algorithms and which résnit in good layouts: in retation to the layouts’ produced 
without the restrictions. Pye TEuKy 


_._s "Phe moidel we have used considers:only the aréa of layouts. In practice, other characteristics of 


circuit layouts, such as circuit ming and power consumption: ‘are 6 éfititerést.-‘Reseaith is needed into 


- ways of modeling these parameters. Very few algorithms hive ‘beet developed which take these | 


parameters into secount in any way. We would also ike to have a less ‘restrictive model of components. 
The extension of the mode? to include logically anid physically equivakent terminals was dlecused in ' 
Chapter 3. We would also like the ability to to-defime variable size componetite ~ ‘cotnponents which could 


sreich oF shrink depending: on the requirements of the tout. A very weft concept is that of - 
semi-restyicted regions for routing. These would i be shine ets aha WOR Bt A ee. whe 


n the thesis, we have € mparated placement and routing in the conventional Hii. The . 
‘interaction of placement-and ‘routing is very complicated. ‘Ott aten of tesearch is the dewopment of 
ses which do ot pee ponent ad rie with ae apa a ting 


an +16)- | 
| gecimetry of components adiita significant complication to the problem. When placement and routing 
‘are separated, % is non cloas-wivether it 4s/enprtant ty get's goo6 geometric: packing or a packing which 
facilitates routing. At one extreme, 4 plavement Which is a gdod packing of components as rectangles ai 
- bo-sought. This, of coune, ia NP-ocinplete, At the otter exthine; We placement can be governed by the 
desired interconncctions between componente with Seabee tr regan for how the rectangles fit together. 
oritseis for judging whether a placement does facilitate the routing are needed. The criteria currently 


used most often estimate wire length, but not routing area. 


We have intended our model to be independent of the particular technology of interest, 


although we have used the design rules of [Me86} for nMOS as our guide. These design rules are 
intended wo scale down as the fabrication technology becomes capable of creating srnaller and smaller 
objects. ‘Therefore, our model should remain -warthd-deupite the fict that actual paremeters of designs 

Ta sel nba, “Wo have: Ma <Goncernea. with. the Quality oF alabiiduna 40: acive. problenss 
encountered in layout. The criseria of complexity theory are not the only criteria useful in judging an 


algorithm to be used in practice. Algorithms should be easy to use and maintain. They should allow 


some interaction with the human designer. The algorithm presented in Chapter 5 finds optimal solutions. 
Therefore, interaction with the designer is not as important as it is for heuristic algorithms, where the 
designer's insight is an important input. Nevertheless, the algorithm of Chapter 5 can be written to allow 
the designer to predetermine the paths of some wires. It will find the optimal routing using these paths. 
Another criterion for judging algorithms is the number of good solutions which an algorithm can produce 
efficiently. An sigorim which can present a designer with several good choices wil be preferable wo one 


_ which presents one good solution, unless the one solution is substantially better than any of the choices. 


-_Bven among algorithms which find optimal solutions, one which obtains more than onc optimal solution 
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in approximately the same ruaaing time may be preferred.-Critesia otherthan the optimization criteria” 
may be used to choose among optimal. sohutions. For-exemple,. the algorithm of Chapter 5 is biased 
towards assigning "0"s to left top terminals, A. dual version of the. algorithm biased towards assigning 
"1"s to bottom right terminals might.also be run.to find an alternate solution. Complexity analysis can 
“assist in the development of better algorithms for layout problems, . However, the usability of an 
. algorithm must also be considered in evaluating new algorithnss for layout design. . 
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- heeeatie: Basic Defiaitions — 


nti the, the once of argh tory wil equal be ued to debe problems All 
definitions can be found in [Liu68}, We review them here. 

A.praph, G =: (¥.E) contin of a fae ent of: nedes for vertices). Vv, and a set of edges, E. Each 
edge, ¢ = (yyy) is a pair of nodes. Edge ¢ is said wo be incident on vertices v, and v5; vertices v, and v, 
are the endpoints of e. Vertices vad var alco 0 ach other. If the edges are ordered pairs, the 
"graph is directed otherwise, it ia wndirected. Ina directed graph, an edge (v,,v,) B08 from (or out of) 
wore: I nade ha edge incident it, then the node of dared 

A path ins raph a quence of get ek: — | 

ODMH WO ates S 


The nodes on the path are the mode vv av The path i of length n. Nodes v, and v, , , are the 


“at+L 
endpoints of the path; he path get rm, 1 If no node appears inthe sequence of edges more 
than once. the path i said to be simple. We wl aways mean “simple path hon we uy “peli”. A path 
is a cycle if v, =v, |). A simple eee is path on which only, = v, 1 appears more than once. A 
graph is acyclic if there are no cycles in the graph. An acyclic ndiecied graph Wicitlede renA 
subgraph, S = (V_Fy), of a graph, G = (V,,E,,), i a graph whose set of nodes, Ve, is a subset of V,, and 
whose set of edges, By, is mubeet of Fz, containing only edges whose endpoints ae ia V | 

A graph, G, is planar if there is a mapping from the nodes of G to points in the plane and from 
edges of G to curves in the plane such thet: . 

(1) The endpoints of the curve corresponding to an edge are the images of the endpoint of the 

edge; 

(2) No two curves intersect at any point other than their endpoints; 

(3) A curve does not intersect the image of a node unless that node is an endpoint of the edge 


corresponding to the curve. 
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Informally, a planar graph is oie whieh Gai bi dete in dhe pine without crossing edges. 
| _ An m by n two-dimensional grid graph is an undirected graph with node set: 
= {(ij) |i and j are positive integer, Asigm and 1SiSn} 
and edge set: | . 
B= {vy = “we ¥, = (ij+D) for 1gigm and 1<j<n-l 
= (ij) and = = at 1j) sok Isigml and Isisn}. 

To indicate the eke rate mot apne when diacussing bene pertarmaace of an algorithm, we 
need the following notation. Given sal auica functions f and'g on the same domain, f is O(g) if there is 
a positive constant, c, such that for al eu a finite ase of domain values, fx) s cat). The O notation 
is used when discussing upper bounds on a function. A statement of the form “the running time is 
Ole(n))" gives an upper bound on the running time. When discussing Souk bounds, the Q notation is 
used. Given real-valued bpacaes fand gon the same e domain, fis Oe) if there is a positive tetera CG, 
such that for all but a finite number of domain values, fs) 2 cpr). ines “the THAIN, time is O(a)" is 

| sang a lower bound on the running time. 

The above defo ae of the concept mos equenty used in his hes ea iit ead 


are prescnzed. as pagriaecal 


[Ab72] 
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